29.19.32 problem 545

Internal problem ID [5141]
Book : Ordinary differential equations and their solutions. By George Moseley Murphy. 1960
Section : Various 19
Problem number : 545
Date solved : Tuesday, March 04, 2025 at 08:12:29 PM
CAS classification : [_rational, _Bernoulli]

\begin{align*} 2 x y y^{\prime }+x^{2} \left (a \,x^{3}+1\right )&=6 y^{2} \end{align*}

✓ Maple. Time used: 0.006 (sec). Leaf size: 45

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

\begin{align*} y \left (x \right ) &= -\frac {\sqrt {4 c_{1} x^{4}+4 a \,x^{3}+1}\, x}{2} \\ y \left (x \right ) &= \frac {\sqrt {4 c_{1} x^{4}+4 a \,x^{3}+1}\, x}{2} \\ \end{align*}

✓ Mathematica. Time used: 0.949 (sec). Leaf size: 59

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

\begin{align*} y(x)\to -\frac {1}{2} \sqrt {4 a x^5+4 c_1 x^6+x^2} \\ y(x)\to \frac {1}{2} \sqrt {4 a x^5+4 c_1 x^6+x^2} \\ \end{align*}

✓ Sympy. Time used: 0.573 (sec). Leaf size: 46

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

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