massode.asv

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)*sin(y(5))-(m1+m2)*g
    y(6)
    -g*L*cos(y(5))
];