lor.m

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

axis([10 50 -50 50 -50 50]);
view(3)
hold on
title('Lonrenz Attractor')
options =odeset('OutputFcn','odephas3') ;
[t,y]=ode23('lorfun',[0,20],[0,0,eps],options);