I compiled the test-program at http://www.gnu.org/software/gsl/manual/html_node/Roots-of-Polynomials-Examples.htmlfor retrieving the roots of a polynomial. I compiled the program using " *g++ -g -Wall -I /usr/include/gsl -L /usr/lib -lgsl -lgslcblas main.cpp -o gsltest*". However the output is as follows:
z0 = -0.809016994374947673 +0.587785252292473359 z1 = -0.809016994374947673 -0.587785252292473359 z2 = +0.309016994374947396 +0.951056516295153753 z3 = +0.309016994374947396 -0.951056516295153753 z4 = +1.000000000000000222 +0.000000000000000000 ...which is only correct up to a precision of 15 decimals, while the example output is correct up to at least 18: z0 = -0.809016994374947451 +0.587785252292473137 z1 = -0.809016994374947451 -0.587785252292473137 z2 = +0.309016994374947451 +0.951056516295153642 z3 = +0.309016994374947451 -0.951056516295153642 z4 = +1.000000000000000000 +0.000000000000000000 I'm using g++ (Ubuntu 4.3.2-1ubuntu11) 4.3.2, libgsl0 v.1.9-3, libgsl0-dev v.1.11-2 and 32-bits Ubuntu with kernel 2.6.24-19-generic. What is my problem? Jostein A. Jacobsen _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
