90.19.13 problem 16

Internal problem ID [25295]
Book : Ordinary Differential Equations. By William Adkins and Mark G Davidson. Springer. NY. 2010. ISBN 978-1-4614-3617-1
Section : Chapter 4. Linear Constant Coefficient Differential Equations. Exercises at page 309
Problem number : 16
Date solved : Sunday, October 12, 2025 at 05:55:39 AM
CAS classification : system_of_ODEs

\begin{align*} \frac {d^{2}}{d t^{2}}y_{1} \left (t \right )+2 y_{1} \left (t \right )&=-3 y_{2} \left (t \right )\\ \frac {d^{2}}{d t^{2}}y_{2} \left (t \right )+2 \frac {d}{d t}y_{2} \left (t \right )-9 y_{2} \left (t \right )&=6 y_{1} \left (t \right ) \end{align*}

With initial conditions

\begin{align*} y_{1} \left (0\right )&=10 \\ D\left (y_{1} \right )\left (0\right )&=0 \\ y_{2} \left (0\right )&=10 \\ D\left (y_{2} \right )\left (0\right )&=0 \\ \end{align*}

✓ Maple. Time used: 0.047 (sec). Leaf size: 41

ode:=[diff(diff(y__1(t),t),t)+2*y__1(t) = -3*y__2(t), diff(diff(y__2(t),t),t)+2*diff(y__2(t),t)-9*y__2(t) = 6*y__1(t)]; 
ic:=[y__1(0) = 10, D(y__1)(0) = 0, y__2(0) = 10, D(y__2)(0) = 0]; 
dsolve([ode,op(ic)]);
 

\begin{align*} y_{1} \left (t \right ) &= -15-{\mathrm e}^{-4 t}+26 \,{\mathrm e}^{t}-30 \,{\mathrm e}^{t} t \\ y_{2} \left (t \right ) &= 6 \,{\mathrm e}^{-4 t}-6 \,{\mathrm e}^{t}+30 \,{\mathrm e}^{t} t +10 \\ \end{align*}

✓ Mathematica. Time used: 0.676 (sec). Leaf size: 510164

ode={D[y1[t],{t,2}]+D[y1[t],{t,1}]+2*y1[t]==-3*y2[t], D[y2[t],{t,2}]+2*D[y2[t],t]-9*y2[t]==6*y1[t]}; 
ic={y1[0]==10,Derivative[1][y1][0] ==0,y2[0]==10,Derivative[1][y2][0] ==0}; 
DSolve[{ode,ic},{y1[t],y2[t]},t,IncludeSingularSolutions->True]
 

Too large to display

✗ Sympy

from sympy import * 
t = symbols("t") 
y1 = Function("y1") 
y2 = Function("y2") 
ode=[Eq(2*y1(t) + 3*y2(t) + Derivative(y1(t), t) + Derivative(y1(t), (t, 2)),0),Eq(-6*y1(t) - 9*y2(t) + 2*Derivative(y2(t), t) + Derivative(y2(t), (t, 2)),0)] 
ics = {y1(0): 10, Subs(Derivative(y1(t), t), t, 0): 0, y2(0): 10, Subs(Derivative(y2(t), t), t, 0): 0} 
dsolve(ode,func=[y1(t),y2(t)],ics=ics)
 

Timed Out