2.29 problem 29

Internal problem ID [8118]
Book : Second order enumerated odes
Section : section 2
Problem number : 29
Date solved : Monday, October 21, 2024 at 04:53:35 PM
CAS classification : [_Lienard]

Solve

\begin{align*} \cos \left (x \right )^{2} y^{\prime \prime }-2 \cos \left (x \right ) \sin \left (x \right ) y^{\prime }+y \cos \left (x \right )^{2}&=0 \end{align*}

2.29.1 Solved as second order ode using change of variable on y method 1

Time used: 0.437 (sec)

In normal form the given ode is written as

\begin{align*} y^{\prime \prime }+p \left (x \right ) y^{\prime }+q \left (x \right ) y&=0 \tag {2} \end{align*}

Where

\begin{align*} p \left (x \right )&=-\frac {\sin \left (2 x \right )}{\cos \left (x \right )^{2}}\\ q \left (x \right )&=1 \end{align*}

Calculating the Liouville ode invariant \(Q\) given by

\begin{align*} Q &= q - \frac {p'}{2}- \frac {p^2}{4} \\ &= 1 - \frac {\left (-\frac {\sin \left (2 x \right )}{\cos \left (x \right )^{2}}\right )'}{2}- \frac {\left (-\frac {\sin \left (2 x \right )}{\cos \left (x \right )^{2}}\right )^2}{4} \\ &= 1 - \frac {\left (-\frac {2 \cos \left (2 x \right )}{\cos \left (x \right )^{2}}-\frac {2 \sin \left (2 x \right ) \sin \left (x \right )}{\cos \left (x \right )^{3}}\right )}{2}- \frac {\left (\frac {\sin \left (2 x \right )^{2}}{\cos \left (x \right )^{4}}\right )}{4} \\ &= 1 - \left (-\frac {\cos \left (2 x \right )}{\cos \left (x \right )^{2}}-\frac {\sin \left (2 x \right ) \sin \left (x \right )}{\cos \left (x \right )^{3}}\right )-\frac {\sin \left (2 x \right )^{2}}{4 \cos \left (x \right )^{4}}\\ &= 2 \end{align*}

Since the Liouville ode invariant does not depend on the independent variable \(x\) then the transformation

\begin{align*} y = v \left (x \right ) z \left (x \right )\tag {3} \end{align*}

is used to change the original ode to a constant coefficients ode in \(v\). In (3) the term \(z \left (x \right )\) is given by

\begin{align*} z \left (x \right )&={\mathrm e}^{-\left (\int \frac {p \left (x \right )}{2}d x \right )}\\ &= e^{-\int \frac {-\frac {\sin \left (2 x \right )}{\cos \left (x \right )^{2}}}{2} }\\ &= \sec \left (x \right )\tag {5} \end{align*}

Hence (3) becomes

\begin{align*} y = v \left (x \right ) \sec \left (x \right )\tag {4} \end{align*}

Applying this change of variable to the original ode results in

\begin{align*} \cos \left (x \right ) \left (v^{\prime \prime }\left (x \right )+2 v \left (x \right )\right ) = 0 \end{align*}

Which is now solved for \(v \left (x \right )\).

The above ode can be simplified to

\begin{align*} v^{\prime \prime }\left (x \right )+2 v \left (x \right ) = 0 \end{align*}

This is second order with constant coefficients homogeneous ODE. In standard form the ODE is

