正交多项式

[zhangjin]

在【0,1】区间的Legendre 多项式怎么在matlab中编程？

[zhangjin]

>> syms n
x=-1:1;
yn=(x.^2-1)^n;
pn=1/(2^n*gamma(n+1))*diff(yn,x,n);
??? Error using ==> sym.diff at 43
Do not recognize argument number 1.
这个错在哪？应该怎么编