预估校正算法解分数阶统一系统

function [ xyz ] = ff（ x0y0z0N ）


x=[0];
y=[0];
z=[0];
x（1）=x0;
y（1）=y0;
z（1）=z0;
h=0.01;
alpha=0.95;
delta=0.9;

%a（jxy）=（25*delta+10）*（y（j）-x（j））
%b（jxyz）=（28-35*delta）*x（j）-x（j）*z（j）+（29*delta-1）*y（j）
%c（jxyz）=x（j）*y（j）-（delta+8）*z（j）/3
%w（jn）=（h^alpha）/alpha*（（n-j+1）^alpha-（n-j）^alpha）
%rs（jn）=n^（alpha+1）-（n-alpha）（n+1）^alpha
%rt（jn）=（n-j+2）^（alpha+1）+（n-j）^（alpha+1）-2*（n-j+1）^（alpha+1）

for n=1:N+1
    Aw=0;
    Bw=0;
    Cw=0;
    Ar=0;
    Br=0;
    Cr=0;
    for j=0:n-1
        w1= w（ jnhalpha ）;
        aw=w1*a（x（j+1）y（j+1）delta）;
        Aw=Aw+aw;
    end
    xp=x（1）+1/gamma（alpha）*Aw;
    for j=0:n-1
        w2= w（ jnhalpha ）;
        bw=w2*b（ x（j+1）y（j+1）z（j+1）delta ）;
        Bw=Bw+bw;
    end
    yp=y（1）+1/gamm