求教一个给多项式赋值的问题

xuexihaohao321

我想求一个挈比雪夫多项式的表达式，并且给表达式赋值，原程序没有error，但是运行了好几天还是一直busy，我觉得可能是陷入了死循环等原因，但是又找不出原因，希望各位高手和大师指点小弟一下，小弟不胜感激。
源程序为：
clc;
clear;
%以下为求递推公式及其导数
syms x
T0(1)=sym(1/sqrt(pi));
T0(2)=(2*x/0.21-1)*sqrt(2/pi);
for i=3:21
    T0(i)=2*(2*x/0.21-1)*T0(i-1)-T0(i-2);
    T1(i)=diff(T0(i),x,1);
    T2(i)=diff(T0(i),x,2);
end
T1(1)=0;
T1(2)=(2/0.21)*sqrt(2/pi);
for i=3:n
    T1(i)=subs(T1(i),x,i/100);
    T2(i)=subs(T2(i),x,i/100);
end
T111=double(T1);
T222=double(T2);

xuexihaohao321

希望各位大虾们帮帮忙啊，

xuexihaohao321

我自己已经把这个问题解决了是不是可以收回自己的这个悬赏，也收回自己的金币啊？