تفاضليه عاديه
Ordinary differential equation

It has been suggested that Strang splitting be merged into this article. (Discuss) Proposed since December 2014. 
Differential equations 

Navier–Stokes differential equations used to simulate airflow around an obstruction.

Classification 
Solution 
In mathematics, an ordinary differential equation or ODE is an equation containing a function of one independent variable and its derivatives. The term "ordinary" is used in contrast with the term partial differential equation which may be with respect to more thanone independent variable.
Linear differential equations, which have solutions that can be added and multiplied by coefficients, are welldefined and understood, and exact closedform solutions are obtained. By contrast, ODEs that lack additive solutions are nonlinear, and solving them is far more intricate, as one can rarely represent them by elementary functions in closed form: Instead, exact and analytic solutions of ODEs are in series or integral form. Graphical and numerical methods, applied by hand or by computer, may approximate solutions of ODEs and perhaps yield useful information, often sufficing in the absence of exact, analytic solutions.
Contents
[hide]
Background[edit]
Ordinary differential equations (ODEs) arise in many different contexts throughout mathematics and science (social and natural) one way or another, because when describing changes mathematically, the most accurate way uses differentials and derivatives (related, though not quite the same). Since various differentials, derivatives, and functions become inevitably related to each other via equations, a differential equation is the result, describing dynamical phenomena, evolution, and variation. Often, quantities are defined as the rate of change of other quantities (time derivatives), or gradients of quantities, which is how they enter differential equations.
Specific mathematical fields include geometry and analytical mechanics. Scientific fields include much of physics and astronomy(celestial mechanics), geology (weather modelling), chemistry (reaction rates),^{[1]} biology (infectious diseases, genetic variation),ecology and population modelling (population competition), economics (stock trends, interest rates and the market equilibrium price changes).
Many mathematicians have studied differential equations and contributed to the field, including Newton, Leibniz, the Bernoulli family, Riccati, Clairaut, d'Alembert, and Euler.
A simple example is Newton's second law of motion — the relationship between the displacement x and the time t of the object under the force F, which leads to the differential equation
for the motion of a particle of constant mass m. In general, F depends on the position x(t) of the particle at time t, and so the unknown function x(t) appears on both sides of the differential equation, as is indicated in the notation F(x(t)).^{[2]}^{[3]}^{[4]}^{[5]}
Definitions[edit]
In what follows, let y be a dependent variable and x an independent variable, so that y = y(x) is an unknown function in x. The notation for differentiation varies depending upon the author and upon which notation is most useful for the task at hand. In this context, the Leibniz's notation (dy/dx,d^{2}y/dx^{2},...d^{n}y/dx^{n}) is useful for differentials and whenintegration is to be done, whereas Newton's and Lagrange's notation (y′,y′′, ... y^{(n)}) is useful for representing derivatives of any order compactly.
General definition of an ODE[edit]
Let F be a given function of x, y, and derivatives of y. Then an equation of the form
is called an explicit ordinary differential equation of order n.^{[6]}^{[7]}
More generally, an implicit ordinary differential equation of order n takes the form:^{[8]}
There are further classifications:
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:
where a_{i}(x) and r(x) are continuous functions in x.^{[9]}^{[10]}^{[11]} Nonlinear equations cannot be written in this form. The function r(x) is called the source term, leading to two further important classifications:^{[12]}^{[13]}
Homogeneous: If r(x) = 0, and consequently one "automatic" solution is the trivial solution, y = 0. The solution of a linear homogeneous equation is a complementary function, denoted here by y_{c}.
Nonhomogeneous (or inhomogeneous): If r(x) ≠ 0. The additional solution to the complementary function is the particular integral, denoted here by y_{p}.
The general solution to a linear equation can be written as y = y_{c} + y_{p}.
System of ODEs[edit]
A number of coupled differential equations form a system of equations. If y is a vector whose elements are functions; y(x) = [y_{1}(x), y_{2}(x),..., y_{m}(x)], and F is a vector valued function of y and its derivatives, then
is an explicit system of ordinary differential equations of order or dimension m. In column vector form:
These are not necessarily linear. The implicit analogue is:
where 0 = (0, 0,... 0) is the zero vector. In matrix form
For a system of the form , some sources also require that the Jacobian matrix be nonsingular in order to call this an implicit ODE [system]; an implicit ODE system satisfying this Jacobian nonsingularity condition can be transformed into an explicit ODE system. In the same sources, implicit ODE systems with a singular Jacobian are termed differential algebraic equations (DAEs). This distinction is not merely one of terminology; DAEs have fundamentally different characteristics and are generally more involved to solve that (nonsigular) ODE systems.^{[14]}^{[15]} Presumably for additional derivatives, the Hessian matrix and so forth are also assumed nonsingular according to this scheme,^{[citation needed]} although note that any ODE of order greater than one can [and usually is] rewritten as system of ODEs of first order,^{[16]} which makes the Jacobian singularity criterion sufficient for this taxonomy to be comprehensive at all orders.
Solutions[edit]
Given a differential equation
a function u: I ⊂ R → R is called the solution or integral curve for F, if u is ntimes differentiable on I, and
Given two solutions u: J ⊂ R → R and v: I ⊂ R → R, u is called an extension of v if I ⊂ J and
A solution that has no extension is called a maximal solution. A solution defined on all of R is called a global solution.
A general solution of an nthorder equation is a solution containing n arbitrary independent 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 conditions or boundary conditions'.^{[17]} A singular solution is a solution that cannot be obtained by assigning definite values to the arbitrary constants in the general solution.^{[18]}
Theories of ODEs[edit]
Singular solutions[edit]
The theory of singular solutions of ordinary and partial differential equations was a subject of research from the time of Leibniz, but only since the middle of the nineteenth century did it receive special attention. A valuable but littleknown work on the subject is that of Houtain (1854). Darboux (starting in 1873) was a leader in the theory, and in the geometric interpretation of these solutions he opened a field worked by various writers, notable ones being Casorati and Cayley. To the latter is due (1872) the theory of singular solutions of differential equations of the first order as accepted circa 1900.
Reduction to quadratures[edit]
The primitive attempt in dealing with differential equations had in view a reduction to quadratures. As it had been the hope of eighteenthcentury algebraists to find a method for solving the general equation of the nth degree, so it was the hope of analysts to find a general method for integrating any differential equation. Gauss (1799) showed, however, that the differential equation meets its limitations very soon unless complex numbers are introduced. Hence, analysts began to substitute the study of functions, thus opening a new and fertile field. Cauchy was the first to appreciate the importance of this view. Thereafter, the real question was to be not whether a solution is possible by means of known functions or their integrals but whether a given differential equation suffices for the definition of a function of the independent variable or variables, and, if so, what are the characteristic properties of this function.
Fuchsian theory[edit]
Two memoirs by Fuchs (Crelle, 1866, 1868), inspired a novel approach, subsequently elaborated by Thomé and Frobenius. Collet was a prominent contributor beginning in 1869, although his method for integrating a nonlinear system was communicated to Bertrand in 1868. Clebsch (1873) attacked the theory along lines parallel to those followed in his theory of Abelian integrals. As the latter can be classified according to the properties of the fundamental curve that remains unchanged under a rational transformation, so Clebsch proposed to classify the transcendent functions defined by the differential equations according to the invariant properties of the corresponding surfaces f = 0 under rational onetoone transformations.
Lie's theory[edit]
From 1870, Sophus Lie's work put the theory of differential equations on a more satisfactory foundation. He showed that the integration theories of the older mathematicians can, by the introduction of what are now called Lie groups, be referred to a common source, and that ordinary differential equations that admit the same infinitesimal transformationspresent comparable difficulties of integration. He also emphasized the subject of transformations of contact.
Lie's group theory of differential equations has been certified, namely: (1) that it unifies the many ad hoc methods known for solving differential equations, and (2) that it provides powerful new ways to find solutions. The theory has applications to both ordinary and partial differential equations.^{[19]}
A general approach to solve DEs uses the symmetry property of differential equations, the continuous infinitesimal transformations of solutions to solutions (Lie theory). Continuous group theory, Lie algebras, and differential geometry are used to understand the structure of linear and nonlinear (partial) differential equations for generating integrable equations, to find its Lax pairs, recursion operators, Bäcklund transform, and finally finding exact analytic solutions to the DE.
Symmetry methods have been recognized to study differential equations, arising in mathematics, physics, engineering, and many other disciplines.
Sturm–Liouville theory[edit]
Sturm–Liouville theory is a theory of eigenvalues and eigenfunctions of linear operators defined in terms of secondorder homogeneous linear equations, and is useful in the analysis of certain partial differential equations.
Existence and uniqueness of solutions[edit]
There are several theorems that establish existence and uniqueness of solutions to initial value problems involving ODEs both locally and globally. The two main theorems are

