Hello, I just released an extension to gsl that computes jacobi polynomials, zeros of jacobi polynomials and several operations related to Gauss-Jacobi quadrature (numerical integration, numerical derivatives and interpolation using the quadrature nodes).
The extension can be found at: http://strouhal.mcef.ep.usp.br/pg/pjabardo/jacobi Now comes the question: I still have not implemented the error estimates. I tried to look at the code used to compute Legendre polynomials but I couldn't figure where did the error come from. Apparently it is related to the number of floating point operations. Does any one have any hints or references on that that I could use? Thanks for any help. Paulo Jabardo _______________________________________________________ Novidade no Yahoo! Mail: receba alertas de novas mensagens no seu celular. Registre seu aparelho agora! http://br.mobile.yahoo.com/mailalertas/ _______________________________________________ Help-gsl mailing list Help-gsl@gnu.org http://lists.gnu.org/mailman/listinfo/help-gsl