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

%主函數
function main()

clear all;
clc;

[t,x1]=ode45(@ly21,[0 100],[0.01 0.01])

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

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

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

grid