29.22.3 problem 609
Internal
problem
ID
[5203]
Book
:
Ordinary
differential
equations
and
their
solutions.
By
George
Moseley
Murphy.
1960
Section
:
Various
22
Problem
number
:
609
Date
solved
:
Sunday, March 30, 2025 at 06:51:58 AM
CAS
classification
:
[_separable]
\begin{align*} y \left (1+y\right ) y^{\prime }&=x \left (1+x \right ) \end{align*}
✓ Maple. Time used: 0.003 (sec). Leaf size: 494
ode:=y(x)*(1+y(x))*diff(y(x),x) = x*(1+x);
dsolve(ode,y(x), singsol=all);
\begin{align*}
y &= \frac {\left (-1+4 x^{3}+6 x^{2}+12 c_1 +2 \sqrt {4 x^{6}+12 x^{5}+24 c_1 \,x^{3}+9 x^{4}+36 c_1 \,x^{2}-2 x^{3}+36 c_1^{2}-3 x^{2}-6 c_1}\right )^{{1}/{3}}}{2}+\frac {1}{2 \left (-1+4 x^{3}+6 x^{2}+12 c_1 +2 \sqrt {4 x^{6}+12 x^{5}+24 c_1 \,x^{3}+9 x^{4}+36 c_1 \,x^{2}-2 x^{3}+36 c_1^{2}-3 x^{2}-6 c_1}\right )^{{1}/{3}}}-\frac {1}{2} \\
y &= -\frac {\left (1+i \sqrt {3}\right ) \left (-1+4 x^{3}+6 x^{2}+12 c_1 +2 \sqrt {\left (2 x^{3}+3 x^{2}+6 c_1 \right ) \left (2 x^{3}+3 x^{2}+6 c_1 -1\right )}\right )^{{2}/{3}}-i \sqrt {3}+2 \left (-1+4 x^{3}+6 x^{2}+12 c_1 +2 \sqrt {\left (2 x^{3}+3 x^{2}+6 c_1 \right ) \left (2 x^{3}+3 x^{2}+6 c_1 -1\right )}\right )^{{1}/{3}}+1}{4 \left (-1+4 x^{3}+6 x^{2}+12 c_1 +2 \sqrt {\left (2 x^{3}+3 x^{2}+6 c_1 \right ) \left (2 x^{3}+3 x^{2}+6 c_1 -1\right )}\right )^{{1}/{3}}} \\
y &= \frac {\left (i \sqrt {3}-1\right ) \left (-1+4 x^{3}+6 x^{2}+12 c_1 +2 \sqrt {\left (2 x^{3}+3 x^{2}+6 c_1 \right ) \left (2 x^{3}+3 x^{2}+6 c_1 -1\right )}\right )^{{2}/{3}}-i \sqrt {3}-2 \left (-1+4 x^{3}+6 x^{2}+12 c_1 +2 \sqrt {\left (2 x^{3}+3 x^{2}+6 c_1 \right ) \left (2 x^{3}+3 x^{2}+6 c_1 -1\right )}\right )^{{1}/{3}}-1}{4 \left (-1+4 x^{3}+6 x^{2}+12 c_1 +2 \sqrt {\left (2 x^{3}+3 x^{2}+6 c_1 \right ) \left (2 x^{3}+3 x^{2}+6 c_1 -1\right )}\right )^{{1}/{3}}} \\
\end{align*}
✓ Mathematica. Time used: 4.352 (sec). Leaf size: 346
ode=y[x]*(1+y[x])*D[y[x],x]==x*(1+x);
ic={};
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
\begin{align*}
y(x)\to \frac {1}{2} \left (\sqrt [3]{4 x^3+6 x^2+\sqrt {-1+\left (4 x^3+6 x^2-1+12 c_1\right ){}^2}-1+12 c_1}+\frac {1}{\sqrt [3]{4 x^3+6 x^2+\sqrt {-1+\left (4 x^3+6 x^2-1+12 c_1\right ){}^2}-1+12 c_1}}-1\right ) \\
y(x)\to \frac {1}{8} \left (2 i \left (\sqrt {3}+i\right ) \sqrt [3]{4 x^3+6 x^2+\sqrt {-1+\left (4 x^3+6 x^2-1+12 c_1\right ){}^2}-1+12 c_1}+\frac {-2-2 i \sqrt {3}}{\sqrt [3]{4 x^3+6 x^2+\sqrt {-1+\left (4 x^3+6 x^2-1+12 c_1\right ){}^2}-1+12 c_1}}-4\right ) \\
y(x)\to \frac {1}{8} \left (-2 \left (1+i \sqrt {3}\right ) \sqrt [3]{4 x^3+6 x^2+\sqrt {-1+\left (4 x^3+6 x^2-1+12 c_1\right ){}^2}-1+12 c_1}+\frac {2 i \left (\sqrt {3}+i\right )}{\sqrt [3]{4 x^3+6 x^2+\sqrt {-1+\left (4 x^3+6 x^2-1+12 c_1\right ){}^2}-1+12 c_1}}-4\right ) \\
\end{align*}
✗ Sympy
from sympy import *
x = symbols("x")
y = Function("y")
ode = Eq(-x*(x + 1) + (y(x) + 1)*y(x)*Derivative(y(x), x),0)
ics = {}
dsolve(ode,func=y(x),ics=ics)
ZeroDivisionError : polynomial division