Will Leckie writes:
> When comparing the "precision" of the two quartic solvers, I find a
> very small (machine-level?) difference between the roots that they
> find, typically 10^-30, using gsl_complex_sub() to calculate the
> difference between roots.
Hello,
The error on the root should be
Hi,
Great work on the GSL, I love using it in my research.
I'm using the GSL in my research to solve arbitrary polynomials of order up to
4. I've written a little wrapper around the GSL to use
gsl_poly_complex_solve_quadratic() for quadratics and
gsl_poly_complex_solve_cubic() for cubics, and