请问二阶模态的rugge-kutta法的matlab格式是什么?