5.5.3 problem 3

Internal problem ID [1508]
Book : Elementary differential equations and boundary value problems, 11th ed., Boyce, DiPrima, Meade
Section : Chapter 6.5, The Laplace Transform. Impulse functions. page 273
Problem number : 3
Date solved : Tuesday, September 30, 2025 at 04:35:01 AM
CAS classification : [[_2nd_order, _linear, _nonhomogeneous]]

\begin{align*} y^{\prime \prime }+3 y^{\prime }+2 y&=\delta \left (t -5\right )+\operatorname {Heaviside}\left (t -10\right ) \end{align*}

Using Laplace method With initial conditions

\begin{align*} y \left (0\right )&=0 \\ y^{\prime }\left (0\right )&={\frac {1}{2}} \\ \end{align*}

✓ Maple. Time used: 0.174 (sec). Leaf size: 70

ode:=diff(diff(y(t),t),t)+3*diff(y(t),t)+2*y(t) = Dirac(t-5)+Heaviside(t-10); 
ic:=[y(0) = 0, D(y)(0) = 1/2]; 
dsolve([ode,op(ic)],y(t),method='laplace');
 

\[ y = \frac {{\mathrm e}^{-t}}{2}-\frac {{\mathrm e}^{-2 t}}{2}+\frac {\operatorname {Heaviside}\left (t -10\right ) {\mathrm e}^{20-2 t}}{2}-\operatorname {Heaviside}\left (t -10\right ) {\mathrm e}^{-t +10}+\frac {\operatorname {Heaviside}\left (t -10\right )}{2}+\operatorname {Heaviside}\left (t -5\right ) {\mathrm e}^{-t +5}-\operatorname {Heaviside}\left (t -5\right ) {\mathrm e}^{10-2 t} \]

✓ Mathematica. Time used: 0.129 (sec). Leaf size: 71

ode=D[y[t],{t,2}]+3*D[y[t],t]+2*y[t]==DiracDelta[t-5]+UnitStep[t-10]; 
ic={y[0]==0,Derivative[1][y][0] ==1/2}; 
DSolve[{ode,ic},y[t],t,IncludeSingularSolutions->True]
 

\begin{align*} y(t)&\to \frac {1}{2} e^{-2 t} \left (2 e^5 \left (e^t-e^5\right ) \theta (t-5)+\left (e^{10}-e^t\right )^2 (-\theta (10-t))+e^t+e^{2 t}-2 e^{t+10}+e^{20}-1\right ) \end{align*}

✓ Sympy. Time used: 0.966 (sec). Leaf size: 99

from sympy import * 
t = symbols("t") 
y = Function("y") 
ode = Eq(-Dirac(t - 5) + 2*y(t) - Heaviside(t - 10) + 3*Derivative(y(t), t) + Derivative(y(t), (t, 2)),0) 
ics = {y(0): 0, Subs(Derivative(y(t), t), t, 0): 1/2} 
dsolve(ode,func=y(t),ics=ics)
 

\[ y{\left (t \right )} = \left (\left (- \int \left (\operatorname {Dirac}{\left (t - 5 \right )} + \theta \left (t - 10\right )\right ) e^{2 t}\, dt + \int \limits ^{0} \operatorname {Dirac}{\left (t - 5 \right )} e^{2 t}\, dt + \int \limits ^{0} e^{2 t} \theta \left (t - 10\right )\, dt - \frac {1}{2}\right ) e^{- t} + \int \left (\operatorname {Dirac}{\left (t - 5 \right )} + \theta \left (t - 10\right )\right ) e^{t}\, dt - \int \limits ^{0} \operatorname {Dirac}{\left (t - 5 \right )} e^{t}\, dt - \int \limits ^{0} e^{t} \theta \left (t - 10\right )\, dt + \frac {1}{2}\right ) e^{- t} \]