上传者: 39331817
|
上传时间: 2022-01-10 09:59:09
|
文件大小: 2KB
|
文件类型: -
matlab分数解方程,预估校正算法h=0.05*pi;alpha=0.5;gem=1/gamma(0.5);
x(1)=0;
for k=1:101
t(k)=h*(k-1);
end
for i=1:100
f=0;
a(i+1)=h^alpha/((alpha+1)*alpha);
for j=1:i
a(j)=h^alpha/(alpha*(alpha+1))*((i-j+2)^(alpha+1)-2*(i-j+1)^(alpha+1)+(i-j)^(alpha+1));
ff(j)=sin(j);