%混沌對初始條件非常敏感，即初始條件的微小差別常常使軌道按指數形式分開（蝴蝶效應）

%主函數
function main()

clear all;
clc;

[t,x1]=ode45(@ly26,[0 1e+8],0.001)

%x11=x1(:,1); 
%x12=x1(:,2);

figure
hold on
plot(t,x1,'r-')

%plot(t,x11,'r-')
%plot(t,x12,'k-')

[t,x1]=ode45(@ly26,[0 1e+8],-0.001)

%x11=x1(:,1); 
%x12=x1(:,2);

hold on
plot(t,x1,'k-')

%plot(t,x11,'r-')
%plot(t,x12,'k-')

grid