lorenzcomet.m

来自「lorenz吸引子的画法」· M 代码 · 共 7 行

% 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));