\[ A v''(x) + B v'(x) + C v(x) = 0 \]

Where in the above \(A=1, B=0, C=2\). Let the solution be \(v \left (x \right )=e^{\lambda x}\). Substituting this into the ODE gives

\[ \lambda ^{2} {\mathrm e}^{x \lambda }+2 \,{\mathrm e}^{x \lambda } = 0 \tag {1} \]

Since exponential function is never zero, then dividing Eq(2) throughout by \(e^{\lambda x}\) gives

\[ \lambda ^{2}+2 = 0 \tag {2} \]

Equation (2) is the characteristic equation of the ODE. Its roots determine the general solution form.Using the quadratic formula

\[ \lambda _{1,2} = \frac {-B}{2 A} \pm \frac {1}{2 A} \sqrt {B^2 - 4 A C} \]

Substituting \(A=1, B=0, C=2\) into the above gives

\begin{align*} \lambda _{1,2} &= \frac {0}{(2) \left (1\right )} \pm \frac {1}{(2) \left (1\right )} \sqrt {0^2 - (4) \left (1\right )\left (2\right )}\\ &= \pm i \sqrt {2} \end{align*}

Hence

\begin{align*} \lambda _1 &= + i \sqrt {2}\\ \lambda _2 &= - i \sqrt {2} \end{align*}

Which simplifies to

\begin{align*} \lambda _1 &= i \sqrt {2} \\ \lambda _2 &= -i \sqrt {2} \\ \end{align*}

Since roots are complex conjugate of each others, then let the roots be

\[ \lambda _{1,2} = \alpha \pm i \beta \]

Where \(\alpha =0\) and \(\beta =\sqrt {2}\). Therefore the final solution, when using Euler relation, can be written as

\[ v \left (x \right ) = e^{\alpha x} \left ( c_1 \cos (\beta x) + c_2 \sin (\beta x) \right ) \]

Which becomes

\[ v \left (x \right ) = e^{0}\left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \]

Or

\[ v \left (x \right ) = c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right ) \]

Will add steps showing solving for IC soon.

Now that \(v \left (x \right )\) is known, then

\begin{align*} y&= v \left (x \right ) z \left (x \right )\\ &= \left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \left (z \left (x \right )\right )\tag {7} \end{align*}

But from (5)

\begin{align*} z \left (x \right )&= \sec \left (x \right ) \end{align*}

Hence (7) becomes

\begin{align*} y = \left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \sec \left (x \right ) \end{align*}

Will add steps showing solving for IC soon.

2.29.2 Solved as second order ode adjoint method

Time used: 2.852 (sec)

In normal form the ode

\begin{align*} \cos \left (x \right )^{2} y^{\prime \prime }-2 \cos \left (x \right ) \sin \left (x \right ) y^{\prime }+y \cos \left (x \right )^{2} = 0 \tag {1} \end{align*}

Becomes

\begin{align*} y^{\prime \prime }+p \left (x \right ) y^{\prime }+q \left (x \right ) y&=r \left (x \right ) \tag {2} \end{align*}

Where

\begin{align*} p \left (x \right )&=-2 \tan \left (x \right )\\ q \left (x \right )&=1\\ r \left (x \right )&=0 \end{align*}

The Lagrange adjoint ode is given by

\begin{align*} \xi ^{''}-(\xi \, p)'+\xi q &= 0\\ \xi ^{''}-\left (-2 \tan \left (x \right ) \xi \left (x \right )\right )' + \left (\xi \left (x \right )\right ) &= 0\\ 2 \tan \left (x \right ) \xi ^{\prime }\left (x \right )+\xi \left (x \right )+\xi ^{\prime \prime }\left (x \right )+2 \sec \left (x \right )^{2} \xi \left (x \right )&= 0 \end{align*}

Which is solved for \(\xi (x)\). In normal form the given ode is written as

\begin{align*} \xi ^{\prime \prime }+p \left (x \right ) \xi ^{\prime }+q \left (x \right ) \xi &=0 \tag {2} \end{align*}

Where

\begin{align*} p \left (x \right )&=2 \tan \left (x \right )\\ q \left (x \right )&=1+2 \sec \left (x \right )^{2} \end{align*}

Calculating the Liouville ode invariant \(Q\) given by

\begin{align*} Q &= q - \frac {p'}{2}- \frac {p^2}{4} \\ &= 1+2 \sec \left (x \right )^{2} - \frac {\left (2 \tan \left (x \right )\right )'}{2}- \frac {\left (2 \tan \left (x \right )\right )^2}{4} \\ &= 1+2 \sec \left (x \right )^{2} - \frac {\left (2+2 \tan \left (x \right )^{2}\right )}{2}- \frac {\left (4 \tan \left (x \right )^{2}\right )}{4} \\ &= 1+2 \sec \left (x \right )^{2} - \left (1+\tan \left (x \right )^{2}\right )-\tan \left (x \right )^{2}\\ &= 2 \end{align*}

Since the Liouville ode invariant does not depend on the independent variable \(x\) then the transformation

\begin{align*} \xi = v \left (x \right ) z \left (x \right )\tag {3} \end{align*}

