32.2.21 problem 21

Internal problem ID [7738]
Book : Engineering Mathematics. By K. A. Stroud. 5th edition. Industrial press Inc. NY. 2001
Section : Program 24. First order differential equations. Further problems 24. page 1068
Problem number : 21
Date solved : Tuesday, September 30, 2025 at 05:04:13 PM
CAS classification : [_Bernoulli]

\begin{align*} y^{\prime }+y&=x y^{3} \end{align*}

✓ Maple. Time used: 0.003 (sec). Leaf size: 39

ode:=diff(y(x),x)+y(x) = x*y(x)^3; 
dsolve(ode,y(x), singsol=all);
 

\begin{align*} y &= -\frac {2}{\sqrt {2+4 \,{\mathrm e}^{2 x} c_1 +4 x}} \\ y &= \frac {2}{\sqrt {2+4 \,{\mathrm e}^{2 x} c_1 +4 x}} \\ \end{align*}

✓ Mathematica. Time used: 9.555 (sec). Leaf size: 78

ode=D[y[x],x]+y[x]==x*y[x]^3; 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 

\begin{align*} y(x)&\to -\frac {1}{\sqrt {e^{2 x} \left (-2 \int _1^xe^{-2 K[1]} K[1]dK[1]+c_1\right )}}\\ y(x)&\to \frac {1}{\sqrt {e^{2 x} \left (-2 \int _1^xe^{-2 K[1]} K[1]dK[1]+c_1\right )}}\\ y(x)&\to 0 \end{align*}

✓ Sympy. Time used: 0.481 (sec). Leaf size: 49

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

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