Lorenz系统仿真动态可视化

%   输入：
%      f1,f2,f3 = y'(t,y) as a string
%            y0 = initial condition
%         inter = interval
%             h = step size

 
% 调用函数求解Lorenz方程
L=LorenzRK('-10\*y1+10\*y2','-y1\*y3+28\*y1-y2','y1\*y2-(8/3)\*y3',\[0,50\],\[5,5,5\],.01)
 
% 绘制Lorenz解
plot3(L(:,2),L(:,3),L(:,4))
     
% 创建GIF：

 
L=LorenzRK('-10\*y1+10\*y2','-y1\*y3+28\*y1-y2','y1\*y2-(8/3)\*y3',\[0,50\],\[5,5,5\],.01)