Separation Of Variables Article Index for
Separation Of
Website Links For
Separation
 

Information About

Separation Of Variables





ORDINARY DIFFERENTIAL EQUATIONS (ODE)


Suppose a differential equation can be written in the form

: rac{d}{dx} f(x) = g(x)h(f(x)),\qquad\qquad (1)

which we can write more simply by letting y = f(x):

: rac{dy}{dx}=g(x)h(y)\qquad\qquad (1).

As long as ''h''(''y'') ≠ 0, we can rearrange terms to obtain:

:{dy \over h(y)} = {g(x)dx},

so that the two variables ''x'' and ''y'' have been separated.

Some who dislike the Leibniz Notation may prefer to write this as

: rac{1}{h(y)} rac{dy}{dx} = g(x),

but that fails to make it quite as obvious why this is called "separation of variables".

Integrating both sides of the equation with respect to x, we have

:\int rac{1}{h(y)} rac{dy}{dx} \, dx = \int g(x) \, dx + C, \qquad\qquad (2)

or equivalently,

:\int rac{1}{h(y)} \, dy = \int g(x) \, dx + C

because of the Substitution Rule For Integrals .

If one can evaluate the two integrals, one can find a solution to the differential equation. Observe that this process effectively allows us to treat the Derivative rac{dy}{dx} as a fraction which can be separated. This allows us to solve separable differential equations more conveniently, as demonstrated in the example below.

(Note that we do not need to use two Constants Of Integration , in equation (2) as in

:\int rac{1}{h(y)} \, dy + C_1 = \int g(x) \, dx + C_2,

because a single constant C = C_2 - C_1 is equivalent.)


Example


The ordinary differential equation

: rac{d}{dx}f(x)=f(x)(1-f(x))

may be written as

: rac{dy}{dx}=y(1-y).

If we let g(x) = 1 and h(y) = y(1-y), we can write the differential equation in the form of equation (1) above. Thus, the differential equation is separable.

As shown above, we can treat dy and dx as separate values, so that both sides of the equation may be multiplied by dx. Subsequently dividing both sides by y(1 - y), we have

: rac{dy}{y(1-y)}=dx.

At this point we have ''separated'' the variables ''x'' and ''y'' from each other, since ''x'' appears only on the right side of the equation and ''y'' only on the left.

Integrating both sides, we get

:\int rac{dy}{y(1-y)}=\int dx,

which, via Partial Fraction s, becomes

:\int rac{1}{y}+ rac{1}{1-y}\,dy=\int dx,

and then

:\ln y -\ln (1-y)=x+C

where ''C'' is the Constant Of Integration . A bit of Algebra gives a solution for ''y'':

:y= rac{1}{1+Be^{-x}}.

One may check our solution by taking the derivative with respect to x of the function we found, where ''B'' is an arbitrary constant. The result should be equal to our original problem.

Note that since we divided by (1 - y) and (y - 0) we must check to see whether the solutions y(x) = 0 and y(x) = 1 solve the s.


PARTIAL DIFFERENTIAL EQUATIONS


Given a Partial Differential Equation of a function

: F(x_1,x_2,\dots,x_n)

of ''n'' variables, it is sometimes useful to guess solution of the form

: F = F_1(x_1) \cdot F_2(x_2) \cdots F_n(x_n)

or

: F = f_1(x_1) + f_2(x_2) + \cdots + f_n(x_n)

which turns the partial differential equation (PDE) into a set of ODE s. Usually, each independent variable creates a ''separation constant'' that cannot be determined only from the equation itself.


Example (I)


Suppose ''F''(''x'', ''y'', ''z'') and the following PDE:

: rac{\partial F}{\partial x} + rac{\partial F}{\partial y} + rac{\partial F}{\partial z} = 0 \qquad\qquad (1)

We shall guess

: F(x,y,z) = X(x) + Y(y) + Z(z)\qquad\qquad (2)

thus making the equation (1) to

: rac{dX}{dx} + rac{dY}{dy} + rac{dZ}{dz} = 0

(since rac{\partial F}{\partial x} = rac{dX}{dx} ).

Now, since ''X'''(''x'') is dependent only at ''x'' and ''Y'''(''y'') is dependent only at ''y'' (so on for ''Z'''(''z'')) and that the equation (1) is true for every ''x, y, z'' it is clear that each one of the term is constant. More precisely,

: rac{dX}{dx} = c_1 \quad rac{dY}{dy} = c_2 \quad rac{dZ}{dz} = c_3\qquad\qquad (3)

were the constants ''c''1, ''c''2, ''c''3 satisfy

: c_1 + c_2 + c_3 = 0\qquad\qquad (4)

Eq. (3) is actually a set of three ODE s. In this case they are trivial and can be solved by simple Integration , giving:

: F(x,y,z) = c_1 x + c_2 y + c_3 z + c_4\qquad\qquad (5)

where the integration constant ''c''4 is determined by initial conditions.


Example (II)

Consider the differential equation

:
abla^2 v + \lambda v = {\partial^2 v \over \partial x^2} + {\partial^2 v \over \partial y^2} + \lambda v = 0.

First we seek solutions of the form

: v = X(x)Y(y).\,

Most solutions are not of that form, but other solutions are sums of (generally infinitely many) solutions of that form.

Substituting,

: {\partial^2\over\partial x^2}(X(x)Y(y))+{\partial^2\over\partial y^2} {Link without Title} +\lambda X(x)Y(y)=

: = X''(x)Y(y)+X(x)Y''(y)+\lambda X(x)Y(y)= 0\,

Divide throughout by ''X''(''x'')

: = {X''(x)Y(y) \over X(x)}+{X(x)Y''(y)\over X(x)}+{\lambda X(x)Y(y)\over X(x)}

: ={X''(x)Y(y) \over X(x)}+Y''(y)+\lambda Y(y) = 0

and then by ''Y''(''y'')

: ={X''(x)\over X(x)}+{Y''(y)+\lambda Y(y)\over Y(y)} = 0

Now ''X''′′(''x'')/''X''(''x'') is a function of ''x'' only, as is (''Y''′′(''y'')+λ''Y''(''y''))/''Y''(''y''), so there are separation constants so

: {X''(x)\over X(x)} = k = {Y''(y)+\lambda Y(y)\over Y(y)}

which splits up into Ordinary Differential Equation s

: {X''(x)\over X(x)} = k

: X''(x) - k X(x)=0\,

and

: {Y''(y)+\lambda Y(y)\over Y(y)} =k

: Y''(y)+(\lambda-k) Y(y) =0\,

which we can solve accordingly. If the equation as posed originally was a Boundary Value Problem , one would use the given boundary values. See that article for an example which uses boundary values.


EXTERNAL LINKS



BIBLIOGRAPHY

  • A. D. Polyanin and V. F. Zaitsev, ''Handbook of Exact Solutions for Ordinary Differential Equations'', Chapman & Hall/CRC Press, Boca Raton, 2003 (2nd edition). ISBN 1-58488-297-2

  • A. D. Polyanin, ''Handbook of Linear Partial Differential Equations for Engineers and Scientists'', Chapman & Hall/CRC Press, Boca Raton, 2002. ISBN 1-58488-299-9