rigid.m

《MATLAB 7.0编程基础》第3章 （源码实例）主要讲解MATLAB的数值计算

% 首先编写方程函数，名为rigid.m
function dy = rigid(t,y)
dy = zeros(3,1);   
% 生成3×1的矩阵
dy(1) = y(2) * y(3);
% 第一个元素是 
dy(2) = -y(1) * y(3);
% 第二个元素是 
dy(3) = -0.51 * y(1) * y(2);
% 第三个元素是