Bernoulli differential equation

Differential equations
Scope
Fields Natural sciences Engineering Astronomy Physics Chemistry Biology Geology Applied mathematics Continuum mechanics Chaos theory Dynamical systems Social sciences Economics Population dynamics List of named differential equations
Classification
Types Ordinary Partial Differential-algebraic Integro-differential Fractional Linear Non-linear By variable type Dependent and independent variables Autonomous Coupled / Decoupled Exact Homogeneous / Nonhomogeneous Features Order Operator Notation
Relation to processes Difference (discrete analogue) Stochastic Stochastic partial Delay
Solution
Existence and uniqueness Picard–Lindelöf theorem Peano existence theorem Carathéodory's existence theorem Cauchy–Kowalevski theorem
General topics Initial conditions Boundary values Dirichlet Neumann Robin Cauchy problem Wronskian Phase portrait Lyapunov / Asymptotic / Exponential stability Rate of convergence Series / Integral solutions Numerical integration Dirac delta function
Solution methods Inspection Method of characteristics Euler Exponential response formula Finite difference (Crank–Nicolson) Finite element Infinite element Finite volume Galerkin Petrov–Galerkin Green's function Integrating factor Integral transforms Perturbation theory Runge–Kutta Separation of variables Undetermined coefficients Variation of parameters
People
List Isaac Newton Gottfried Leibniz Jacob Bernoulli Leonhard Euler Józef Maria Hoene-Wroński Joseph Fourier Augustin-Louis Cauchy George Green Carl David Tolmé Runge Martin Kutta Rudolf Lipschitz Ernst Lindelöf Émile Picard Phyllis Nicolson John Crank
v t e

In mathematics, an ordinary differential equation is called a Bernoulli differential equation if it is of the form

${\displaystyle y'+P(x)y=Q(x)y^{n},}$

where ${\displaystyle n}$ is a real number. Some authors allow any real ${\displaystyle n}$,^[1][2] whereas others require that ${\displaystyle n}$ not be 0 or 1.^[3][4] The equation was first discussed in a work of 1695 by Jacob Bernoulli, after whom it is named. The earliest solution, however, was offered by Gottfried Leibniz, who published his result in the same year and whose method is the one still used today.^[5]

Bernoulli equations are special because they are nonlinear differential equations with known exact solutions. A notable special case of the Bernoulli equation is the logistic differential equation.

Transformation to a linear differential equation

When ${\displaystyle n=0}$, the differential equation is linear. When ${\displaystyle n=1}$, it is separable. In these cases, standard techniques for solving equations of those forms can be applied. For ${\displaystyle n\neq 0}$ and ${\displaystyle n\neq 1}$, the substitution ${\displaystyle u=y^{1-n}}$ reduces any Bernoulli equation to a linear differential equation

${\displaystyle {\frac {du}{dx}}-(n-1)P(x)u=-(n-1)Q(x).}$

For example, in the case ${\displaystyle n=2}$, making the substitution ${\displaystyle u=y^{-1}}$ in the differential equation ${\displaystyle {\frac {dy}{dx}}+{\frac {1}{x}}y=xy^{2}}$ produces the equation ${\displaystyle {\frac {du}{dx}}-{\frac {1}{x}}u=-x}$, which is a linear differential equation.

Solution

Let ${\displaystyle x_{0}\in (a,b)}$ and

${\displaystyle {\begin{cases}z:(a,b)\rightarrow (0,\infty ),&{\text{if }}\alpha \in \mathbb {R} \smallsetminus \{1,2\},\\[4pt]z:(a,b)\rightarrow \mathbb {R} \smallsetminus \{0\},&{\text{if }}\alpha =2,\end{cases}}}$

be a solution of the linear differential equation

${\displaystyle z'(x)=(1-\alpha )P(x)z(x)+(1-\alpha )Q(x).}$

Then we have that ${\displaystyle y(x):=[z(x)]^{1/(1-\alpha )}}$ is a solution of

${\displaystyle y'(x)=P(x)y(x)+Q(x)y^{\alpha }(x)\ ,\ y(x_{0})=y_{0}:=[z(x_{0})]^{1/(1-\alpha )}.}$

And for every such differential equation, for all ${\displaystyle \alpha >0}$ we have ${\displaystyle y\equiv 0}$ as solution for ${\displaystyle y_{0}=0}$.

Example

Consider the Bernoulli equation

${\displaystyle y'-{\frac {2y}{x}}=-x^{2}y^{2}}$

(in this case, more specifically a Riccati equation). The constant function ${\displaystyle y=0}$ is a solution. Division by ${\displaystyle y^{2}}$ yields

${\displaystyle y'y^{-2}-{\frac {2}{x}}y^{-1}=-x^{2}}$

Changing variables gives the equations

${\displaystyle {\begin{aligned}u={\frac {1}{y}}\;&,~u'={\frac {-y'}{y^{2}}}\\[5pt]-u'-{\frac {2}{x}}u&=-x^{2}\\[5pt]u'+{\frac {2}{x}}u&=x^{2}\end{aligned}}}$

which can be solved using the integrating factor

${\displaystyle M(x)=e^{2\int {\frac {1}{x}}\,dx}=e^{2\ln x}=x^{2}.}$

Multiplying by ${\displaystyle M(x)}$,

${\displaystyle u'x^{2}+2xu=x^{4}.}$

The left side can be represented as the derivative of ${\displaystyle ux^{2}}$ by reversing the product rule. Applying the chain rule and integrating both sides with respect to ${\displaystyle x}$ results in the equations

${\displaystyle {\begin{aligned}\int \left(ux^{2}\right)'dx&=\int x^{4}\,dx\\[5pt]ux^{2}&={\frac {1}{5}}x^{5}+C\\[5pt]{\frac {1}{y}}x^{2}&={\frac {1}{5}}x^{5}+C\end{aligned}}}$

The solution for ${\displaystyle y}$ is

${\displaystyle y={\frac {x^{2}}{{\frac {1}{5}}x^{5}+C}}.}$

Notes

References

• Bernoulli, Jacob (1695), "Explicationes, Annotationes & Additiones ad ea, quae in Actis sup. de Curva Elastica, Isochrona Paracentrica, & Velaria, hinc inde memorata, & paratim controversa legundur; ubi de Linea mediarum directionum, alliisque novis", Acta Eruditorum. Cited in Hairer, Nørsett & Wanner (1993).
• Hairer, Ernst; Nørsett, Syvert Paul; Wanner, Gerhard (1993), Solving ordinary differential equations I: Nonstiff problems, Berlin, New York: Springer-Verlag, ISBN 978-3-540-56670-0.

External links

• Index of differential equations