function?adams（x1y1z1hn）
a=35;b=3;c=28;
[xyz]=crk（x1y1z1h3）;
for?i=4:n
fxi_0=a*（y（i）-x（i））;
fxi_1=a*（y（i-1）-x（i-1））;
fxi_2=a*（y（i-2）-x（i-2））;
fxi_3=a*（y（i-3）-x（i-3））;

fyi_0=（c-a）*x（i）-x（i）*z（i）+c*y（i）;
fyi_1=（c-a）*x（i-1）-x（i-1）*z（i-1）+c*y（i-1）;
fyi_2=（c-a）*x（i-2）-x（i-2）*z（i-2）+c*y（i-2）;
fyi_3=（c-a）*x（i-3）-x（i-3）*z（i-3）+c*y（i-3）;

fzi_0=x（i）*y（i）-b*z（i）;
fzi_1=x（i-1）*y（i-1）-b*z（i-1）;
fzi_2=x（i-2）*y（i-2）-b*z（i-2）;
fzi_3=x（i-3）*y（i-3）-b*z（i-3）;

x（i+1）=x（i）+h/24*（55*fxi_0-59*fxi_1+37*fxi_2-9*fxi_3）;
y（i+1）=y（i）+h/24*（55*fyi_0-59*fyi_1+37*fyi_2-9*fyi_3）;
z（i+1）=z（i）+h/24*（55*fzi_0-59*fzi_1+37*fzi_2-9*fzi_3）;

fxi=a*（y（i+1）-x（i+1））;
fyi=（c-a）*x（i+1）-x（i+1）*z（i+1）+c*y（i+1）;
fzi=x（i+1）*y（i+1）-b*z（i+1）;

x（i+1）=x（i）+h/24*（9*fxi+19*fxi_0-5*fxi_1+fxi_2）;
y（i+1）=y（i）+h/24*（9*fyi+19*fyi_0-5*fyi_1+fyi_2）;
z（i+1）=z（i）+h/24*（9*fzi+19*fzi_0-5*fzi_1+fzi_2）;
end

figure（1）;
plot3（xyz）;
title（‘Chen等式的三維相空間吸引子‘）;
%plot（x）;
grid;
figure（2）;
plot（xy）;
title（‘Chen等式的xy平面投影圖‘）;
%plot（y）;
grid;
figure（3）;
plot（xz）;
title（‘Chen等式的xz平面投影圖‘）;
%plot（z）;
grid;
figure（4）;
plot（yz）;
title（‘Chen等式的yz平面投影圖‘）;
grid;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1199??2012-07-11?20:14??adams\adams.m

?????文件????????793??2012-07-11?20:15??adams\crk.m

?????目錄??????????0??2012-07-11?20:14??adams

-----------?---------??----------?-----??----

?????????????????1992????????????????????3