INTEGRATION BY PARTS


In calculus, and more generally in mathematical analysis, 'integration by parts' is a rule that transforms the integral of products of functions into other, possibly simpler, integrals. The rule arises from the product rule of differentiation.

Contents
The rule
Examples
The ILATE rule
Recursive integration by parts
Tabular integration by parts
Higher dimensions
Cultural references
References
External links

The rule


Suppose ''f''(''x'') and ''g''(''x'') are two continuously differentiable functions. Then the integration by parts rule states that given an interval with endpoints ''a'', ''b'', one has
:int_a^b f(x) g'(x),dx = left[ f(x) g(x)
ight]_{a}^{b} - int_a^b f'(x) g(x),dx
where we use the common notation
:left[ f(x) g(x)
ight]_{a}^{b} = f(b) g(b) - f(a) g(a).
The rule is shown to be true by using the product rule for derivatives and the fundamental theorem of calculus. Thus
:{|
|-
| f(b)g(b) - f(a)g(a),
|= int_a^b rac{d}{dx} ( f(x) g(x) ) , dx
|-
|
|=int_a^b f'(x) g(x) , dx + int_a^b f(x) g'(x) , dx.
|}
In the traditional calculus curriculum, this rule is often stated using indefinite integrals in the form
:int f(x) g'(x),dx = f(x) g(x) - int f'(x) g(x),dx,
or in an even shorter form, if we let ''u'' = ''f''(''x''), ''v'' = ''g''(''x'') and the differentials ''du'' = ''f'' ′(''x'') ''dx'' and ''dv'' = ''g''′(''x'') ''dx'', then it is in the form in which it is most often seen:
:int u,dv = u v - int v,du.
Note that the original integral contains the derivative of ''g''; in order to be able to apply the rule, the antiderivative ''g'' must be found, and then the resulting integral ∫''g'' ''f''′ d''x'' must be evaluated.
One can also formulate a discrete analogue for sequences, called summation by parts.
An alternative notation has the advantage that the factors of the original expression are identified as ''f'' and ''g'', but the drawback of a nested integral:
:int f g,dx = f int g,dx - int left ( f' int g,dx
ight ),dx.
This formula is valid whenever ''f'' is continuously differentiable and ''g'' is continuous.
More general formulations of integration by parts exist for the Riemann-Stieltjes integral and Lebesgue-Stieltjes integral.
Note: More complicated forms such as the one below are also valid:
:int u v,dw = u v w - int u w,dv - int v w,du.

Examples


In order to calculate:
:int xcos (x) ,dx
Let:
:''u'' = ''x'', so that ''du'' = ''dx'',
:''dv'' = cos(''x'') ''dx'', so that ''v'' = sin(''x'').
Then:
:
egin{align}
int xcos (x) ,dx & = int u ,dv \
& = uv - int v ,du \
& = xsin (x) - int sin (x) ,dx \
& = xsin (x) + cos (x) + C
end{align}

