求教一个给多项式赋值的问题
我想求一个挈比雪夫多项式的表达式,并且给表达式赋值,原程序没有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); 希望各位大虾们帮帮忙啊, 我自己已经把这个问题解决了是不是可以收回自己的这个悬赏,也收回自己的金币啊?
页:
[1]