ode45求解微分方程

%% ode45求解微分方程
 %无附加参数
 %Lorenz 模型状态方程
F=@(t,x)[-8/3*x(1)+x(2)*x(3); %dx1
    -10*x(2)+10*x(3);
    -x(1)*x(2)+28*x(2)-x(3)];%匿名函数表示方程
tn=100;%终止时刻
x0=[0;0;1e-10];
%1e-10=10^(-10);
[t,x]=ode45(F,[0,tn],x0);

figure()
plot(t,x,'Linewidth',1.3)

figure()
plot3(x(:,1),x(:,2),x(:,3),'Linewidth',1.3)

figure()
comet3(x(:,1),x(:,2),x(:,3))%三维彗星图