Theorem Assumption Conclusion Peano existence theorem F continuous local existence only Picard–Lindelöf theorem F Lipschitz continuous local existence and uniqueness
which are both local results.
Note that uniqueness theorems like the Lipschitz one above do not apply to DAE systems, which may have multiple solutions stemming from their (nonlinear) algebraic part alone.^{[20]}
Local existence and uniqueness theorem simplified[edit]
The theorem can be stated simply as follows.^{[21]} For the equation and initial value problem:
if F and ∂F/∂y are continuous in a closed rectangle
in the xy plane, where a and b are real (symbolically: a, b ∈ ℝ) and × denotes the cartesian product, square brackets denote closed intervals, then there is an interval
for some h ∈ ℝ where the solution to the above equation and initial value problem can be found. That is, there is a solution and it is unique. Since there is no restriction on F to be linear, this applies to nonlinear equations that take the form F(x, y), and it can also be applied to systems of equations.
Global uniqueness and maximum domain of solution[edit]
When the hypotheses of the Picard–Lindelöf theorem are satisfied, then local existence and uniqueness can be extended to a global result. More precisely:^{[22]}
For each initial condition (x_{0}, y_{0}) there exists a unique maximum (possibly infinite) open interval
such that any solution that satisfies this initial condition is a restriction of the solution that satisfies this initial condition with domain I_{max}.
In the case that , there are exactly two possibilities
 explosion in finite time:
 leaves domain of definition:
