global?x?y?z
a=1;?%parameter?a
q1=0.9;q2=0.9;q3=0.9;?%the?fractional?order
x0=3;y0=1;z0=5;
h=0.01;N=1000;
x（N+1）=[0];y（N+1）=[0];z（N+1）=[0];
x1（N+1）=[0];y1（N+1）=[0];z1（N+1）=[0];
M1=0;M2=0;M3=0;N1=0;N2=0;N3=0;
?
x1（1）=x0+h^q1*（25*a+10）*（y0-x0）/（q1*gamma（q1））;
y1（1）=y0+h^q2*（（28-35*a-z0）*x0+（29*a-1）*y0）/（q2*gamma（q2））;
z1（1）=z0+h^q3*（x0*y0-z0*（8+a）/3）/（q3*gamma（q3））;
?
x1=x0+h^q1*（25*a+10）*（y1（1）-x1（1）+q1*（y0-x0））/gamma（q1+2）;
y1=y0+h^q2*（（28-35*a-z1（1））*x1（1）+（29*a-1）*y1（1）+q2*（（28-35*a-z0）*x0+（29*a-1）*y0））/gamma（q2+2）;
z1=z0+h^q3*（x1（1）*y1（1）-z1（1）*（8+a）/3+q3*（x0*y0-z0*（8+a）/3））/gamma（q3+2）;
?
for?n=1:N
????M1=（n^（q1+1）-（n-q1）*（n+1）^q1）*（y0-x0）;
????M2=（n^（q2+1）-（n-q2）*（n+1）^q2）*（（28-35*a-z0）*x0+（29*a-1）*y0）;
????M3=（n^（q3+1）-（n-q3）*（n+1）^q3）*（x0*y0-z0*（8+a）/3）;
????
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1814??2014-10-30?11:15??上傳\BIANLIANG.m
?????文件?????????826??2014-02-24?12:31??上傳\chen_lyapunov.m
?????文件?????????480??2014-10-30?11:43??上傳\HANSHU.m
?????文件?????????370??2014-11-06?09:06??上傳\solvechen.m
?????文件????????2496??2014-09-27?18:36??上傳\Unified.m