Artificial Intelligence 🤖
Differential equations
Ordinary differential equations
Nonlinear second order ODEs

Nonlinear second order ODEs

So far we have only looked at second order LODEs. The techniques we have used to solve the ODEs rely on the fact that the principle of superposition may be used to find a general solution. When we depart from linearity, finding general solutions to differential equations becomes a far more difficult task. In general, the second-order, nonlinear ODE (NLODE), may be expressed as,

y=f(x,y,y).y^{\prime \prime}=f\left(x, y, y^{\prime}\right) .

We can identify some obvious differences between a LODE and NLODE:

  1. A NLODE does not possess the property that a linear combination of two (or more, if higher order) linearly independent solutions is also a solution;
  2. It is much rarer that a solution to a NLODE may be expressed in exact form. It is therefore necessary to know how to use alternative methods such as:

(a) Qualitative methods (like the method in Section 12.6);

(b) Numerical methods.

  1. Generally, if we can solve a NLODE to obtain a solution, the resulting solution is not a general solution; in other words, NLODEs may possess singular solutions. The latter is a solution that cannot be described by the general integral obtained through solving the ODE.

Special cases

Very little is known about second order nonlinear ODEs but, in a few cases, a second order NLODE can be solved as a first order ODE through the use of a substitution. In particular, we look at two cases:

(i) y=f(x,y)y^{\prime \prime}=f\left(x, y^{\prime}\right), where yy is missing;

(ii) y=f(y,y)y^{\prime \prime}=f\left(y, y^{\prime}\right), where xx is missing.

Here, y=y(x)y=y(x) where xx is the independent variable and yy is the dependent variable. Case (ii) represents second order autonomous ODEs where ff does not depend explicitly on the independent variable.

Case 1: yy is missing

Such ODEs take the form,

y=f(x,y)y^{\prime \prime}=f\left(x, y^{\prime}\right)

where the substitution v=yv=y^{\prime} clearly transforms Eq. (13.149) to a first order ODE, namely,

v=f(x,v).v^{\prime}=f(x, v) .

Equation (13.150) can be solved to find a solution describing v(x)v(x). Then, y(x)y(x) is a solution to,

y=v(x)y^{\prime}=v(x)

Consider the following ODE,

y=2x(y)2.y^{\prime \prime}=2 x\left(y^{\prime}\right)^{2} .

We define v=yv=y^{\prime} which means v=yv^{\prime}=y^{\prime \prime}. In (13.152),

v=2xv2.v^{\prime}=2 x v^{2} .

Equation (13.153) is separable and can be solved to give,

v(x)=1x2+c12.v(x)=-\frac{1}{x^{2}+c_{1}^{2}} .

Using a uu-substitution, x=c1tanux=c_{1} \tan u, we obtain the following integral describing the solution to the NLODE, y(x)y(x),

y(x)=1c12tan2u+1)sec2uc1sec2udu.y(x)=-\int \frac{1}{c_{1}^{2} \underbrace{\left.\tan ^{2} u+1\right)}_{\sec ^{2} u}} c_{1} \sec ^{2} u d u .

Simplifying (13.155) and integrating gives,

y(x)=1c1tan1(xc1)+c2.y(x)=-\frac{1}{c_{1}} \tan ^{-1}\left(\frac{x}{c_{1}}\right)+c_{2} .

Case 2: xx is missing

In the case where the independent variable is missing, i.e.

y=f(y,y)y^{\prime \prime}=f\left(y, y^{\prime}\right)

we can still solve the NLODE using the same substitution as before but the substitution is used in a different way. We let y=vy^{\prime}=v and differentiate,

y=dvdx=dvdydydxv=vdvdy,y^{\prime \prime}=\frac{d v}{d x}=\frac{d v}{d y} \underbrace{\frac{d y}{d x}}_{v}=v \frac{d v}{d y},

where in (13.158), we have made use of the chain rule. The NLODE (13.157) becomes,

vdvdy=f(y,v)v \frac{d v}{d y}=f(y, v)

Now, we can think of (13.159) as an ODE for v(y)v(y), i.e. yy is the independent variable and vv is the dependent variable. If we can solve for v(y)v(y) then y(x)y(x) is a solution to,

y=v(y).y^{\prime}=v(y) .

Consider the following NLODE,

yy=(y)2.y y^{\prime \prime}=\left(y^{\prime}\right)^{2} .

Let y=vy^{\prime}=v; from (13.158) we have y=vdvdyy^{\prime \prime}=v \frac{d v}{d y}. Then, Eq. (13.161) is transformed to,

y(vdvdy)=v2.y\left(v \frac{d v}{d y}\right)=v^{2} .

Separating variables and integrating gives us:

v=c1yv=c_{1} y

But, v=yv=y^{\prime} so upon integrating once more, the solution to the NLODE (13.161) is:

y(x)=c2ec1x.y(x)=c_{2} e^{c_{1} x} .

Table 13.1: 'Simple' forms of f(x)f(x) and trial solution for ypy_{p}