lorfun.m

来自「四面体有限元程序,单元刚度的计算及总刚度矩征集成」· M 代码 · 共 11 行

function ydot=lorfun(t,y)
ydot=[-8/3*y(1)+y(2)*y(3);
     -10*y(2)+10*y(3);
     -y(2)*y(1)+35*y(2)-y(3)];


%dy1/dt=-8/3*y(1)+y(2)*y(3)
%dy2/dt=-10*y(2)+10*y(3)
%dy3/dt=-y(2)*y(1)+35*y(2)-y(3)