where Ω is the open set in which F is defined, and is its boundary.
Note that the maximum domain of the solution
 is always an interval (to have uniqueness)
 may be smaller than ℝ
 may depend on the specific choice of (x_{0}, y_{0}).
 Example
This means that F(x, y) = y^{2}, which is C^{1} and therefore locally Lipschitz continuous, satisfying the Picard–Lindelöf theorem.
Even in such a simple setting, the maximum domain of solution cannot be all ℝ, since the solution is
which has maximum domain:
This shows clearly that the maximum interval may depend on the initial conditions. The domain of y could be taken as being , but this would lead to a domain that is not an interval, so that the side opposite to the initial condition would be disconnected from the initial condition, and therefore not uniquely determined by it.
The maximum domain is not ℝ because
which is one of the two possible cases according to the above theorem.
Reduction of order[edit]
Differential equations can usually be solved more easily if the order of the equation can be reduced.
Reduction to a firstorder system[edit]
Any differential equation of order n,
can be written as a system of n firstorder differential equations by defining a new family of unknown functions
for i = 1, 2,... n. The ndimensional system of firstorder coupled differential equations is then
more compactly in vector notation:
where
Summary of exact solutions[edit]
Some differential equations have solutions that can be written in an exact and closed form. Several important classes are given here.
In the table below, P(x), Q(x), P(y), Q(y), and M(x,y), N(x,y) are any integrable functions of x, y, and b and c are real given constants, and C_{1}, C_{2},... are arbitrary constants (complex in general). The differential equations are in their equivalent and alternative forms that lead to the solution through integration.
In the integral solutions, λ and ε are dummy variables of integration (the continuum analogues of indices in summation), and the notation ∫^{x}F(λ)dλ just means to integrate F(λ) with respect to λ, then after the integration substitute λ = x, without adding constants (explicitly stated).

Differential equation Solution method General solution Separable equations Firstorder, separable in x and y (general case, see below for special cases)^{[23]} Separation of variables (divide byP_{2}Q_{1}). Firstorder, separable in x^{[21]} Direct integration. Firstorder, autonomous, separable in y^{[21]} Separation of variables (divide byF). Firstorder, separable in x and y^{[21]} Integrate throughout. General firstorder equations Firstorder, homogeneous^{[21]} Set y = ux, then solve by separation of variables in u and x. Firstorder, separable^{[23]} Separation of variables (divide byxy). If N = M, the solution is xy = C.
Exact differential, firstorder^{[21]} where
Integrate throughout. where Y(y) and X(x) are functions from the integrals rather than constant values, which are set to make the final function F(x, y) satisfy the initial equation.
Inexact differential, firstorder^{[21]} where
Integration factor μ(x, y) satisfying If μ(x, y) can be found: General secondorder equations Secondorder, autonomous^{[24]} Multiply equation by 2dy/dx, substitute , then integrate twice. Linear equations (up to nth order) Firstorder, linear, inhomogeneous, function coefficients^{[21]} Integrating factor: . Secondorder, linear, inhomogeneous, constant coefficients^{[25]} Complementary function y_{c}: assume y_{c} = e^{αx}, substitute and solve polynomial in α, to find thelinearly independent functions . Particular integral y_{p}: in general themethod of variation of parameters, though for very simple r(x) inspection may work.^{[21]}
If b^{2} > 4c, then:
If b^{2} = 4c, then:
If b^{2} < 4c, then:
nthorder, linear, inhomogeneous, constant coefficients^{[25]} Complementary function y_{c}: assume y_{c} = e^{αx}, substitute and solve polynomial in α, to find thelinearly independent functions . Particular integral y_{p}: in general themethod of variation of parameters, though for very simple r(x) inspection may work.^{[21]}
Since α_{j} are the solutions of the polynomial of degree n: , then:
for α_{j} all different,
for each root α_{j} repeated k_{j} times,
for some α_{j} complex, then setting α = χ_{j} + iγ_{j}, and using Euler's formula, allows some terms in the previous results to be written in the form
where ϕ_{j} is an arbitrary constant (phase shift).
Software for ODE solving[edit]
 Maxima computer algebra system (GPL)
 COPASI a free (Artistic License 2.0) software package for the integration and analysis of ODEs.
 MATLAB a Technical Computing Software (MATrix LABoratory)
 GNU Octave a highlevel language, primarily intended for numerical computations.
 Scilab an open source software for numerical computation.
 Maple
 Mathematica
 Julia (programming language)
 SciPy a Python package that includes an ODE integration module.