3.27 problem 81

3.27.1 Solving as quadrature ode
3.27.2 Maple step by step solution

Internal problem ID [3345]
Internal file name [OUTPUT/2837_Sunday_June_05_2022_08_41_27_AM_65249478/index.tex]

Book: Ordinary differential equations and their solutions. By George Moseley Murphy. 1960
Section: Various 3
Problem number: 81.
ODE order: 1.
ODE degree: 1.

The type(s) of ODE detected by this program : "quadrature"

Maple gives the following as the ode type

[_quadrature]

\[ \boxed {y^{\prime }-\operatorname {a1} y-\operatorname {a2} y^{2}-\operatorname {a3} y^{3}=\operatorname {a0}} \]

3.27.1 Solving as quadrature ode

Integrating both sides gives \begin {align*} \int _{}^{y}\frac {1}{\textit {\_a}^{3} \operatorname {a3} +\textit {\_a}^{2} \operatorname {a2} +\textit {\_a} \operatorname {a1} +\operatorname {a0}}d \textit {\_a} = x +c_{1} \end {align*}

Summary

The solution(s) found are the following \begin{align*} \tag{1} \int _{}^{y}\frac {1}{\textit {\_a}^{3} \operatorname {a3} +\textit {\_a}^{2} \operatorname {a2} +\textit {\_a} \operatorname {a1} +\operatorname {a0}}d \textit {\_a} &= x +c_{1} \\ \end{align*}

Verification of solutions

\[ \int _{}^{y}\frac {1}{\textit {\_a}^{3} \operatorname {a3} +\textit {\_a}^{2} \operatorname {a2} +\textit {\_a} \operatorname {a1} +\operatorname {a0}}d \textit {\_a} = x +c_{1} \] Verified OK.

3.27.2 Maple step by step solution

\[ \begin {array}{lll} & {} & \textrm {Let's solve}\hspace {3pt} \\ {} & {} & y^{\prime }-\mathit {a1} y-\mathit {a2} y^{2}-\mathit {a3} y^{3}=\mathit {a0} \\ \bullet & {} & \textrm {Highest derivative means the order of the ODE is}\hspace {3pt} 1 \\ {} & {} & y^{\prime } \\ \bullet & {} & \textrm {Solve for the highest derivative}\hspace {3pt} \\ {} & {} & y^{\prime }=\mathit {a0} +\mathit {a1} y+\mathit {a2} y^{2}+\mathit {a3} y^{3} \\ \bullet & {} & \textrm {Separate variables}\hspace {3pt} \\ {} & {} & \frac {y^{\prime }}{\mathit {a0} +\mathit {a1} y+\mathit {a2} y^{2}+\mathit {a3} y^{3}}=1 \\ \bullet & {} & \textrm {Integrate both sides with respect to}\hspace {3pt} x \\ {} & {} & \int \frac {y^{\prime }}{\mathit {a0} +\mathit {a1} y+\mathit {a2} y^{2}+\mathit {a3} y^{3}}d x =\int 1d x +c_{1} \\ \bullet & {} & \textrm {Evaluate integral}\hspace {3pt} \\ {} & {} & \munderset {\textit {\_R} =\mathit {RootOf}\left (\mathit {a3} \,\textit {\_Z}^{3}+\mathit {a2} \,\textit {\_Z}^{2}+\mathit {a1} \textit {\_Z} +\mathit {a0} \right )}{\sum }\frac {\ln \left (y-\textit {\_R} \right )}{3 \textit {\_R}^{2} \mathit {a3} +2 \textit {\_R} \mathit {a2} +\mathit {a1}}=x +c_{1} \end {array} \]

Maple trace

`Methods for first order ODEs: 
--- Trying classification methods --- 
trying a quadrature 
trying 1st order linear 
trying Bernoulli 
trying separable 
<- separable successful`
 

Solution by Maple

Time used: 0.0 (sec). Leaf size: 30

dsolve(diff(y(x),x) = a0+a1*y(x)+a2*y(x)^2+a3*y(x)^3,y(x), singsol=all)
 

\[ x -\left (\int _{}^{y \left (x \right )}\frac {1}{\textit {\_a}^{3} \operatorname {a3} +\textit {\_a}^{2} \operatorname {a2} +\textit {\_a} \operatorname {a1} +\operatorname {a0}}d \textit {\_a} \right )+c_{1} = 0 \]

Solution by Mathematica

Time used: 0.112 (sec). Leaf size: 54

DSolve[y'[x]==a0+a1 y[x]+a2 y[x]^2+ a3 y[x]^3,y[x],x,IncludeSingularSolutions -> True]
 

\[ \text {Solve}\left [\text {RootSum}\left [\text {$\#$1}^3 \text {a3}+\text {$\#$1}^2 \text {a2}+\text {$\#$1} \text {a1}+\text {a0}\&,\frac {\log (y(x)-\text {$\#$1})}{3 \text {$\#$1}^2 \text {a3}+2 \text {$\#$1} \text {a2}+\text {a1}}\&\right ]=x+c_1,y(x)\right ] \]