Question

Find the general solution of the given differential equation on$\left(0,\infty \right)$

$4x^2{y}^{\text{'}\text{'}}+\left(16x^2+1\right)y=0$

Short answer

The general solution of the given differential equation using the Forbinous Method is

$y=x^{\frac{1}{2}}{C}_1{J}_0\left(2x\right)+C_2{Y}_0\left(2x\right)$

Step-by-step solution

Definition of Forbinous method for differential equations.

The Forbinous method is a method for finding an infinite series solution to a second-order ordinary differential equation with and near the regular singular point.

We have the differential equation $4x^2{y}^{\text{'}\text{'}}+\left(16x^2+1\right)y=0$

Which has a regular singular point at and has the form

$y^{\text{'}\text{'}}+\frac{1-2a}{x}{y}^{\text{'}}+\left(b^2{c}^2{x}^{2c-2}+\frac{a^2-p^2{c}^2}{x^2}\right)y=0\to \left(1\right)$

Which yields

$$\begin{gathered} \frac{4x}{4x}{y}^{\text{'}\text{'}}+\frac{16x^2+1}{4x}\times \frac{x}{x}y=0 \\ {y}^{\text{'}\text{'}}+\left(4+\frac{\frac{1}{4}}{x^2}\right)y=0 \\ {y}^{\text{'}\text{'}}+\frac{0}{x}{y}^{\text{'}}+\left(4x^0+\frac{\frac{1}{4}}{x^2}\right)y=0\to \left(2\right) \end{gathered}$$

Assuming the form of the general solution

Our aim is to find the general solution for the given differential equation. We can use Forbinous Method to solve this differential equation.

Assume the solution, according to Forbinous, to be in the form,

$y=\sum _{n=0}^{\infty }{c}_n{x}^{n+r}$

Using Bessel Function of First Kind and Second Kind

Find the first and second derivative, substitute it into the given differential equation and you would end up with two roots, which is r₁ = p and r₂= -p. Then you can find the two series solution for each case, based on the recurrence relation. Then you will identify the solution to have the form of Bessel Function of First Kind with a little modification and the series solutions are,

$$\begin{gathered} y_1=x^a{J}_p\left(b{x}^c\right) \\ {y}_2=x^a{J}_{-p}\left(b{x}^c\right) \end{gathered}$$

And you will find that the two solutions are linearly independent for integer. For , we can obtain another solution, which is linearly independent with , using Bessel Function of Second Kind, which is,

$$\begin{gathered} y_3=x^a{Y}_b\left(b{x}^c\right) \\ =x^a\frac{\cos p\pi J_p\left(b{x}^c\right)-J_{-p}\left(b{x}^c\right)}{\sin p\pi }\to \left(3\right) \end{gathered}$$

Generalized equation of the differential equation

Regarding the given differential equation, by comparing we have,

$$\begin{gathered} 1-2a=0\to a=\frac{1}{2} \\ 2c-2=0\to c=1 \\ {b}^2{c}^2=4\to b=2 \\ {a}^2-p^2{c}^2=\frac{1}{4}\to p=0 \end{gathered}$$

For p = 0, we have only a single series solution, which has the form

$y_1=x^{\frac{1}{2}}{J}_0\left(2x\right)$

We can obtain another solution which are linearly independent with using Bessel Function of Second Kind which is,

$y_2\left(x\right)=x^{\frac{1}{2}}{Y}_0\left(2x\right)$

By using the Superposition Principle, the general solution is

$$\begin{gathered} y=x^{\frac{1}{2}}{C}_1{J}_0\left(2x\right)+x^{\frac{1}{2}}{C}_2{Y}_0\left(2x\right) \\ y=x^{\frac{1}{2}}\left[C_1{J}_0\left(2x\right)+C_2{Y}_0\left(2x\right)\right] \end{gathered}$$

Where C₁C₂ are arbitrary constants.

Therefore the general solution of the given differential equation is

$y=x^{\frac{1}{2}}\left[C_1{J}_0\left(2x\right)+C_2{Y}_0\left(2x\right)\right]$