Boundary Value Problems/Ordinary-Differential-Equations

In mathematics, an ordinary differential equation (or ODE) is a relation that contains functions of only one independent variable, and one or more of its derivatives with respect to that variable.

A simple example is Newton's second law of motion, which leads to the differential equation

${\displaystyle m{\frac {d^{2}x(t)}{dt^{2}}}=F(x(t)),\,}$

for the motion of a particle of constant mass m. In general, the force F depends upon the position of the particle x(t) at time t, and thus the unknown function x(t) appears on both sides of the differential equation, as is indicated in the notation F(x(t)).

Ordinary differential equations are distinguished from partial differential equations, which involve partial derivatives of several variables.

Ordinary differential equations arise in many different contexts including geometry, mechanics, astronomy and population modelling. Many famous mathematicians have studied differential equations and contributed to the field, including Isaac Newton|Newton, Leibniz, the Bernoulli family, Riccati, Alexis Claude Clairaut, d'Alembert and Euler.

Much study has been devoted to the solution of ordinary differential equations. In the case where the equation is linear, it can be solved by analytical methods. Unfortunately, most of the interesting differential equations are non-linear and, with a few exceptions, cannot be solved exactly. Approximate solutions are arrived at using computer approximations (see numerical ordinary differential equations).

The trajectory of a projectile launched from a cannon follows a curve determined by an ordinary differential equation that is derived from Newton's second law.

Ordinary differential equation

Let y be an unknown function

${\displaystyle y:\mathbb {R} \to \mathbb {R} }$

in x with ${\displaystyle y^{(n)}}$ the nth derivative of y, then an equation of the form

