lor.m

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

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);