matlab分数阶

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);