lorfun.m

四面体有限元程序,单元刚度的计算及总刚度矩征集成

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)