is used to change the original ode to a constant coefficients ode in \(v\). In (3) the term \(z \left (x \right )\) is given by

\begin{align*} z \left (x \right )&={\mathrm e}^{-\left (\int \frac {p \left (x \right )}{2}d x \right )}\\ &= e^{-\int \frac {2 \tan \left (x \right )}{2} }\\ &= \cos \left (x \right )\tag {5} \end{align*}

Hence (3) becomes

\begin{align*} \xi = v \left (x \right ) \cos \left (x \right )\tag {4} \end{align*}

Applying this change of variable to the original ode results in

\begin{align*} \left (2 v \left (x \right )+v^{\prime \prime }\left (x \right )\right ) \cos \left (x \right ) = 0 \end{align*}

Which is now solved for \(v \left (x \right )\).

The above ode can be simplified to

\begin{align*} 2 v \left (x \right )+v^{\prime \prime }\left (x \right ) = 0 \end{align*}

This is second order with constant coefficients homogeneous ODE. In standard form the ODE is

\[ A v''(x) + B v'(x) + C v(x) = 0 \]

Where in the above \(A=1, B=0, C=2\). Let the solution be \(v \left (x \right )=e^{\lambda x}\). Substituting this into the ODE gives

\[ \lambda ^{2} {\mathrm e}^{x \lambda }+2 \,{\mathrm e}^{x \lambda } = 0 \tag {1} \]

Since exponential function is never zero, then dividing Eq(2) throughout by \(e^{\lambda x}\) gives

\[ \lambda ^{2}+2 = 0 \tag {2} \]

Equation (2) is the characteristic equation of the ODE. Its roots determine the general solution form.Using the quadratic formula

\[ \lambda _{1,2} = \frac {-B}{2 A} \pm \frac {1}{2 A} \sqrt {B^2 - 4 A C} \]

Substituting \(A=1, B=0, C=2\) into the above gives

\begin{align*} \lambda _{1,2} &= \frac {0}{(2) \left (1\right )} \pm \frac {1}{(2) \left (1\right )} \sqrt {0^2 - (4) \left (1\right )\left (2\right )}\\ &= \pm i \sqrt {2} \end{align*}

Hence

\begin{align*} \lambda _1 &= + i \sqrt {2}\\ \lambda _2 &= - i \sqrt {2} \end{align*}

Which simplifies to

\begin{align*} \lambda _1 &= i \sqrt {2} \\ \lambda _2 &= -i \sqrt {2} \\ \end{align*}

Since roots are complex conjugate of each others, then let the roots be

\[ \lambda _{1,2} = \alpha \pm i \beta \]

Where \(\alpha =0\) and \(\beta =\sqrt {2}\). Therefore the final solution, when using Euler relation, can be written as

\[ v \left (x \right ) = e^{\alpha x} \left ( c_1 \cos (\beta x) + c_2 \sin (\beta x) \right ) \]

Which becomes

\[ v \left (x \right ) = e^{0}\left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \]

Or

\[ v \left (x \right ) = c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right ) \]

Will add steps showing solving for IC soon.

Now that \(v \left (x \right )\) is known, then

\begin{align*} \xi &= v \left (x \right ) z \left (x \right )\\ &= \left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \left (z \left (x \right )\right )\tag {7} \end{align*}

But from (5)

\begin{align*} z \left (x \right )&= \cos \left (x \right ) \end{align*}

Hence (7) becomes

\begin{align*} \xi = \left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \cos \left (x \right ) \end{align*}

Will add steps showing solving for IC soon.

The original ode (2) now reduces to first order ode

\begin{align*} \xi \left (x \right ) y^{\prime }-y \xi ^{\prime }\left (x \right )+\xi \left (x \right ) p \left (x \right ) y&=\int \xi \left (x \right ) r \left (x \right )d x\\ y^{\prime }+y \left (p \left (x \right )-\frac {\xi ^{\prime }\left (x \right )}{\xi \left (x \right )}\right )&=\frac {\int \xi \left (x \right ) r \left (x \right )d x}{\xi \left (x \right )}\\ y^{\prime }+y \left (-2 \tan \left (x \right )-\frac {\left (-c_1 \sqrt {2}\, \sin \left (\sqrt {2}\, x \right )+c_2 \sqrt {2}\, \cos \left (\sqrt {2}\, x \right )\right ) \cos \left (x \right )-\left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \sin \left (x \right )}{\left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \cos \left (x \right )}\right )&=0 \end{align*}

