lorenzcomet.m

lorenz吸引子的画法

% lorenzcomet funtion
[t,y]=ode45('lorenz',[0 30],[12,9,7]);
clf;
axis([-20,20, -25, 25 ,10,50 ]);
view([20,42]);
hold on;
comet3(y(:,1),y(:,2),y(:,3));