massode.m

来自「matlab从入门到精通第6章书后源码对初学者很有帮助」· M 代码 · 共 8 行

function dydt=massode(t,y,m1,m2,L,g)
% 微分方程矩阵
dydt=[y(2);
    m2*L*y(6)^2*cos(y(5));
    y(4);
    m2*L*y(6)^2*sin(y(5))-(m1+m2)*g;
    y(6);
    -g*L*cos(y(5))];