${\displaystyle F(x,y,y',\ \dots ,\ y^{(n-1)})=y^{(n)}}$

is called an ordinary differential equation (ODE) of order n; for vector valued functions,

${\displaystyle y:\mathbb {R} \to \mathbb {R} ^{m}}$,

it is called a system of ordinary differential equations of dimension m.

When a differential equation of order n has the form

${\displaystyle F\left(x,y,y',y'',\ \dots ,\ y^{(n)}\right)=0}$

it is called an implicit differential equation whereas the form

${\displaystyle F\left(x,y,y',y'',\ \dots ,\ y^{(n-1)}\right)=y^{(n)}}$

is called an explicit differential equation.

A differential equation not depending on x is called autonomous.

A differential equation is said to be linear if F can be written as a linear combination of the derivatives of y

${\displaystyle y^{(n)}=\sum _{i=0}^{n-1}a_{i}(x)y^{(i)}+r(x)}$

with ai(x) and r(x) continuous functions in x. The function r(x) is called the source term; if r(x)=0 then the linear differential equation is called homogeneous, otherwise it is called non-homogeneous or inhomogeneous.

Solutions

Given a differential equation

${\displaystyle F(x,y,y',\dots ,y^{(n)})=0}$

a function u: IRR, is called the solution or integral curve for F, if u is n-times differentiable on I, and

${\displaystyle F(x,u,u',\ \dots ,\ u^{(n)})=0\quad x\in I.}$

Given two solutions u: JRR and |v: IRR, u is called an extension of v ifIJ and

${\displaystyle u(x)=v(x)\quad x\in I.\,}$

A solution which has no extension is called a global solution.

A general solution of an n-th order equation is a solution containing n arbitrary variables, corresponding to n constant of integration|constants of integration. A particular solution is derived from the general solution by setting the constants to particular values, often chosen to fulfill set 'initial value problem|initial conditions or boundary conditions. A singular solution is a solution that can't be derived from the general solution.

Reduction to a first order system

Any differential equation of order n can be written as a system of n first-order differential equations. Given an explicit ordinary differential equation of order n and dimension 1,

${\displaystyle F\left(x,y,y',y'',\ \dots ,\ y^{(n-1)}\right)=y^{(n)}}$

we define a new family of unknown functions

${\displaystyle y_{n}:=y^{(n-1)}.\!}$

We can then rewrite the original differential equation as a system of differential equations with order 1 and dimension n.

${\displaystyle y_{1}^{'}=y_{2}}$
${\displaystyle y_{2}^{'}=y_{3}}$
${\displaystyle \vdots }$
${\displaystyle y_{n-1}^{'}=y_{n}}$
${\displaystyle y_{n}^{'}=F(y_{n},\dots ,y_{1},x).}$

which can be written concisely in vector notation as

${\displaystyle \mathbf {y} ^{'}=\mathbf {F} (\mathbf {y} ,x)}$

with

${\displaystyle \mathbf {y} :=(y,\ldots ,y_{n}).}$

Linear ordinary differential equations

A well understood particular class of differential equations is linear differential equations. We can always reduce an explicit linear differential equation of any order to a system of differential equation of order 1

${\displaystyle y_{i}'(x)=\sum _{j=1}^{n}a_{i,j}(x)y_{j}+b_{i}(x)\,\mathrm {,} \quad i=1,\ldots ,n}$

which we can write concisely using matrix and vector notation as

${\displaystyle \mathbf {y} ^{'}(x)=\mathbf {A} (x)\mathbf {y} (x)+\mathbf {b} (x)}$

with

${\displaystyle \mathbf {y} (x):=(y_{1}(x),\ldots ,y_{n}(x))}$
${\displaystyle \mathbf {b} (x):=(b_{1}(x),\ldots ,b_{n}(x))}$
${\displaystyle \mathbf {A} (x):=(a_{i,j}(x))\,\mathrm {,} \quad i,j=1,\ldots ,n.}$

Homogeneous equations

The set of solutions for a system of homogeneous linear differential equations of order 1 and dimension n

${\displaystyle \mathbf {y} ^{'}(x)=\mathbf {A} (x)\mathbf {y} (x)}$

forms an n-dimensional vector space. Given a basis for this vector space ${\displaystyle \mathbf {z} _{1}(x),\ldots ,\mathbf {z} _{n}(x)}$, which is called a fundamental system, every solution ${\displaystyle \mathbf {s} (x)}$ can be written as

${\displaystyle \mathbf {s} (x)=\sum _{i=1}^{n}c_{i}\mathbf {z} _{i}(x).}$

The n × n matrix

${\displaystyle \mathbf {Z} (x):=(\mathbf {z} _{1}(x),\ldots ,\mathbf {z} _{n}(x))}$

is called fundamental matrix. In general there is no method to explicitly construct a fundamental system, but if one solution is known d'Alembert reduction can be used to reduce the dimension of the differential equation by one.

Nonhomogeneous equations

The set of solutions for a system of inhomogeneous linear differential equations of order 1 and dimension n

${\displaystyle \mathbf {y} ^{'}(x)=\mathbf {A} (x)\mathbf {y} (x)+\mathbf {b} (x)}$

can be constructed by finding the fundamental system ${\displaystyle \mathbf {z} _{1}(x),\ldots ,\mathbf {z} _{n}(x)}$ to the corresponding homogeneous equation and one particular solution ${\displaystyle \mathbf {p} (x)}$ to the inhomogeneous equation. Every solution ${\displaystyle \mathbf {s} (x)}$ to nonhomogeneous equation can then be written as

${\displaystyle \mathbf {s} (x)=\sum _{i=1}^{n}c_{i}\mathbf {z} _{i}(x)+\mathbf {p} (x).}$

A particular solution to the nonhomogeneous equation can be found by the method of undetermined coefficients or the method of variation of parameters.

Fundamental systems for homogeneous equations with constant coefficients

If a system of homogeneous linear differential equations has constant coefficients

${\displaystyle \mathbf {y} ^{'}(x)=\mathbf {A} \mathbf {y} (x)}$

then we can explicitly construct a fundamental system. The fundamental system can be written as a matrix differential equation

${\displaystyle \mathbf {Y} ^{'}=\mathbf {A} \mathbf {Y} }$

with solution as a matrix exponential

${\displaystyle e^{x\mathbf {A} }}$

which is a fundamental matrix for the original differential equation. To explicitly calculate this expression we first transform A into Jordan normal form

${\displaystyle e^{x\mathbf {A} }=e^{x\mathbf {C} ^{-1}\mathbf {J} \mathbf {C} ^{1}}=\mathbf {C} ^{-1}e^{x\mathbf {J} }\mathbf {C} ^{1}}$

and then evaluate the Jordan blocks

${\displaystyle J_{i}={\begin{bmatrix}\lambda _{i}&1&\;&\;\\\;&\ddots &\ddots &\;\\\;&\;&\ddots &1\\\;&\;&\;&\lambda _{i}\end{bmatrix}}}$

of J separately as

${\displaystyle e^{x\mathbf {J_{i}} }=e^{\lambda _{i}x}{\begin{bmatrix}1&x&{\frac {x^{2}}{2}}&\dots &{\frac {x^{n-1}}{(n-1)!}}\\\;&\ddots &\ddots &\ddots &\vdots \\\;&\;&\ddots &\ddots &{\frac {x^{2}}{2}}\\\;&\;&\;&\ddots &x\\\;&\;&\;&\;&1\end{bmatrix}}.}$

Bibliography

• A. D. Polyanin and V. F. Zaitsev, Handbook of Exact Solutions for Ordinary Differential Equations (2nd edition), Chapman & Hall/CRC Press, Boca Raton, 2003. ISBN 1-58488-297-2
• A. D. Polyanin, V. F. Zaitsev, and A. Moussiaux, Handbook of First Order Partial Differential Equations, Taylor & Francis, London, 2002. ISBN 0-415-27267-X
• D. Zwillinger, Handbook of Differential Equations (3rd edition), Academic Press, Boston, 1997.
• Hartman, Philip, Ordinary Differential Equations, 2nd Ed., Society for Industrial & Applied Math, 2002. ISBN 0-89871-510-5.
• W. Johnson, A Treatise on Ordinary and Partial Differential Equations, John Wiley and Sons, 1913, in University of Michigan Historical Math Collection
• E.L. Ince, Ordinary Differential Equations, Dover Publications, 1958, ISBN 0486603490
• Witold Hurewicz Lectures on Ordinary Differential Equations, Dover Publications, ISBN 0-486-49510-8