where ''C'' is an arbitrary constant of integration.
By repeatedly using integration by parts, integrals such as
:int x^{3} sin (x) ,dx quad mbox{and} quad int x^{2} e^{x} ,dx
can be computed in the same fashion: each application of the rule lowers the power of ''x'' by one.
An interesting example that is commonly seen is:
:int e^{x} cos (x) ,dx
where, strangely enough, in the end, the actual integration does not need to be performed.
This example uses integration by parts twice. First let:
:''u'' = cos(''x''); thus ''du'' = −sin(''x'') ''dx''
:d''v'' = e''x'' ''dx''; thus ''v'' = e''x''
Then:
:int e^{x} cos (x) ,dx = e^{x} cos (x) + int e^{x} sin (x) ,dx
Now, to evaluate the remaining integral, we use integration by parts again, with:
:''u'' = sin(''x''); d''u'' = cos(''x'') ''dx''
:''v'' = e''x''; ''dv'' = e''x'' ''dx''
Then:
:{|
|-
|int e^{x} sin (x) ,dx
|= e^{x} sin (x) - int e^{x} cos (x) ,dx
|-
|}
Putting these together, we get
:int e^{x} cos (x) ,dx = e^{x} cos (x) + e^x sin (x) - int e^{x} cos (x) ,dx
Notice that the same integral shows up on both sides of this equation. So we can simply add the integral to both sides to get:
:2 int e^{x} cos (x) ,dx = e^{x} ( sin (x) + cos (x) ) + C
:int e^{x} cos (x) ,dx = {e^{x} ( sin (x) + cos (x) ) over 2} + C
where, again, ''C'' is an arbitrary constant of integration.
A similar trick is used to find the integral of secant cubed.
Two other well-known examples are when integration by parts is applied to a function expressed as a product of 1 and itself. This works if the derivative of the function is known, and the integral of this derivative times ''x'' is also known.
The first example is ∫ ln(''x'') d''x''. We write this as:
:int ln (x) cdot 1 ,dx
Let:
:''u'' = ln(''x''); d''u'' = 1/''x'' d''x''
:''v'' = ''x''; d''v'' = 1·d''x''
Then:
:{|
|-
|int ln (x) ,dx
|= x ln (x) - int rac{x}{x} ,dx
|-
|
|= x ln (x) - int 1 ,dx
|}
:int ln (x) ,dx = x ln (x) - {x} + {C}
:int ln (x) ,dx = x ( ln (x) - 1 ) + C
where, again, C is the arbitrary constant of integration
The second example is ∫ arctan(''x'') d''x'', where arctan(''x'') is the inverse tangent function. Re-write this as:
:int rctan (x) cdot 1 ,dx
Now let:
:''u'' = arctan(''x''); d''u'' = 1/(1+''x''2) d''x''
:''v'' = ''x''; d''v'' = 1·d''x''
Then:
:{|
|-
|int rctan (x) ,dx
|= x rctan (x) - int rac{x}{1 + x^2} ,dx
|-
|
|= x rctan (x) - {1 over 2} ln left( 1 + x^2
ight) + C
|}
using a combination of the inverse chain rule method and the natural logarithm integral condition.

The ILATE rule


A rule of thumb for choosing which of two functions is to be ''u'' and which is to be ''dv'' is to choose ''u'' by whichever function comes first in this list:
:'I': inverse trigonometric functions: arctan ''x'', arcsec ''x'', etc.
:'L': logarithmic functions: ln ''x'', log_2(x), etc.
:'A': algebraic functions: x^2, 3x^{50}, etc.
:'T': trigonometric functions: sin ''x'', tan ''x'', etc.
:'E': exponential functions: e^x, 13^x, etc.
Then make ''dv'' the other function. You can remember the list by the mnemonic ILATE. The reason for this is that functions longer down in the list have easier antiderivatives than the functions above them.
To demonstrate this rule, consider the integral
:int xcos x ,dx.,
Following the ILATE rule, ''u'' = ''x'' and ''dv'' = cos ''x'' dx , hence ''du'' = ''dx'' and ''v'' = sin ''x'' , which makes the integral become
: xsin x - int 1sin x ,dx,,
which equals
: xsin x + cos x+C. ,
In general, one tries to choose ''u'' and ''dv'' such that ''du'' is simpler than ''u'' and ''dv'' is easy to integrate. If instead cos ''x'' was chosen as ''u'' and ''x'' as ''dv'', we would have the integral
: rac{x^2}2cos x + int rac{x^2}2sin x,dx,,
which, after recursive application of the integration by parts formula, would clearly result in an infinite recursion and lead nowhere.
Although a useful rule of thumb, there are exceptions to the ILATE rule. A common alternative is to consider the rules in the "LIATE" order instead. Also, in some cases, polynomial terms need to be split in non-trivial ways. For example, to integrate
:int x^3e^{x^2},dx,
we would set
:u=x^2,quad dv=xe^{x^2},dx.
This results in
:int x^3e^{x^2},dx= rac12e^{x^2}(x^2-1)+C.

Recursive integration by parts


Integration by parts can often be applied recursively on the int v,du term to provide the following formula
:int uv = u v_1 - u' v_2 + u'' v_3 - cdots + (-1)^{n} u^{(n)} v_{n+1}
Here, u' is the first derivative of u and u'' is the second derivative of u . Further, u^{(n)} is a notation to describe its ''n''th derivative (with respect to the variable ''u'' and ''v'' are functions of). Another notation has been adopted:
: v_{n+1}(x)=int! int cdots int v (dx)^{n+1}.
There are ''n'' + 1 integrals.
Note that the integrand above ( uv ) differs from the previous equation. The dv factor has been written as v purely for convenience.
The above mentioned form is convenient because it can be evaluated by differentiating the first term and integrating the second (with a sign reversal each time), starting out with u v_1 . It is very useful especially in cases when u^{(k+1)} becomes zero for some ''k'' + 1. Hence, the integral evaluation can stop once the u^{(k)} term has been reached.

