massode.m

matlab从入门到精通第6章书后源码对初学者很有帮助

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