4.1 problem 5

4.1.1 Existence and uniqueness analysis
4.1.2 Solving as quadrature ode

Internal problem ID [12941]
Internal file name [OUTPUT/11594_Tuesday_November_07_2023_11_27_52_PM_60285105/index.tex]

Book: DIFFERENTIAL EQUATIONS by Paul Blanchard, Robert L. Devaney, Glen R. Hall. 4th edition. Brooks/Cole. Boston, USA. 2012
Section: Chapter 1. First-Order Differential Equations. Exercises section 1.5 page 71
Problem number: 5.
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 }-y \left (y-1\right ) \left (y-3\right )=0} \] With initial conditions \begin {align*} [y \left (0\right ) = 4] \end {align*}

4.1.1 Existence and uniqueness analysis

This is non linear first order ODE. In canonical form it is written as \begin {align*} y^{\prime } &= f(t,y)\\ &= y \left (y -1\right ) \left (y -3\right ) \end {align*}

The \(y\) domain of \(f(t,y)\) when \(t=0\) is \[ \{-\infty

The \(y\) domain of \(\frac {\partial f}{\partial y}\) when \(t=0\) is \[ \{-\infty

4.1.2 Solving as quadrature ode

Integrating both sides gives \begin {align*} \int \frac {1}{y \left (y -1\right ) \left (y -3\right )}d y &= \int {dt}\\ \frac {\ln \left (y -3\right )}{6}+\frac {\ln \left (y \right )}{3}-\frac {\ln \left (y -1\right )}{2}&= t +c_{1} \end {align*}

Initial conditions are used to solve for \(c_{1}\). Substituting \(t=0\) and \(y=4\) in the above solution gives an equation to solve for the constant of integration. \begin {align*} \frac {2 \ln \left (2\right )}{3}-\frac {\ln \left (3\right )}{2} = c_{1} \end {align*}

The solutions are \begin {align*} c_{1} = \frac {2 \ln \left (2\right )}{3}-\frac {\ln \left (3\right )}{2} \end {align*}

Trying the constant \begin {align*} c_{1} = \frac {2 \ln \left (2\right )}{3}-\frac {\ln \left (3\right )}{2} \end {align*}

Substituting \(c_{1}\) found above in the general solution gives \begin {align*} \frac {\ln \left (y -3\right )}{6}+\frac {\ln \left (y \right )}{3}-\frac {\ln \left (y -1\right )}{2} = t +\frac {2 \ln \left (2\right )}{3}-\frac {\ln \left (3\right )}{2} \end {align*}

The constant \(c_{1} = \frac {2 \ln \left (2\right )}{3}-\frac {\ln \left (3\right )}{2}\) gives valid solution.

Summary

The solution(s) found are the following \begin{align*} \tag{1} \frac {\ln \left (y-3\right )}{6}+\frac {\ln \left (y\right )}{3}-\frac {\ln \left (y-1\right )}{2} &= t +\frac {2 \ln \left (2\right )}{3}-\frac {\ln \left (3\right )}{2} \\ \end{align*}

Verification of solutions

\[ \frac {\ln \left (y-3\right )}{6}+\frac {\ln \left (y\right )}{3}-\frac {\ln \left (y-1\right )}{2} = t +\frac {2 \ln \left (2\right )}{3}-\frac {\ln \left (3\right )}{2} \] Verified OK.

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: 1.828 (sec). Leaf size: 133

dsolve([diff(y(t),t)=y(t)*(y(t)-1)*(y(t)-3),y(0) = 4],y(t), singsol=all)
 

\[ y \left (t \right ) = \frac {48 \left (\frac {{\mathrm e}^{6 t}}{3}-\frac {9}{16}\right ) \left (27-32 \,{\mathrm e}^{6 t}+8 \sqrt {16 \,{\mathrm e}^{12 t}-27 \,{\mathrm e}^{6 t}}\right )^{\frac {2}{3}}+48 \left (\left (27-32 \,{\mathrm e}^{6 t}+8 \sqrt {16 \,{\mathrm e}^{12 t}-27 \,{\mathrm e}^{6 t}}\right )^{\frac {1}{3}}+3\right ) \left ({\mathrm e}^{6 t}-\frac {\sqrt {16 \,{\mathrm e}^{12 t}-27 \,{\mathrm e}^{6 t}}}{4}-\frac {27}{16}\right )}{\left (27-32 \,{\mathrm e}^{6 t}+8 \sqrt {16 \,{\mathrm e}^{12 t}-27 \,{\mathrm e}^{6 t}}\right )^{\frac {2}{3}} \left (16 \,{\mathrm e}^{6 t}-27\right )} \]

Solution by Mathematica

Time used: 0.172 (sec). Leaf size: 132

DSolve[{y'[t]==y[t]*(y[t]-1)*(y[t]-3),{y[0]==4}},y[t],t,IncludeSingularSolutions -> True]
 

\[ y(t)\to \frac {3 i \left (\sqrt {3}+i\right ) \sqrt [3]{4 \sqrt {e^{6 t} \left (16 e^{6 t}-27\right )^3}+864 e^{6 t}-256 e^{12 t}-729}}{32 e^{6 t}-54}+\frac {9 \left (1+i \sqrt {3}\right )}{2 \sqrt [3]{4 \sqrt {e^{6 t} \left (16 e^{6 t}-27\right )^3}+864 e^{6 t}-256 e^{12 t}-729}}+1 \]