Maths encyclopedia and lessons  
Search

Mathematics Encyclopedia and Lessons

 
     
 

Lessons

Popular
Subjects

algebra
arithmetic
calculus
equations
geometry
differential equations
trigonometry
number theory
probability theory
more
 

References

applied mathematics
mathematical games
mathematicians
more
 
 

Separation of variables

In mathematics, separation of variables is any of several methods of solving ordinary and partial differential equations.

Contents

Ordinary differential equations (ODE)

Occasionally a differential equation allows a separation of variables, which we here exemplify rather than define. The differential equation

\frac{df(x)}{dx}=f(x)(1-f(x))

may be written as

\frac{dy}{dx}=y(1-y).

Pretend that dy and dx are numbers, so that both sides of the equation may be multiplied by dx. Also divide both sides by y(1 − y). We get

\frac{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\frac{dy}{y(1-y)}=\int dx,

which, via partial fractions, becomes

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

and then

logey - loge(1 - y) = x + C.

A bit of algebra gives a solution for y:

y=\frac{1}{1+Be^{-x}}.

One may check that if B is any positive constant, this function satisfies the differential equation.

This process also exemplifies the utility of the Leibniz notation, in which dy and dx are thought of as infinitely small increments of y and x respectively.

Note that once dividing in (1 − y) and (y − 0) one must check if the solutions y(x) = 0 and y(x) = 1 solve the differential equation. See also: singular solutions.

Partial differential equations

Given a partial differential equation of a function

F(x1,x2,...,xn)

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 ODEs. 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:

\frac{\partial F}{\partial x} + \frac{\partial F}{\partial y} + \frac{\partial F}{\partial z} = 0 (1)

We shall guess

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

thus making the equation (1) to

\frac{dX}{dx} + \frac{dY}{dy} + \frac{dZ}{dz} = 0

(since \frac{\partial F}{\partial x} = \frac{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,

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

were the constants c1, c2, c3 satisfy

c1 + c2 + c3 = 0 (4)

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

F(x,y,z) = c1x + c2y + c3z + c4 (5)

where the integration constant c4 is determined by initial conditions.

Example (II)

Consider the differential equation

\nabla^2 v + \lambda v = {\partial^2 v \over \partial x^2} + {\partial^2 v \over \partial y^2} + \lambda v = 0

We suppose the solution is of the form

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

substituting,

{\partial^2\over\partial x^2}(X(x)Y(y))+{\partial^2\over\partial y}(X(x)Y(y))+\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 equations

{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.

01-04-2007 01:18:14
The contents of this article are licensed from Wikipedia.org
under the GNU Free Documentation License. How to see transparent copy