second-order differential equation

在用ode23函数求解动力学方程时出错，请问如何修正。

function dydt=mytest1(t,x)
dydt(1)=x(4);
dydt(2)=x(5);
dydt(3)=x(6);
exp1=vpa(inv(Mthea)*equction,6);

 dydt(4)=vpa(subs(exp1(1,1),[var1 var2 var3 dotvar1 dotvar2 dotvar3],[x(1) x(2) x(3) x(4) x(5) x(6)]),3);
 dydt(5)=vpa(subs(exp1(2,1),[var1 var2 var3 dotvar1 dotvar2 dotvar3],[x(1) x(2) x(3) x(4) x(5) x(6)]),3);
 dydt(6)=vpa(subs(exp1(3,1),[var1 var2 var3 dotvar1 dotvar2 dotvar3],[x(1) x(2) x(3) x(4) x(5) x(6)]),3);

这里的exp1是我另外求出的一个动力学公式。