工学解析って何？！　日本オンリーワン株式会社

ローレンツ方程式(その２）

前回に、余計な項が含まれていたので、解きなおした。

>function xdot=f(x,t)
xdot(1)=-10*x(1)+10*x(2)
xdot(2)=-x(1)*x(3)+28*x(1)-x(2)
xdot(3)=x(1)*x(2)-8/3*x(3)
endfunction
>x0=[1;1;1];
>t=linspace(0,500,1000);
>x=lsode("f",x0.t);
>plot(t,x);