6Â¥: Originally posted by zÂÜÀò at 2012-09-22 20:46:03
u=zeros(100,1);
a=zeros(100,1);
u(1)=1;
for i=1:100
a(i)=(-1)^(i-1)*gamma(i)/gamma(2.5-i);
end
for j=2:100
for i=1:j-1
u(j)=u(i)+a(i)*u(j-i);
end
end

ÔÚmatlabÀïÃæÊýÁеÄÏÂ±í²»ÄÜ´Ó0¿ªÊ¼ ...

2Â¥: Originally posted by dingd at 2012-09-22 09:36:09
r()±íʾʲôº¯Êý°¡£¿