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