10.1.24 problem 24

Internal problem ID [1121]
Book : Elementary differential equations and boundary value problems, 10th ed., Boyce and DiPrima
Section : Section 2.1. Page 40
Problem number : 24
Date solved : Tuesday, March 04, 2025 at 12:09:44 PM
CAS classification : [_linear]

\begin{align*} \left (1+t \right ) y+t y^{\prime }&=2 t \,{\mathrm e}^{-t} \end{align*}

With initial conditions

\begin{align*} y \left (1\right )&=a \end{align*}

Maple. Time used: 0.339 (sec). Leaf size: 21
ode:=(t+1)*y(t)+t*diff(y(t),t) = 2*t/exp(t); 
ic:=y(1) = a; 
dsolve([ode,ic],y(t), singsol=all);
 
\[ y = \frac {\left (t^{2}+a \,{\mathrm e}-1\right ) {\mathrm e}^{-t}}{t} \]
Mathematica. Time used: 0.076 (sec). Leaf size: 22
ode=(1+t)*y[t]+t*D[y[t],t] == 2*t/Exp[t]; 
ic=y[1]==a; 
DSolve[{ode,ic},y[t],t,IncludeSingularSolutions->True]
 
\[ y(t)\to \frac {e^{-t} \left (e a+t^2-1\right )}{t} \]
Sympy. Time used: 0.278 (sec). Leaf size: 15
from sympy import * 
t = symbols("t") 
y = Function("y") 
ode = Eq(t*Derivative(y(t), t) - 2*t*exp(-t) + (t + 1)*y(t),0) 
ics = {y(1): a} 
dsolve(ode,func=y(t),ics=ics)
 
\[ y{\left (t \right )} = \left (t + \frac {e a - 1}{t}\right ) e^{- t} \]