連立二回微分方程式の解き方を教えてください

d^2y/dt = a11*dy/dt-a11*u*θ+(a12+u)*dθ/dt+b11*u

d^2θ/dt = a21*dy/dt-a21*u*θ+a22*dθ/dt+b21*u

a11,a12,a21,a22,b11,b21,uは定数です。この連立二回微分方程式の解き方教えていただきたいです。ode45で解こうとしたのですがエラーが出てしまいました。

[Akira Agata]

ode45 を使って解くには、1階の連立微分方程式に書き直す必要があります。今回のケースですと、 とおいて、 の 4変数に対する 1階の連立微分方程式に書き直したうえで ode45 を適用して下さい。