LegendrePoly

This program returns the coefficients of the Legendre polynomial P_n, given n.

Description

This program returns the coefficients of the Legendre polynomial P_n, given n. The result is a row vector with powers of x decreasing from left to right (standard MATLAB notation). Like for other polynomials, LegendrePoly(n) can be evaluated at x by typing

polyval(LegendrePoly(n),x).

• Acknowledgements: This submission has inspired the following:; LegendreShiftPoly
• MATLAB release: MATLAB 6.0 (R12)

[Joerg Enderlein]

 

[anas baha]

thank you very much

[Thijs Stegeman]

Hi,
I find your function interesting and usefull, but I'm afraid summing the vectors produced reveals an error; while the sum of the coefficients should yield 1 (as Pn(z)=1 for every n), in fact this goes wrong (at least on my computer) at n=25 because of roudof-error (I think), to deviate to something of the order of 1.0e20 at n=99. Do you know of any kind of solution?
Other than that it works very well!
Thanks.

[Liang G]

I think there is runoff error when n>39 in my computer.Can anyone fix this problem?