2.3.13 Problem 13

Internal problem ID [9070]
Book : First order enumerated odes
Section : section 3. First order odes solved using Laplace method
Problem number : 13
Date solved : Wednesday, March 05, 2025 at 07:19:00 AM
CAS classification : [_separable]

Solve

$$\begin{array}{rl}{y}^{\prime }+(at+bt)y& =0\end{array}$$

With initial conditions

$$\begin{array}{rl}y\left(0\right)& =0\end{array}$$

We will now apply Laplace transform to each term in the ode. Since this is time varying, the following Laplace transform property will be used

$$\begin{array}{rl}{t}^{n}f\left(t\right)& \stackrel{\mathcal{L}}{\to }(-1{)}^n\frac{d^n}{d{s}^n}F(s)\end{array}$$

Where in the above $F(s)$ is the laplace transform of $f\left(t\right)$. Applying the above property to each term of the ode gives

$$\begin{array}{rl}(at+bt)y& \stackrel{\mathcal{L}}{\to }-a\left(\frac{d}{ds}Y\left(s\right)\right)-b\left(\frac{d}{ds}Y\left(s\right)\right)\\ y^{\prime }& \stackrel{\mathcal{L}}{\to }Y\left(s\right)s-y\left(0\right)\end{array}$$

Collecting all the terms above, the ode in Laplace domain becomes

$$Ys-y\left(0\right)-a{Y}^{\prime }-b{Y}^{\prime }=0$$

Replacing $y\left(0\right)=0$ in the above results in

$$Ys-a{Y}^{\prime }-b{Y}^{\prime }=0$$

The above ode in Y(s) is now solved.

In canonical form a linear first order is

$$\begin{array}{rl}{Y}^{\prime }+q(s)Y& =p(s)\end{array}$$

Comparing the above to the given ode shows that

$$\begin{array}{rl}q(s)& =-\frac{s}{a+b}\\ p(s)& =0\end{array}$$

The integrating factor $\mu$ is

$$\begin{array}{rl}\mu & =e^{\int qds}\\ & ={\mathrm{e}}^{\int -\frac{s}{a+b}ds}\\ & ={\mathrm{e}}^{-\frac{s^2}{2a+2b}}\end{array}$$

The ode becomes

$$\begin{array}{rl}\frac{\mathrm{d}}{\mathrm{d}s}\mu Y& =0\\ \frac{\mathrm{d}}{\mathrm{d}s}\left(Y{\mathrm{e}}^{-\frac{s^2}{2a+2b}}\right)& =0\end{array}$$

Integrating gives

$$\begin{array}{rl}Y{\mathrm{e}}^{-\frac{s^2}{2a+2b}}& =\int 0ds+c_1\\ & =c_1\end{array}$$

Dividing throughout by the integrating factor ${\mathrm{e}}^{-\frac{s^2}{2a+2b}}$ gives the final solution

$$Y=c_1{\mathrm{e}}^{\frac{s^2}{2a+2b}}$$

Applying inverse Laplace transform on the above gives.

$$\begin{array}{}\text{(1)}& y=c_1{\mathcal{L}}^{-1}({\mathrm{e}}^{\frac{s^2}{2a+2b}},s,t)\end{array}$$

Substituting initial conditions $y\left(0\right)=0$ and $y^{\prime }\left(0\right)=0$ into the above solution Gives

$$0=c_1{\mathcal{L}}^{-1}({\mathrm{e}}^{\frac{s^2}{2a+2b}},s,t)$$

Solving for the constant $c_1$ from the above equation gives

$$\begin{array}{r}{c}_1=0\end{array}$$

Substituting the above back into the solution (1) gives

$$y=0$$

Figure 2.79: Solution $y=0$

✓ Maple. Time used: 1.399 (sec). Leaf size: 5

ode:=diff(y(t),t)+(a*t+b*t)*y(t) = 0; 
ic:=y(0) = 0; 
dsolve([ode,ic],y(t),method='laplace');
 

$$y=0$$

Maple trace

`Methods for first order ODEs: 
--- Trying classification methods --- 
trying a quadrature 
trying 1st order linear 
<- 1st order linear successful`
 

Maple step by step

$$\begin{array}{lll}& & \text{Let’s solve}\\ & & [y^{\prime }+(at+bt)y=0,y\left(0\right)=0]\\ \bullet & & \text{Highest derivative means the order of the ODE is}1\\ & & y^{\prime }\\ \bullet & & \text{Solve for the highest derivative}\\ & & y^{\prime }=-(at+bt)y\\ \bullet & & \text{Separate variables}\\ & & \frac{y^{\prime }}{y}=-at-bt\\ \bullet & & \text{Integrate both sides with respect to}t\\ & & \int \frac{y^{\prime }}{y}dt=\int (-at-bt)dt+\mathit{C}\mathit{1}\\ \bullet & & \text{Evaluate integral}\\ & & \ln \left(y\right)=-\frac{t^2(a+b)}{2}+\mathit{C}\mathit{1}\\ \bullet & & \text{Solve for}y\\ & & y={\mathrm{e}}^{-\frac{1}{2}{t}^{2}a-\frac{1}{2}{t}^{2}b+\mathit{C}\mathit{1}}\\ \bullet & & \text{Use initial condition}y\left(0\right)=0\\ & & 0={\mathrm{e}}^{\mathit{C}\mathit{1}}\\ \bullet & & \text{Solve for}\mathit{\text{\_C1}}\\ & & \mathit{C}\mathit{1}=\left(\right)\\ \bullet & & \text{Solution does not satisfy initial condition}\end{array}$$

✓ Mathematica. Time used: 0.002 (sec). Leaf size: 6

ode=D[y[t],t]+(a*t+b*t)*y[t]==0; 
ic=y[0]==0; 
DSolve[{ode,ic},y[t],t,IncludeSingularSolutions->True]
 

$$y(t)\to 0$$

✓ Sympy. Time used: 0.340 (sec). Leaf size: 3

from sympy import * 
t = symbols("t") 
a = symbols("a") 
b = symbols("b") 
y = Function("y") 
ode = Eq((a*t + b*t)*y(t) + Derivative(y(t), t),0) 
ics = {y(0): 0} 
dsolve(ode,func=y(t),ics=ics)
 

$$y\left(t\right)=0$$