MATLAB: Math Legendre problem

legendremathematicspolynomial

Initial conditions p(0)=1 and p(1)=x , we also have that p(n+1)=((2*n+1)*x*p(n)-n*p(n-1))/(n+1).
As we all know, matlab cannt start from p(0), therefor i do this :
p(1)=1, p(2)=x and p(n)=((2*n-1)*x*p(n-1)-(n-1)*p(n-2))/n
but i DON'T take the same results and i don't know the reason !! As you see i am working with Legendre Polynomials !!
Thanks !!

Best Answer

That's almost right, but wherever you had p(n) etc, you shouldn't adjust n. You should only have subtracted 1 from the instances of n that are on their own. So change all your p(??) bits back to their original values and it should work.
Related Question