29.9.15 problem 255

Internal problem ID [4855]
Book : Ordinary differential equations and their solutions. By George Moseley Murphy. 1960
Section : Various 9
Problem number : 255
Date solved : Sunday, March 30, 2025 at 04:04:24 AM
CAS classification : [_linear]

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

✓ Maple. Time used: 0.002 (sec). Leaf size: 27

ode:=x^2*diff(y(x),x)+x*(x+2)*y(x) = x*(1-exp(-2*x))-2; 
dsolve(ode,y(x), singsol=all);
 

\[ y = \frac {-3+{\mathrm e}^{-2 x}+x \,{\mathrm e}^{-2 x}+x +{\mathrm e}^{-x} c_1}{x^{2}} \]

✓ Mathematica. Time used: 0.177 (sec). Leaf size: 32

ode=x^2*D[y[x],x]+x*(2+x)*y[x]==x*(1-Exp[-2*x])-2; 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 

\[ y(x)\to \frac {e^{-2 x} \left (e^{2 x} (x-3)+x+c_1 e^x+1\right )}{x^2} \]

✓ Sympy. Time used: 0.425 (sec). Leaf size: 31

from sympy import * 
x = symbols("x") 
y = Function("y") 
ode = Eq(x**2*Derivative(y(x), x) - x*(1 - exp(-2*x)) + x*(x + 2)*y(x) + 2,0) 
ics = {} 
dsolve(ode,func=y(x),ics=ics)
 

\[ y{\left (x \right )} = \frac {\frac {C_{1} e^{- x}}{x} + 1 + e^{- 2 x} - \frac {3}{x} + \frac {e^{- 2 x}}{x}}{x} \]