Tabular integration by parts


While the aforementioned recursive definition is correct, it is often tedious to remember and implement. A much easier visual representation of this process is often taught to students and is dubbed either "the tabular method" or "the tic-tac-toe method". This method works when one of the two functions in the product is a polynomial, that is, after differentiating it several times one obtains zero.
For example, consider the integral
:int x^3 cos x ,dx.
Let u=x^3. Begin with this function and list in a column all the subsequent derivatives until zero is reached. Secondly, begin with the function ''v'' (in this case cos x) and list each integral of ''v'' until the size of the column is the same as that of ''u''. The result should appear as follows.
{| class="wikitable" style="text-align:center"
! Derivatives of ''u'' (Column A) !! Integrals of ''v'' (Column B)
|-
| x^3 , || cos x ,
|-
| 3x^2 , || sin x ,
|-
| 6x , || -cos x ,
|-
| 6 , || -sin x ,
|-
| 0 , || cos x ,
|}
Now simply pair the 1st entry of column A with the 2nd entry of column B, the 2nd entry of column A with the 3rd entry of column B, etc... with 'alternating' signs (beginning with the positive sign). Do so until further pairing is impossible. The result is the following (notice the alternating signs in each term):
:(+)(x^3)(sin x) + (-)(3x^2)(-cos x) + (+)(6x)(-sin x) + (-)(6)(cos x) + C ,.
Which, with simplification, leads to the result
:x^3sin x + 3x^2cos x - 6xsin x - 6cos x + C. ,

Higher dimensions


The formula for integration by parts can be extended to functions of several variables. Instead of an interval one needs to integrate over a ''n''-dimensional set. Also, one replaces the derivative with a partial derivative.
More specifically, suppose Ω is an open bounded subset of mathbb{R}^n with a piecewise smooth boundary ∂Ω. If ''u'' and ''v'' are two continuously differentiable functions on the closure of Ω, then the formula for integration by parts is
:int_{Omega} rac{partial u}{partial x_i} v ,dx = int_{partialOmega} u v ,
u_i ,dsigma - int_{Omega} u rac{partial v}{partial x_i} , dx
where mathbf{
u} is the outward unit surface normal to ∂Ω, ν''i'' is its ''i''-th component, and ''i'' ranges from 1 to ''n''. Replacing ''v'' in the above formula with ''v''''i'' and summing over ''i'' gives the vector formula
: int_{Omega}
abla u cdot mathbf{v}, dx = int_{partialOmega} u, mathbf{v}cdot
u, dsigma - int_Omega u,
ablacdot mathbf{v}, dx
where 'v' is a vector-valued function with components ''v''1, ..., ''v''''n''.
Setting ''u'' equal to the constant function 1 in the above formula gives the divergence theorem. For mathbf{v}=
abla v where vin C^2(ar{Omega}), one gets
: int_{Omega}
abla u cdot
abla v, dx = int_{partialOmega} u,
abla vcdot
u, dsigma - int_Omega u, Delta v, dx
which is the first Green's identity.
The regularity requirements of the theorem can be relaxed. For instance, the boundary ∂Ω need only be Lipschitz continuous. In the first formula above, only u,vin H^1(Omega) is necessary (where ''H''1 is a Sobolev space); the other formulas have similarly relaxed requirements.
For reference, consult Appendix C of Evans or the applied math notes of Arbogast and Bona.

Cultural references



★ The method of tabular integration by parts is featured in the 1988 film ''Stand and Deliver''.

References



Partial Differential Equations, , Lawrence C., Evans, American Mathematical Society, 1998, ISBN 0-8218-0772-2

Methods of Applied Mathematics, , Todd, Arbogast, , 2005,

Tabular Integration by Parts, , David, Horowitz, The College Mathematics Journal,

External links



Integration by Parts - From MathWorld

Tabular Integration by Parts

Tabular Integration by Parts Demonstrated

This article provided by Wikipedia. To edit the contents of this article, click here for original source.

psst.. try this: add to faves