Which is now a first order ode. This is now solved for \(y\). In canonical form a linear first order is

\begin{align*} y^{\prime } + q(x)y &= p(x) \end{align*}

Comparing the above to the given ode shows that

\begin{align*} q(x) &=-\frac {2 \tan \left (x \right ) \cos \left (x \right ) \sin \left (\sqrt {2}\, x \right ) c_2 +2 \tan \left (x \right ) \cos \left (x \right ) \cos \left (\sqrt {2}\, x \right ) c_1 -\cos \left (x \right ) \sqrt {2}\, \sin \left (\sqrt {2}\, x \right ) c_1 +\cos \left (x \right ) \sqrt {2}\, \cos \left (\sqrt {2}\, x \right ) c_2 -\sin \left (x \right ) \sin \left (\sqrt {2}\, x \right ) c_2 -\sin \left (x \right ) \cos \left (\sqrt {2}\, x \right ) c_1}{\left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \cos \left (x \right )}\\ p(x) &=0 \end{align*}

The integrating factor \(\mu \) is

\begin{align*} \mu &= e^{\int {q\,dx}}\\ &= {\mathrm e}^{\int -\frac {2 \tan \left (x \right ) \cos \left (x \right ) \sin \left (\sqrt {2}\, x \right ) c_2 +2 \tan \left (x \right ) \cos \left (x \right ) \cos \left (\sqrt {2}\, x \right ) c_1 -\cos \left (x \right ) \sqrt {2}\, \sin \left (\sqrt {2}\, x \right ) c_1 +\cos \left (x \right ) \sqrt {2}\, \cos \left (\sqrt {2}\, x \right ) c_2 -\sin \left (x \right ) \sin \left (\sqrt {2}\, x \right ) c_2 -\sin \left (x \right ) \cos \left (\sqrt {2}\, x \right ) c_1}{\left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \cos \left (x \right )}d x}\\ &= \frac {\cos \left (x \right )}{c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )} \end{align*}

The ode becomes

\begin{align*} \frac {\mathop {\mathrm {d}}}{ \mathop {\mathrm {d}x}} \mu y &= 0 \\ \frac {\mathop {\mathrm {d}}}{ \mathop {\mathrm {d}x}} \left (\frac {y \cos \left (x \right )}{c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )}\right ) &= 0 \end{align*}

Integrating gives

\begin{align*} \frac {y \cos \left (x \right )}{c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )}&= \int {0 \,dx} + c_3 \\ &=c_3 \end{align*}

Dividing throughout by the integrating factor \(\frac {\cos \left (x \right )}{c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )}\) gives the final solution

\[ y = \left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \sec \left (x \right ) c_3 \]

Hence, the solution found using Lagrange adjoint equation method is

\begin{align*} y &= \left (c_1 \cos \left (\sqrt {2}\, x \right )+c_2 \sin \left (\sqrt {2}\, x \right )\right ) \sec \left (x \right ) c_3 \\ \end{align*}

Will add steps showing solving for IC soon.

2.29.3 Maple step by step solution

2.29.4 Maple trace

Methods for second order ODEs:
 

2.29.5 Maple dsolve solution

Solving time : 0.004 (sec)
Leaf size : 24

dsolve(cos(x)^2*diff(diff(y(x),x),x)-2*cos(x)*sin(x)*diff(y(x),x)+y(x)*cos(x)^2 = 0, 
       y(x),singsol=all)
 

\[ y = \sec \left (x \right ) \left (c_1 \sin \left (\sqrt {2}\, x \right )+c_2 \cos \left (\sqrt {2}\, x \right )\right ) \]

2.29.6 Mathematica DSolve solution

Solving time : 0.096 (sec)
Leaf size : 51

DSolve[{Cos[x]^2*D[y[x],{x,2}]-2*Cos[x]*Sin[x]*D[y[x],x]+y[x]*Cos[x]^2==0,{}}, 
       y[x],x,IncludeSingularSolutions->True]
 

\[ y(x)\to \frac {1}{4} e^{-i \sqrt {2} x} \left (4 c_1-i \sqrt {2} c_2 e^{2 i \sqrt {2} x}\right ) \sec (x) \]