guide:266cbfc1ad: Difference between revisions

From Stochiki
No edit summary
 
No edit summary
Line 1: Line 1:
<div class="d-none"><math>
\newcommand{\ex}[1]{\item }
\newcommand{\sx}{\item}
\newcommand{\x}{\sx}
\newcommand{\sxlab}[1]{}
\newcommand{\xlab}{\sxlab}
\newcommand{\prov}[1] {\quad #1}
\newcommand{\provx}[1] {\quad \mbox{#1}}
\newcommand{\intext}[1]{\quad \mbox{#1} \quad}
\newcommand{\R}{\mathrm{\bf R}}
\newcommand{\Q}{\mathrm{\bf Q}}
\newcommand{\Z}{\mathrm{\bf Z}}
\newcommand{\C}{\mathrm{\bf C}}
\newcommand{\dt}{\textbf}
\newcommand{\goesto}{\rightarrow}
\newcommand{\ddxof}[1]{\frac{d #1}{d x}}
\newcommand{\ddx}{\frac{d}{dx}}
\newcommand{\ddt}{\frac{d}{dt}}
\newcommand{\dydx}{\ddxof y}
\newcommand{\nxder}[3]{\frac{d^{#1}{#2}}{d{#3}^{#1}}}
\newcommand{\deriv}[2]{\frac{d^{#1}{#2}}{dx^{#1}}}
\newcommand{\dist}{\mathrm{distance}}
\newcommand{\arccot}{\mathrm{arccot\:}}
\newcommand{\arccsc}{\mathrm{arccsc\:}}
\newcommand{\arcsec}{\mathrm{arcsec\:}}
\newcommand{\arctanh}{\mathrm{arctanh\:}}
\newcommand{\arcsinh}{\mathrm{arcsinh\:}}
\newcommand{\arccosh}{\mathrm{arccosh\:}}
\newcommand{\sech}{\mathrm{sech\:}}
\newcommand{\csch}{\mathrm{csch\:}}
\newcommand{\conj}[1]{\overline{#1}}
\newcommand{\mathds}{\mathbb}
</math></div>
===Partial Fractions.===
A rational function is by definition one which can be expressed as the ratio
of two polynomials. A simple example is the function <math>f</math>, defined by


<math display="block">
f(x) = \frac{1}{(x^2 + 1)(x - 2)} =  \frac{1}{x^3 - 2x^2 + x - 2},
</math>
for every real value of <math>x</math> except 2. At present, we have no way of integrating this function. However, in this section we shall develop a method of integration which is applicable to any rational function. It is called the method of partial fractions.
To illustrate the method, consider the equation
<math display="block">
\frac{1}{x - 2} - \frac{x + 2}{x^2 + 1} = \frac{5}{(x^2 + 1)(x - 2)},
</math>
which is easily seen to be true for all real values of <math>x</math>
except 2. It follows that
<math display="block">
\begin{eqnarray*}
\frac{1}{(x^2 + 1)(x - 2)}
&=& \frac{1}{5} \frac{1}{x - 2} - \frac{1}{5} \frac{x + 2}{x^2 + 1}\\
&=& \frac{1}{5} \frac{1}{x - 2} - \frac{1}{5} \frac{x}{x^2 + 1} - \frac{2}{5} \frac{1}{x^2 + 1}.
\end{eqnarray*}
</math>
Hence
<math display="block">
\begin{eqnarray*}
\int \frac{dx}{(x^2 + 1)(x - 2)}
&=& \frac{1}{5} \int \frac{dx}{x - 2} - \frac{1}{5} \int \frac{x dx}{x^2 + 1} - \frac{2}{5} \int \frac{dx}{x^2 + 1}\\
&=& \frac{1}{5} \ln |x - 2| - \frac{1}{10} \ln(x^2 + 1) - \frac{2}{5} \arctan x + c. 
\end{eqnarray*}
</math>
Thus <math>\frac{1}{(x^2 + 1)(x -2)}</math> can be integrated, since it can be written as a sum of simpler rational functions, each of which can be integrated separately. The method of integration by partial fractions is based on the fact that such a decomposition exists for every rational function. In this example, we have given no indication of how the decomposition is to be found. However, the general method, which we now describe, consists of just such a prescription.
We begin with the result from algebra that it is always possible by means
of division to express any given rational function as the sum of a polynomial
and a rational function in which the degree of the numerator is less than the degree of the denominator. Stated formally, this theorem says that, if <math>N(x)</math> and <math>D(x)</math> are any two polynomials and if <math>D(x)</math> is not the zero function, then there exist uniquely determined polynomials <math>Q(x)</math> and <math>R(x)</math> such that
<span id{{=}}"eq7.4.1"/>
<math display="block">
\begin{equation}
\frac{N(x)}{D(x)} = Q(x) + \frac{R(x)}{D(x)},
\label{eq7.4.1}
\end{equation}
</math>
and such that the degree of <math>R(x)</math> is less than the degree of <math>D(x)</math>. (The letters <math>N</math>, <math>D</math>, <math>Q</math>, and <math>R</math> have been used to suggest, respectively, the words “numerator,” “denominator,” “quotient,” and
“remainder.”) The first step in the method of integration by partial fractions is to write the given rational function <math>\frac{N(x)}{D(x)}</math> in the form of equation (1). Since we can obviously integrate the polynomial <math>Q(x)</math>, we need next consider only rational functions in which the degree of the numerator is less than the degree of the denominator. If we start with such a function, then no division is necessary.
'''Example'''
Write the function <math>\frac{x^4 - 4x^3 + 8x^2 - 7x +3}{x^3 - 2x^2 + 3x - 4}</math> as the sum of a polynomial and a rational function in which the degree of the numerator is less than the degree of the denominator. Dividing, we have
<math display="block">
\begin{array}{rl}
                              &\;\;x - 2  \\     
x^3 - 2x^2 + 3x - 4&)\;\overline {x^4 - 4x^3 + 8x^2 - 7x + 3}  \\
                            &x^4 - 2x^3 + 3x^2 - 4x \\   
                            &\overline{\;\;\;\;\;\;\;\, -2x^3 + 5x^2 - 3x + 3} \\ 
                              &\underline{-2x^3 + 4x^2 - 6x + 8} \\ 
                              &x^2 + 3x - 5.
\end{array}
</math>
It follows that
<math display="block">
\frac{x^4 - 4x^3 + 8x^2 - 7x + 3}{x^3 - 2x^2 + 3x - 4}
= (x - 2) + \frac{x^2 + 3x - 5}{x^3 - 2x^2 + 3x - 4},
</math>
which gives the required sum.
Another algebraic fact about polynomials, which we shall not prove, but shall assume and use, is that any nonconstant polynomial (i.e., of degree at least 1) with real coefficients can be written as a product of linear and quad
ratic factors, each with real coefficients. By a linear factor we mean a polynomial <math>L(x)</math> of degree 1; that is, <math>L(x) = ax + b</math> and <math>a \neq 0</math>. Similarly, a quadratic factor is a polynomial <math>Q(x)</math> of degree 2; thus
<math>Q(x) = cx^2 + dx + e</math> and <math>c \neq 0</math>. The theorem states that, for any polynomial
<math display="block">
f(x) = a_{n}x^n + ... + a_{1}x + a_{0}
</math>
with real coefficients <math>a_i</math> and with <math>n \geq 1</math> and <math>a_n \neq 0</math>, there exist linear factors <math>L_{1}(x), \cdots, L_{p}(x)</math> and quadratic factors <math>Q_1(x), \cdots, Q_q(x)</math> with real coefficients such that
<math display="block">
f(x)= L_1(x) \cdots L_p(x)Q_1(x) \cdots Q_q(x).
</math>
Note that either <math>p</math> or <math>q</math> may be zero. In actual practice, such a factorization of <math>f(x)</math> may be very difficult to find, but the theorem assures us that it exists.
A polynomial is said to be '''irreducible''' if it cannot be written as the product of two polynomials each of degree greater than or equal to 1. The degree of the product of two polynomials is equal to the sum of the degrees of the factors, and it therefore follows that every linear polynomial is irreducible. It was pointed out in Section 3 that a quadratic polynomial <math>cx^2 + dx + e</math> is irreducible over the reals if and only if its discriminant <math>d^2 - 4ce</math> is negative. For example, the polynomials <math>x^2 + 1</math> and <math>x^2 + x + 1</math> are irreducible, whereas <math>x^2 + 2x + 1</math> and <math>x^2 + 2x - 1</math> are not. If a quadratic polynomial is not irreducible, it can be factored and written as the product of two linear polynomials. Hence the factorization of an arbitrary nonconstant polynomial into linear and quadratic factors, as described in the preceding paragraph, can always be done so that all the factors are irreducible.
Returning specifically to the method of integration by partial fractions, we consider a rational function <math>\frac{N(x)}{D(x)}</math>, with the degree of <math>D(x)</math> greater than the degree of <math>N(x)</math>. The second step is to write the denominator <math>D(x)</math> as a product of irreducible factors. Having done so, we have
<span id{{=}}"eq7.4.2"/>
<math display="block">
\begin{equation}
D(x) = L_1(x) \cdots L_p(x)Q_1(x) \cdots Q_q(x), 
\label{eq7.4.2}
\end{equation}
</math>
where, for each <math>i = 1, . . ., p</math>,
<math display="block">
L_i(x) = a_{i}x + b_{i}, \;\;\; a_i \neq 0,
</math>
and, for each <math>j = 1, . . ., q</math>,
<math display="block">
Q_j(x) = c_{j}x^2 + d_{j} x + e_{j}, \;\;\; d_{j}^2 - 4c_{j}e_{j}  <  0.
</math>
There is no reason to suppose that the factors which appear in equation (2) will all be distinct, and it may very well happen that
<math>L_{1}(x) = L_{2}(x)</math>, etc.
However, the theory is simpler if no repetitions occur, and we shall consider that case first.
''Case 1. The irreducible factors of $D(x)$ are all distinct.'' 
The algebraic theory of partial fractions, which we shall assume, tells us that we can write <math>\frac{N(x)}{D(x)}</math> as a sum of rational functions each of which has one of the factors of <math>D(x)</math> as its denominator and such that in each term the degree of the numerator is less than the degree of the denominator. Moreover, given the factorization of <math>D(x)</math> into irreducibles, this decomposition is unique except for the order in which the terms are written. Thus
<math display="block">
\begin{eqnarray*}
\frac{N(x)}{D(x)} &=& \frac{A_1}{a_{1}x + b_{1}} + \frac{A_2}{a_{2}x + b_{2}} + \cdots + \frac{A_p}{a_{p}x + b_{p}}\\
&& + \frac{B_{1}x + C_1} {c_{1}x^2 + d_{1}x + e_{1}} + \frac{B_{2}x + C_2}{c_{2}x^2 + d_{2}x + e_2}\\
&& + \cdots
+ \frac{B_{q}x + C_q}{c_{q}x^2 + d_{q}x + e_q},
\end{eqnarray*}
</math>
where each of the letters <math>A_i</math>, <math>B_j</math> and <math>C_k</math> represents a uniquely determined real constant. The rational functions which appear on the right side are called the '''partial fractions''' of the decomposition of <math>\frac{N(x)}{D(x)}</math>.
We shall show by means of examples how the constants in the partial fractions decomposition are determined. Consider the rational function
<math>\frac{1}{(x^2 + 1)(x - 2)}</math> discussed at the beginning of the section. The degree of the numerator, zero, is already less than 3, the degree of the denominator. Moreover, the denominator is already factored into irreducibles. Hence, we seek constants <math>A</math>, <math>B</math>, and <math>C</math> such that
<math display="block">
\frac{1}{(x^2 + 1)(x - 2)} = \frac{A}{x - 2} + \frac{Bx + C}{x^2 + 1}.
</math>
Adding the two fractions on the right side, we have
<math display="block">
\frac{1}{(x^2 + 1)(x - 2)} = \frac{A(x^2 + 1) + (Bx + C)(x - 2)}{(x^2 + 1)(x - 2)}.
</math>
The fact that a nonzero polynomial of degree <math>n</math> has at most <math>n</math> distinct roots implies that two rational functions with the same denominator are equal if and only if their numerators are equal (see Problem 8 at the end of this section). Hence the equation
<span id{{=}}"eq7.4.3"/>
<math display="block">
\begin{equation}
1 = A(x^2 + 1) + (Bx + C)(x - 2)
\label{eq7.4.3}
\end{equation}
</math>
is true for all real values of <math>x</math>. There are two common ways to find the constants <math>A</math>, <math>B</math>, and <math>C</math>. One is to multiply and regroup the terms in (3) to obtain the equation
<math display="block">
0 = (A + B)x^2 + (C - 2B)x + (A - 2C - 1),
</math>
which also holds for all real values of <math>x</math>. The only polynomial with infinitely many roots is the zero polynomial, i.e., the polynomial with no nonzero coefficients. Hence the three coefficients on the right side of the preceding equation are all equal to zero, and we can therefore find <math>A</math>, <math>B</math>, and <math>C</math> by solving the system of equations
<math display="block">
\begin{array}{ccrcrc}
A  &  + &    B &    &    & = 0,\\
  &    - &  2B & +  &  C & = 0,\\
A  &      &      &  - & 2C & = 1.
\end{array}
</math>
Usually simpler is the technique in which we take advantage of the fact that (3) must be true for all values of <math>x</math>, and we choose values cleverly to help evaluate the constants. For example, letting <math>x = 2</math> in (3), we have
<math display="block">
\begin{eqnarray*}
1 &=& A(2^2 + 1) + (B \cdot 2 + C)(2 - 2)\\
  &=& 5A.
\end{eqnarray*}
</math>
Hence <math>A = \frac{1}{5}</math>. lf we then let <math>x = 0</math>, we have 
<math display="block">
\begin{eqnarray*}
1 &=& \frac{1}{5}(0^2 + 1) + (B \cdot 0 + C)(0 - 2) \\
  &=& \frac{1}{5} - 2C,
\end{eqnarray*}
</math>
or, equivalently, <math>2C = \frac{1}{5} - 1</math>, and so <math>C = - \frac{2}{5}</math>. Finally, choosing <math>x = 1</math>, we get
<math display="block">
\begin{eqnarray*}
1 &=&\frac{1}{5}(1^2 + 1) + (B \cdot 1 - \frac{2}{5})(1 - 2)\\
  &=& \frac{2}{5} - B + \frac{2}{5},
\end{eqnarray*}
</math>
from which we conclude that <math>B = - \frac{1}{5}</math>. Whichever method we use, we have <math>A = \frac{1}{5}</math>, <math>B = - \frac{1}{5}</math>, and <math>C = -\frac{2}{5}</math>, from which it follows that
<math display="block">
\frac{1}{(x^2+ 1)(x - 2)}
= \frac{1}{ 5}\frac{1}{ x-2} - \frac{1}{ 5} \frac{x + 2}{x^2+ 1},
</math>
the form which we integrated at the beginning of the section.
'''Example'''
Integrate <math>\int \frac{dx}{a^2 - x^2}</math>. Since <math>a^2 - x^2 = (a - x)(a + x)</math>, we
decompose <math>\frac{1}{a^2 - x^2}</math> into partial fractions <math>\frac{A}{a + x}</math> and <math>\frac{B}{a - x}</math>. Thus
<math display="block">
\frac{1}{a^2 - x^2} =  \frac{A}{a + x} + \frac{B}{a - x}
= \frac{A(a - x) + B(a + x)}{a^2 - x^2}.
</math>
Equating numerators on the left and right, we get
<math display="block">
1 = A(a - x) + B(a + x).
</math>
Letting <math>x = a</math>, we obtain the equation <math>1 = A \cdot 0 + B \cdot 2a</math>, and so <math>B = \frac{1}{2a}</math>. Similarly, setting <math>x = - a</math>, we get <math>1 = A \cdot 2a
+ B \cdot 0</math>, from which it follows that <math>A = \frac{1}{2a}</math>. Thus 
<math display="block">
\frac{1}{a^2 - x^2} = \frac{1}{2a} \frac{1}{a + x} + \frac{1}{2a} \frac{1}{a - x},
</math>
and, therefore,
<math display="block">
\begin{eqnarray*}
\int \frac{ dx}{a^2 - x^2}
&=& \frac{1}{2a} \int \frac{dx}{a + x} + \frac{1}{2a} \int \frac{dx}{a - x}\\
&=& \frac{1}{2a} \ln |a + x| - \frac{1}{2a} \ln |a - x| + c \\
&=& \frac{1}{2a} \ln \Big| \frac{a + x}{a - x} \Big|  + c.
\end{eqnarray*}
</math>
Thus the third step in applying this method of integration is the decomposition into partial fractions, and the fourth and final step is the integration of the partial fractions. We shall show later in the section that it is always possible to carry out the last step, but, as the next example shows, doing so can be tedious.
<span id="eq7.4.4"/>
'''Example'''
Integrate <math>\int \frac{3x^3 + x^2 -14x + 46}{(x^2 + x + 1)(x^2 - 5x - 14)} dx</math>. The degree of the numerator is 3 and that of the denominator is 4, so we proceed to the factorization of the denominator into irreducibles. It is already written as the product of two quadratics, of which <math>x^2 + x + 1</math> is irreducible but <math>x^2 - 5x -14</math> is not, since <math>x^2 - 5x - 14 = (x - 7)(x + 2)</math>. Hence the form of the partial fractions decomposition is
<math display="block">
\frac{3x^3 + x^2 -14x + 46}{(x^2 + x + 1)(x - 7)(x + 2)}
= \frac{Ax + B}{x^2 + x + 1} + \frac{C}{x - 7} + \frac{D}{x + 2}.
</math>
The sum of the three fractions on the right side is
<math display="block">
\begin{array}{l}
(Ax + B)(x - 7)(x + 2)\\
\underline{+ C(x^2 + x + 1)(x + 2) + D(x^2 + x + 1)(x - 7)} \\
(x^2 + x + 1)(x - 7)(x + 2)
\end{array}.
</math>
Equating numerators, we have
<math display="block">
\begin{eqnarray*}
&&(Ax + B)(x-7)(x+ 2) + C(x^2 + x + 1)(x + 2) + D(x^2 + x + 1)(x-7) \\
&& = 3x^3 + x^2 - 14x + 46.
\end{eqnarray*}
</math>
If we set <math>x = 7</math> in this equation, then
<math display="block">
C \cdot 513= 1026 \;\;\;\mbox{or}\;\;\; C = 2.
</math>
Letting <math>x = - 2</math>, we obtain
<math display="block">
D \cdot (- 27) = 54 \;\;\; \mbox{or} \;\;\; D = - 2.
</math>
If we let <math>x = 0</math>, then
<math display="block">
B \cdot (- 14) + C \cdot 2 + D \cdot ( - 7) = 46
</math>
or
<math display="block">
-14B = 46 - 2C + 7D = 46 - 4 - 14 = 28 \;\;\;\mbox{or} \;\;\; B = - 2.
</math>
Finally, letting <math>x = - 1</math>, we have
<math display="block">
(-A + B) \cdot (-8) + C + D \cdot (-8) = 58
</math>
or
<math display="block">
8A = 58 + 8B - C + 8D = 58 - 16 - 2 - 16 = 24 \;\;\; \mbox{or} \;\;\; A = 3.
</math>
The partial fractions decomposition is, therefore,
<span id{{=}}"eq7.4.4"/>
<math display="block">
\begin{equation}
\frac{3x^3 + x^2 - 14x + 46}{(x^2 + x + 1)(x^2 - 5x -14)} = \frac{3x - 2}{x^2 + x + 1} + \frac{2}{x - 7} - \frac{2}{x + 2}.
\label{eq7.4.4}
\end{equation}
</math>
Except for the first, the terms on the right side are easily integrated. The first term can be integrated by writing it as the sum of two fractions. We use the identity
<math display="block">
\frac{Bx + C}{cx^2+ dx+e} = \frac{B}{2c} \frac{2cx + d}{cx^2+ dx + e}
+ \Bigl(C - \frac{dB}{2c} \Bigr) \frac{1}{cx^2+ dx + e}.
</math>
Note that the numerator <math>2cx + d</math> is the derivative of <math>cx^2 + dx + e</math>. Both of these fractions can be integrated. In this example, we have
<math display="block">
\frac{3x - 2}{x^2 + x + 1} = \frac{3}{2} \frac{2x + 1}{x^2 + x + 1}
- \frac{7}{2} \frac{1}{x^2 + x + 1}.
</math>
Then
<math display="block">
\int \frac{2x + 1}{x^2 + x + 1} dx = \ln (x^2 + x + 1) + c.
</math>
For the second fraction, we complete the square in the denominator. The result is
<math display="block">
\frac{1}{x^2 + x + 1} = \frac{1}{\Bigl(x + \frac{1}{2}\Bigr)^2 + \Bigl(\frac{\sqrt 3}{2}\Bigr)^2},
</math>
and, since <math>\int \frac{dy}{y^2 + a^2} = \frac{1}{|a|} \arctan \frac{y}{|a|}</math>, it follows that
<math display="block">
\begin{eqnarray*}
\int \frac{dx}{x^2 + x + 1} = \int \frac{dx}{\Bigl(x + \frac{1}{2}\Bigr)^2 + \Bigl(\frac{\sqrt 3}{2}\Bigr)^2}
&=& \frac{2}{\sqrt 3} \arctan \Bigl(\frac{x + \frac{1}{2}}{\frac{\sqrt 3}{2}} \Bigr) + c\\
&=& \frac{2}{\sqrt 3} \arctan \Bigl(\frac{2x + 1}{\sqrt 3} \Bigr) + c.
\end{eqnarray*}
</math>
Hence
<math display="block">
\int \frac{3x - 2}{x^2 + x + 1} dx = \frac{3}{2} \ln(x^2 + x + 1) - \frac{7}{\sqrt 3} \arctan \Bigl(\frac{2x + 1}{\sqrt 3} \Bigr) + c.
</math>
Returning to equation (4), we therefore get the final integral
<math display="block">
\begin{eqnarray*}
\int \frac{3x^3 + x^2 -14x + 46}{(x^2 + x + 1)(x^2 - 5x - 14)} dx
&=& \frac{3}{2} \ln(x^2 + x + 1) - \frac{7}{\sqrt 3} \arctan \Bigl(\frac{2x+ 1}{\sqrt 3} \Bigr)\\
&& + 2 \ln |x - 7| - 2 \ln |x + 2| + c,
\end{eqnarray*}
</math>
and this completes the example.
We consider next the situation in which the factorization of the denominator of <math>\frac{N(x)}{D(x)}</math> as shown in equation (2), contains repeated factors.
\medskip
''Case 2. The irreducible factors of $D(x)$ are not all distinct.''
We assume, as in Case 1, that the degree of <math>N(x)</math> is less than the degree of <math>D(x)</math>. There is still a unique decomposition of <math>\frac{N(x)}{D(x)}</math> into the sum of partial fractions, but now it is more complicated. By regrouping, we may write the factorization of <math>D(x)</math> into irreducibles as
<span id{{=}}"eq7.4.5"/>
<math display="block">
\begin{equation}
D(x) = [L_{1}(x)]^{m_1} \cdots  [L_{r}(x)]^{m_r} [Q1(x)]^{n_1} \cdots [Q_{s}(x)]^{n_s},
\label{eq7.4.5}
\end{equation}
</math>
where <math>m_1, ... , m_r</math> and <math>n_1, ... , n_s</math> are positive integers,
the factors <math>L_{i}(x) = a_{i}x + b_{i}</math> are all distinct, and the factors
<math>Q_{j}(x) = c_{j}x^2 + d_{j}x + e_{j}</math> are all distinct In this case, <math>\frac{N(x)}{D(x)}</math> is the total sum of the following individual sums
of partial fractions: For each <math>i = 1, . . ., r</math>, there is the sum
<math display="block">
\frac{A_{i1}}{a_{i}x + b_{i}} + \frac{A_{i2}}{(a_{i}x + b_{i})^2} + \cdots + \frac{A_{im_i}}{(a_{i}x + b_{i})^{m_i}}, 
</math>
in which the <math>Aik</math> are uniquely determined real constants. Similarly, for each <math>j = 1, ... , s</math>, there is the sum
<math display="block">
\frac{B_{j1}x + C_{j1}}{c_{j}x^2+ d_{j}x + e_{j}} 
+ \frac{B_{j2}x + C_{j2}}{(c_{j}x^2 + d_{j}x + e_{j})^2}
+ \cdots
+ \frac{B_{jn_j}x + C_{jn_j}} {(c_{j}x^2 + d_{j}x + e_{j})^{n_j}},
</math>
in which the <math>B_{jk}</math> and <math>C_{jk}</math> are uniquely
determined real constants.
<span id="eq7.4.6"/>
'''Example'''
Integrate the rational function <math>\frac{2x^2 + x + 2}{x(x - 1)^3}</math> by the method of partial fractions. The degree of the numerator, 2, is less than that of the denominator, 4. So we turn at once to the decomposition into partial fractions. Since the irreducible factor <math>x - 1</math> is repeated three times, the decomposition is of the form
<math display="block">
\frac{2x^2 + x + 2}{x(x - 1)^3}
= \frac{A}{x - 1} + \frac{B}{(x - 1)^2} + \frac{C}{(x - 1)^3} + \frac{D}{x}.
</math>
The sum on the right side is equal to
<math display="block">
\frac{Ax(x - 1)^2 + Bx(x - 1) + Cx + D(x -1)^3}{x(x - 1)^3}.
</math>
Equating numerators, we obtain the equation
<span id{{=}}"eq7.4.6"/>
<math display="block">
\begin{equation}
Ax(x - 1)^2 + Bx(x - 1) + Cx + D(x - 1)^3 = 2x^2 + x + 2,
\label{eq7.4.6}
\end{equation}
</math>
which holds for all real values of <math>x</math>. Setting <math>x = 1</math>, we obtain
<math display="block">
C \cdot 1 = 2 \cdot 1^2 + 1 + 2, \;\;\;\mbox{whence}\; C = 5.
</math>
If we let <math>x = 0</math>, then
<math display="block">
D \cdot ( - 1)^3 = 2, \;\;\; \mbox{whence}\; D = - 2.
</math>
Thus, equation (6) has become
<math display="block">
Ax(x - 1)^2 + Bx(x - 1) + 5x -2(x - 1)^3 = 2x^2 + x + 2.
</math>
In this equation we let <math>x = 2</math>, getting
<math display="block">
A \cdot 2 + B \cdot 2 + 10 - 2 = 8 + 2 + 2
</math>
or
<span id{{=}}"eq7.4.7"/>
<math display="block">
\begin{equation}
2A + 2B = 4. 
\label{eq7.4.7}
\end{equation}
</math>
Next, setting <math>x = - 1</math>, we have
<math display="block">
A \cdot (- 1)(- 2)^2 + B \cdot (- 1)(- 2) - 5 + 16 = 2 - 1 + 2
</math>
or
<span id{{=}}"eq7.4.8"/>
<math display="block">
\begin{equation}
-4A + 2B = - 8.
\label{eq7.4.8}
\end{equation}
</math>
Subtracting (8) from (7), we get <math>6A = 12</math> and so <math>A = 2</math>. It follows that <math>B = 0</math>, and we have therefore found the partial fractions decomposition to be
<math display="block">
\frac{2x^2 + x + 2}{x(x - 1)^3} = \frac{2}{x - 1} + \frac{5}{(x - 1)^3} - \frac{2}{x}.
</math>
Hence
<math display="block">
\begin{eqnarray*}
\int \frac{2x^2 + x + 2}{x(x -1)^3} dx
&=& 2\int \frac{dx}{x - 1} + 5 \int \frac{dx}{(x - 1)^3} - 2 \int \frac{dx}{x}\\
&=& 2 \ln | x - 1| - \frac{5}{2} \frac{1}{(x - 1)^2} - 2 \ln |x| + c,
\end{eqnarray*}
</math>
and this completes the example.
Since any rational function can be written as the sum of a polynomial and a series of partial fractions, the general problem of integrating a rational function reduces to three integration problems: (1) integration of a polynomial; (2) integration of functions of the form <math>\frac{A}{(ax + b)^m}</math>, where <math>m</math> is a positive integer and <math>a \neq 0</math>; and (3) integration of functions of the form <math>\frac{Bx + C}{(cx^2 + dx + e)^n}</math> where <math>n</math> is a positive integer and <math>d^2 - 4ce  <  0</math>. The first, of course, offers no difficulties whatever. The second is also simple, since <math>\frac{A}{(ax + b)^m} = \frac{A}{a} \frac{a}
{(ax + b)^m}</math>, and so 
<math display="block">
\int \frac{A}{(ax + b)^m} dx = \frac{A}{a} \int \frac{a}{(ax + b)^m}dx\\
= \left \{
              \begin{array}{ll}
\frac{A}{a(1 - m)} \frac{1}{(ax + b)^{m-1}} + c  &\mbox{if}\; m \neq 1,\\
\frac{A}{a} \ln | ax + b | + c                              &\mbox{if}\; m = 1.
              \end{array}
\right.
</math>
The third problem can be solved, but we have seen in Example 3 that it is complicated even with <math>n = 1</math>. It is attacked by writing
<math display="block">
\frac{Bx + C} {(cx^2 + dx + e)^n} = \frac{B}{2c} \frac{2cx + d}{(cx^2 + dx + e)^n} + \Bigl(C - \frac{dB}{2c} \Bigr) \frac{1}{ (cx^2 + dx +
e)^n}.
</math>
The integra <math>\int \frac{2cx + d}{(cx^2 + dx + e)^n} dx</math> is easily found, since it is of the form <math>\int \frac{du}{u^n}</math> with <math>u = cx^2 + dx + e</math>.
The problem therefore reduces to finding <math>\int \frac{dx}{(cx^2 + dx + e)^n}</math>. However it was demonstrated at the end of Section 3 that this integral can be evaluated by trigonometric substitution, which reduces it to an integral of the form <math>\int \cos^{2n-2} \theta d\theta</math>. An alternative method is to use the reduction formula given in Problem 9, page 384.
Thus all possible partial fractions resulting from the decomposition of a rational function can be integrated. It follows that every rational function can be integrated. The factoring of the denominator into irreducibles may be difficult, and the decomposition into partial fractions and the resulting integrations may be tedious, but the following important result has been established.
{{proofcard|Theorem|theorem-1|Every rational function can be integrated by the method of partialfractions.|}}
A table of integrals will show how to integrate many of the functions which are the partial fractions of a rational fraction. For this reason, there is no need to memorize the formulas for integration. However, it is necessary to know the technique of separating a rational function into its partial fractions in order to replace an apparently nonintegrable function by a sum of obviously integrable functions.
\end{exercise}
==General references==
{{cite web |title=Crowell and Slesnick’s Calculus with Analytic Geometry|url=https://math.dartmouth.edu/~doyle/docs/calc/calc.pdf |last=Doyle |first=Peter G.|date=2008 |access-date=Oct 29, 2024}}

Revision as of 00:08, 3 November 2024

[math] \newcommand{\ex}[1]{\item } \newcommand{\sx}{\item} \newcommand{\x}{\sx} \newcommand{\sxlab}[1]{} \newcommand{\xlab}{\sxlab} \newcommand{\prov}[1] {\quad #1} \newcommand{\provx}[1] {\quad \mbox{#1}} \newcommand{\intext}[1]{\quad \mbox{#1} \quad} \newcommand{\R}{\mathrm{\bf R}} \newcommand{\Q}{\mathrm{\bf Q}} \newcommand{\Z}{\mathrm{\bf Z}} \newcommand{\C}{\mathrm{\bf C}} \newcommand{\dt}{\textbf} \newcommand{\goesto}{\rightarrow} \newcommand{\ddxof}[1]{\frac{d #1}{d x}} \newcommand{\ddx}{\frac{d}{dx}} \newcommand{\ddt}{\frac{d}{dt}} \newcommand{\dydx}{\ddxof y} \newcommand{\nxder}[3]{\frac{d^{#1}{#2}}{d{#3}^{#1}}} \newcommand{\deriv}[2]{\frac{d^{#1}{#2}}{dx^{#1}}} \newcommand{\dist}{\mathrm{distance}} \newcommand{\arccot}{\mathrm{arccot\:}} \newcommand{\arccsc}{\mathrm{arccsc\:}} \newcommand{\arcsec}{\mathrm{arcsec\:}} \newcommand{\arctanh}{\mathrm{arctanh\:}} \newcommand{\arcsinh}{\mathrm{arcsinh\:}} \newcommand{\arccosh}{\mathrm{arccosh\:}} \newcommand{\sech}{\mathrm{sech\:}} \newcommand{\csch}{\mathrm{csch\:}} \newcommand{\conj}[1]{\overline{#1}} \newcommand{\mathds}{\mathbb} [/math]

Partial Fractions.

A rational function is by definition one which can be expressed as the ratio of two polynomials. A simple example is the function [math]f[/math], defined by

[[math]] f(x) = \frac{1}{(x^2 + 1)(x - 2)} = \frac{1}{x^3 - 2x^2 + x - 2}, [[/math]]

for every real value of [math]x[/math] except 2. At present, we have no way of integrating this function. However, in this section we shall develop a method of integration which is applicable to any rational function. It is called the method of partial fractions. To illustrate the method, consider the equation

[[math]] \frac{1}{x - 2} - \frac{x + 2}{x^2 + 1} = \frac{5}{(x^2 + 1)(x - 2)}, [[/math]]

which is easily seen to be true for all real values of [math]x[/math] except 2. It follows that

[[math]] \begin{eqnarray*} \frac{1}{(x^2 + 1)(x - 2)} &=& \frac{1}{5} \frac{1}{x - 2} - \frac{1}{5} \frac{x + 2}{x^2 + 1}\\ &=& \frac{1}{5} \frac{1}{x - 2} - \frac{1}{5} \frac{x}{x^2 + 1} - \frac{2}{5} \frac{1}{x^2 + 1}. \end{eqnarray*} [[/math]]


Hence

[[math]] \begin{eqnarray*} \int \frac{dx}{(x^2 + 1)(x - 2)} &=& \frac{1}{5} \int \frac{dx}{x - 2} - \frac{1}{5} \int \frac{x dx}{x^2 + 1} - \frac{2}{5} \int \frac{dx}{x^2 + 1}\\ &=& \frac{1}{5} \ln |x - 2| - \frac{1}{10} \ln(x^2 + 1) - \frac{2}{5} \arctan x + c. \end{eqnarray*} [[/math]]


Thus [math]\frac{1}{(x^2 + 1)(x -2)}[/math] can be integrated, since it can be written as a sum of simpler rational functions, each of which can be integrated separately. The method of integration by partial fractions is based on the fact that such a decomposition exists for every rational function. In this example, we have given no indication of how the decomposition is to be found. However, the general method, which we now describe, consists of just such a prescription. We begin with the result from algebra that it is always possible by means of division to express any given rational function as the sum of a polynomial and a rational function in which the degree of the numerator is less than the degree of the denominator. Stated formally, this theorem says that, if [math]N(x)[/math] and [math]D(x)[/math] are any two polynomials and if [math]D(x)[/math] is not the zero function, then there exist uniquely determined polynomials [math]Q(x)[/math] and [math]R(x)[/math] such that

[[math]] \begin{equation} \frac{N(x)}{D(x)} = Q(x) + \frac{R(x)}{D(x)}, \label{eq7.4.1} \end{equation} [[/math]]

and such that the degree of [math]R(x)[/math] is less than the degree of [math]D(x)[/math]. (The letters [math]N[/math], [math]D[/math], [math]Q[/math], and [math]R[/math] have been used to suggest, respectively, the words “numerator,” “denominator,” “quotient,” and “remainder.”) The first step in the method of integration by partial fractions is to write the given rational function [math]\frac{N(x)}{D(x)}[/math] in the form of equation (1). Since we can obviously integrate the polynomial [math]Q(x)[/math], we need next consider only rational functions in which the degree of the numerator is less than the degree of the denominator. If we start with such a function, then no division is necessary. Example

Write the function [math]\frac{x^4 - 4x^3 + 8x^2 - 7x +3}{x^3 - 2x^2 + 3x - 4}[/math] as the sum of a polynomial and a rational function in which the degree of the numerator is less than the degree of the denominator. Dividing, we have

[[math]] \begin{array}{rl} &\;\;x - 2 \\ x^3 - 2x^2 + 3x - 4&)\;\overline {x^4 - 4x^3 + 8x^2 - 7x + 3} \\ &x^4 - 2x^3 + 3x^2 - 4x \\ &\overline{\;\;\;\;\;\;\;\, -2x^3 + 5x^2 - 3x + 3} \\ &\underline{-2x^3 + 4x^2 - 6x + 8} \\ &x^2 + 3x - 5. \end{array} [[/math]]

It follows that

[[math]] \frac{x^4 - 4x^3 + 8x^2 - 7x + 3}{x^3 - 2x^2 + 3x - 4} = (x - 2) + \frac{x^2 + 3x - 5}{x^3 - 2x^2 + 3x - 4}, [[/math]]

which gives the required sum.

Another algebraic fact about polynomials, which we shall not prove, but shall assume and use, is that any nonconstant polynomial (i.e., of degree at least 1) with real coefficients can be written as a product of linear and quad ratic factors, each with real coefficients. By a linear factor we mean a polynomial [math]L(x)[/math] of degree 1; that is, [math]L(x) = ax + b[/math] and [math]a \neq 0[/math]. Similarly, a quadratic factor is a polynomial [math]Q(x)[/math] of degree 2; thus [math]Q(x) = cx^2 + dx + e[/math] and [math]c \neq 0[/math]. The theorem states that, for any polynomial

[[math]] f(x) = a_{n}x^n + ... + a_{1}x + a_{0} [[/math]]

with real coefficients [math]a_i[/math] and with [math]n \geq 1[/math] and [math]a_n \neq 0[/math], there exist linear factors [math]L_{1}(x), \cdots, L_{p}(x)[/math] and quadratic factors [math]Q_1(x), \cdots, Q_q(x)[/math] with real coefficients such that

[[math]] f(x)= L_1(x) \cdots L_p(x)Q_1(x) \cdots Q_q(x). [[/math]]

Note that either [math]p[/math] or [math]q[/math] may be zero. In actual practice, such a factorization of [math]f(x)[/math] may be very difficult to find, but the theorem assures us that it exists. A polynomial is said to be irreducible if it cannot be written as the product of two polynomials each of degree greater than or equal to 1. The degree of the product of two polynomials is equal to the sum of the degrees of the factors, and it therefore follows that every linear polynomial is irreducible. It was pointed out in Section 3 that a quadratic polynomial [math]cx^2 + dx + e[/math] is irreducible over the reals if and only if its discriminant [math]d^2 - 4ce[/math] is negative. For example, the polynomials [math]x^2 + 1[/math] and [math]x^2 + x + 1[/math] are irreducible, whereas [math]x^2 + 2x + 1[/math] and [math]x^2 + 2x - 1[/math] are not. If a quadratic polynomial is not irreducible, it can be factored and written as the product of two linear polynomials. Hence the factorization of an arbitrary nonconstant polynomial into linear and quadratic factors, as described in the preceding paragraph, can always be done so that all the factors are irreducible. Returning specifically to the method of integration by partial fractions, we consider a rational function [math]\frac{N(x)}{D(x)}[/math], with the degree of [math]D(x)[/math] greater than the degree of [math]N(x)[/math]. The second step is to write the denominator [math]D(x)[/math] as a product of irreducible factors. Having done so, we have

[[math]] \begin{equation} D(x) = L_1(x) \cdots L_p(x)Q_1(x) \cdots Q_q(x), \label{eq7.4.2} \end{equation} [[/math]]


where, for each [math]i = 1, . . ., p[/math],

[[math]] L_i(x) = a_{i}x + b_{i}, \;\;\; a_i \neq 0, [[/math]]

and, for each [math]j = 1, . . ., q[/math],

[[math]] Q_j(x) = c_{j}x^2 + d_{j} x + e_{j}, \;\;\; d_{j}^2 - 4c_{j}e_{j} \lt 0. [[/math]]

There is no reason to suppose that the factors which appear in equation (2) will all be distinct, and it may very well happen that [math]L_{1}(x) = L_{2}(x)[/math], etc. However, the theory is simpler if no repetitions occur, and we shall consider that case first. Case 1. The irreducible factors of $D(x)$ are all distinct. The algebraic theory of partial fractions, which we shall assume, tells us that we can write [math]\frac{N(x)}{D(x)}[/math] as a sum of rational functions each of which has one of the factors of [math]D(x)[/math] as its denominator and such that in each term the degree of the numerator is less than the degree of the denominator. Moreover, given the factorization of [math]D(x)[/math] into irreducibles, this decomposition is unique except for the order in which the terms are written. Thus

[[math]] \begin{eqnarray*} \frac{N(x)}{D(x)} &=& \frac{A_1}{a_{1}x + b_{1}} + \frac{A_2}{a_{2}x + b_{2}} + \cdots + \frac{A_p}{a_{p}x + b_{p}}\\ && + \frac{B_{1}x + C_1} {c_{1}x^2 + d_{1}x + e_{1}} + \frac{B_{2}x + C_2}{c_{2}x^2 + d_{2}x + e_2}\\ && + \cdots + \frac{B_{q}x + C_q}{c_{q}x^2 + d_{q}x + e_q}, \end{eqnarray*} [[/math]]


where each of the letters [math]A_i[/math], [math]B_j[/math] and [math]C_k[/math] represents a uniquely determined real constant. The rational functions which appear on the right side are called the partial fractions of the decomposition of [math]\frac{N(x)}{D(x)}[/math]. We shall show by means of examples how the constants in the partial fractions decomposition are determined. Consider the rational function [math]\frac{1}{(x^2 + 1)(x - 2)}[/math] discussed at the beginning of the section. The degree of the numerator, zero, is already less than 3, the degree of the denominator. Moreover, the denominator is already factored into irreducibles. Hence, we seek constants [math]A[/math], [math]B[/math], and [math]C[/math] such that

[[math]] \frac{1}{(x^2 + 1)(x - 2)} = \frac{A}{x - 2} + \frac{Bx + C}{x^2 + 1}. [[/math]]

Adding the two fractions on the right side, we have

[[math]] \frac{1}{(x^2 + 1)(x - 2)} = \frac{A(x^2 + 1) + (Bx + C)(x - 2)}{(x^2 + 1)(x - 2)}. [[/math]]

The fact that a nonzero polynomial of degree [math]n[/math] has at most [math]n[/math] distinct roots implies that two rational functions with the same denominator are equal if and only if their numerators are equal (see Problem 8 at the end of this section). Hence the equation

[[math]] \begin{equation} 1 = A(x^2 + 1) + (Bx + C)(x - 2) \label{eq7.4.3} \end{equation} [[/math]]

is true for all real values of [math]x[/math]. There are two common ways to find the constants [math]A[/math], [math]B[/math], and [math]C[/math]. One is to multiply and regroup the terms in (3) to obtain the equation

[[math]] 0 = (A + B)x^2 + (C - 2B)x + (A - 2C - 1), [[/math]]

which also holds for all real values of [math]x[/math]. The only polynomial with infinitely many roots is the zero polynomial, i.e., the polynomial with no nonzero coefficients. Hence the three coefficients on the right side of the preceding equation are all equal to zero, and we can therefore find [math]A[/math], [math]B[/math], and [math]C[/math] by solving the system of equations

[[math]] \begin{array}{ccrcrc} A & + & B & & & = 0,\\ & - & 2B & + & C & = 0,\\ A & & & - & 2C & = 1. \end{array} [[/math]]

Usually simpler is the technique in which we take advantage of the fact that (3) must be true for all values of [math]x[/math], and we choose values cleverly to help evaluate the constants. For example, letting [math]x = 2[/math] in (3), we have

[[math]] \begin{eqnarray*} 1 &=& A(2^2 + 1) + (B \cdot 2 + C)(2 - 2)\\ &=& 5A. \end{eqnarray*} [[/math]]

Hence [math]A = \frac{1}{5}[/math]. lf we then let [math]x = 0[/math], we have


[[math]] \begin{eqnarray*} 1 &=& \frac{1}{5}(0^2 + 1) + (B \cdot 0 + C)(0 - 2) \\ &=& \frac{1}{5} - 2C, \end{eqnarray*} [[/math]]

or, equivalently, [math]2C = \frac{1}{5} - 1[/math], and so [math]C = - \frac{2}{5}[/math]. Finally, choosing [math]x = 1[/math], we get

[[math]] \begin{eqnarray*} 1 &=&\frac{1}{5}(1^2 + 1) + (B \cdot 1 - \frac{2}{5})(1 - 2)\\ &=& \frac{2}{5} - B + \frac{2}{5}, \end{eqnarray*} [[/math]]

from which we conclude that [math]B = - \frac{1}{5}[/math]. Whichever method we use, we have [math]A = \frac{1}{5}[/math], [math]B = - \frac{1}{5}[/math], and [math]C = -\frac{2}{5}[/math], from which it follows that

[[math]] \frac{1}{(x^2+ 1)(x - 2)} = \frac{1}{ 5}\frac{1}{ x-2} - \frac{1}{ 5} \frac{x + 2}{x^2+ 1}, [[/math]]

the form which we integrated at the beginning of the section. Example

Integrate [math]\int \frac{dx}{a^2 - x^2}[/math]. Since [math]a^2 - x^2 = (a - x)(a + x)[/math], we decompose [math]\frac{1}{a^2 - x^2}[/math] into partial fractions [math]\frac{A}{a + x}[/math] and [math]\frac{B}{a - x}[/math]. Thus

[[math]] \frac{1}{a^2 - x^2} = \frac{A}{a + x} + \frac{B}{a - x} = \frac{A(a - x) + B(a + x)}{a^2 - x^2}. [[/math]]

Equating numerators on the left and right, we get

[[math]] 1 = A(a - x) + B(a + x). [[/math]]

Letting [math]x = a[/math], we obtain the equation [math]1 = A \cdot 0 + B \cdot 2a[/math], and so [math]B = \frac{1}{2a}[/math]. Similarly, setting [math]x = - a[/math], we get [math]1 = A \cdot 2a + B \cdot 0[/math], from which it follows that [math]A = \frac{1}{2a}[/math]. Thus


[[math]] \frac{1}{a^2 - x^2} = \frac{1}{2a} \frac{1}{a + x} + \frac{1}{2a} \frac{1}{a - x}, [[/math]]

and, therefore,

[[math]] \begin{eqnarray*} \int \frac{ dx}{a^2 - x^2} &=& \frac{1}{2a} \int \frac{dx}{a + x} + \frac{1}{2a} \int \frac{dx}{a - x}\\ &=& \frac{1}{2a} \ln |a + x| - \frac{1}{2a} \ln |a - x| + c \\ &=& \frac{1}{2a} \ln \Big| \frac{a + x}{a - x} \Big| + c. \end{eqnarray*} [[/math]]


Thus the third step in applying this method of integration is the decomposition into partial fractions, and the fourth and final step is the integration of the partial fractions. We shall show later in the section that it is always possible to carry out the last step, but, as the next example shows, doing so can be tedious. Example

Integrate [math]\int \frac{3x^3 + x^2 -14x + 46}{(x^2 + x + 1)(x^2 - 5x - 14)} dx[/math]. The degree of the numerator is 3 and that of the denominator is 4, so we proceed to the factorization of the denominator into irreducibles. It is already written as the product of two quadratics, of which [math]x^2 + x + 1[/math] is irreducible but [math]x^2 - 5x -14[/math] is not, since [math]x^2 - 5x - 14 = (x - 7)(x + 2)[/math]. Hence the form of the partial fractions decomposition is

[[math]] \frac{3x^3 + x^2 -14x + 46}{(x^2 + x + 1)(x - 7)(x + 2)} = \frac{Ax + B}{x^2 + x + 1} + \frac{C}{x - 7} + \frac{D}{x + 2}. [[/math]]

The sum of the three fractions on the right side is

[[math]] \begin{array}{l} (Ax + B)(x - 7)(x + 2)\\ \underline{+ C(x^2 + x + 1)(x + 2) + D(x^2 + x + 1)(x - 7)} \\ (x^2 + x + 1)(x - 7)(x + 2) \end{array}. [[/math]]

Equating numerators, we have

[[math]] \begin{eqnarray*} &&(Ax + B)(x-7)(x+ 2) + C(x^2 + x + 1)(x + 2) + D(x^2 + x + 1)(x-7) \\ && = 3x^3 + x^2 - 14x + 46. \end{eqnarray*} [[/math]]

If we set [math]x = 7[/math] in this equation, then

[[math]] C \cdot 513= 1026 \;\;\;\mbox{or}\;\;\; C = 2. [[/math]]

Letting [math]x = - 2[/math], we obtain

[[math]] D \cdot (- 27) = 54 \;\;\; \mbox{or} \;\;\; D = - 2. [[/math]]

If we let [math]x = 0[/math], then

[[math]] B \cdot (- 14) + C \cdot 2 + D \cdot ( - 7) = 46 [[/math]]

or

[[math]] -14B = 46 - 2C + 7D = 46 - 4 - 14 = 28 \;\;\;\mbox{or} \;\;\; B = - 2. [[/math]]

Finally, letting [math]x = - 1[/math], we have

[[math]] (-A + B) \cdot (-8) + C + D \cdot (-8) = 58 [[/math]]

or

[[math]] 8A = 58 + 8B - C + 8D = 58 - 16 - 2 - 16 = 24 \;\;\; \mbox{or} \;\;\; A = 3. [[/math]]

The partial fractions decomposition is, therefore,

[[math]] \begin{equation} \frac{3x^3 + x^2 - 14x + 46}{(x^2 + x + 1)(x^2 - 5x -14)} = \frac{3x - 2}{x^2 + x + 1} + \frac{2}{x - 7} - \frac{2}{x + 2}. \label{eq7.4.4} \end{equation} [[/math]]


Except for the first, the terms on the right side are easily integrated. The first term can be integrated by writing it as the sum of two fractions. We use the identity

[[math]] \frac{Bx + C}{cx^2+ dx+e} = \frac{B}{2c} \frac{2cx + d}{cx^2+ dx + e} + \Bigl(C - \frac{dB}{2c} \Bigr) \frac{1}{cx^2+ dx + e}. [[/math]]

Note that the numerator [math]2cx + d[/math] is the derivative of [math]cx^2 + dx + e[/math]. Both of these fractions can be integrated. In this example, we have

[[math]] \frac{3x - 2}{x^2 + x + 1} = \frac{3}{2} \frac{2x + 1}{x^2 + x + 1} - \frac{7}{2} \frac{1}{x^2 + x + 1}. [[/math]]

Then

[[math]] \int \frac{2x + 1}{x^2 + x + 1} dx = \ln (x^2 + x + 1) + c. [[/math]]

For the second fraction, we complete the square in the denominator. The result is

[[math]] \frac{1}{x^2 + x + 1} = \frac{1}{\Bigl(x + \frac{1}{2}\Bigr)^2 + \Bigl(\frac{\sqrt 3}{2}\Bigr)^2}, [[/math]]

and, since [math]\int \frac{dy}{y^2 + a^2} = \frac{1}{|a|} \arctan \frac{y}{|a|}[/math], it follows that

[[math]] \begin{eqnarray*} \int \frac{dx}{x^2 + x + 1} = \int \frac{dx}{\Bigl(x + \frac{1}{2}\Bigr)^2 + \Bigl(\frac{\sqrt 3}{2}\Bigr)^2} &=& \frac{2}{\sqrt 3} \arctan \Bigl(\frac{x + \frac{1}{2}}{\frac{\sqrt 3}{2}} \Bigr) + c\\ &=& \frac{2}{\sqrt 3} \arctan \Bigl(\frac{2x + 1}{\sqrt 3} \Bigr) + c. \end{eqnarray*} [[/math]]


Hence

[[math]] \int \frac{3x - 2}{x^2 + x + 1} dx = \frac{3}{2} \ln(x^2 + x + 1) - \frac{7}{\sqrt 3} \arctan \Bigl(\frac{2x + 1}{\sqrt 3} \Bigr) + c. [[/math]]

Returning to equation (4), we therefore get the final integral

[[math]] \begin{eqnarray*} \int \frac{3x^3 + x^2 -14x + 46}{(x^2 + x + 1)(x^2 - 5x - 14)} dx &=& \frac{3}{2} \ln(x^2 + x + 1) - \frac{7}{\sqrt 3} \arctan \Bigl(\frac{2x+ 1}{\sqrt 3} \Bigr)\\ && + 2 \ln |x - 7| - 2 \ln |x + 2| + c, \end{eqnarray*} [[/math]]


and this completes the example.

We consider next the situation in which the factorization of the denominator of [math]\frac{N(x)}{D(x)}[/math] as shown in equation (2), contains repeated factors. \medskip Case 2. The irreducible factors of $D(x)$ are not all distinct. We assume, as in Case 1, that the degree of [math]N(x)[/math] is less than the degree of [math]D(x)[/math]. There is still a unique decomposition of [math]\frac{N(x)}{D(x)}[/math] into the sum of partial fractions, but now it is more complicated. By regrouping, we may write the factorization of [math]D(x)[/math] into irreducibles as

[[math]] \begin{equation} D(x) = [L_{1}(x)]^{m_1} \cdots [L_{r}(x)]^{m_r} [Q1(x)]^{n_1} \cdots [Q_{s}(x)]^{n_s}, \label{eq7.4.5} \end{equation} [[/math]]

where [math]m_1, ... , m_r[/math] and [math]n_1, ... , n_s[/math] are positive integers, the factors [math]L_{i}(x) = a_{i}x + b_{i}[/math] are all distinct, and the factors [math]Q_{j}(x) = c_{j}x^2 + d_{j}x + e_{j}[/math] are all distinct In this case, [math]\frac{N(x)}{D(x)}[/math] is the total sum of the following individual sums of partial fractions: For each [math]i = 1, . . ., r[/math], there is the sum

[[math]] \frac{A_{i1}}{a_{i}x + b_{i}} + \frac{A_{i2}}{(a_{i}x + b_{i})^2} + \cdots + \frac{A_{im_i}}{(a_{i}x + b_{i})^{m_i}}, [[/math]]

in which the [math]Aik[/math] are uniquely determined real constants. Similarly, for each [math]j = 1, ... , s[/math], there is the sum

[[math]] \frac{B_{j1}x + C_{j1}}{c_{j}x^2+ d_{j}x + e_{j}} + \frac{B_{j2}x + C_{j2}}{(c_{j}x^2 + d_{j}x + e_{j})^2} + \cdots + \frac{B_{jn_j}x + C_{jn_j}} {(c_{j}x^2 + d_{j}x + e_{j})^{n_j}}, [[/math]]

in which the [math]B_{jk}[/math] and [math]C_{jk}[/math] are uniquely determined real constants.

Example

Integrate the rational function [math]\frac{2x^2 + x + 2}{x(x - 1)^3}[/math] by the method of partial fractions. The degree of the numerator, 2, is less than that of the denominator, 4. So we turn at once to the decomposition into partial fractions. Since the irreducible factor [math]x - 1[/math] is repeated three times, the decomposition is of the form

[[math]] \frac{2x^2 + x + 2}{x(x - 1)^3} = \frac{A}{x - 1} + \frac{B}{(x - 1)^2} + \frac{C}{(x - 1)^3} + \frac{D}{x}. [[/math]]
The sum on the right side is equal to

[[math]] \frac{Ax(x - 1)^2 + Bx(x - 1) + Cx + D(x -1)^3}{x(x - 1)^3}. [[/math]]
Equating numerators, we obtain the equation

[[math]] \begin{equation} Ax(x - 1)^2 + Bx(x - 1) + Cx + D(x - 1)^3 = 2x^2 + x + 2, \label{eq7.4.6} \end{equation} [[/math]]
which holds for all real values of [math]x[/math]. Setting [math]x = 1[/math], we obtain

[[math]] C \cdot 1 = 2 \cdot 1^2 + 1 + 2, \;\;\;\mbox{whence}\; C = 5. [[/math]]
If we let [math]x = 0[/math], then

[[math]] D \cdot ( - 1)^3 = 2, \;\;\; \mbox{whence}\; D = - 2. [[/math]]
Thus, equation (6) has become

[[math]] Ax(x - 1)^2 + Bx(x - 1) + 5x -2(x - 1)^3 = 2x^2 + x + 2. [[/math]]
In this equation we let [math]x = 2[/math], getting

[[math]] A \cdot 2 + B \cdot 2 + 10 - 2 = 8 + 2 + 2 [[/math]]
or

[[math]] \begin{equation} 2A + 2B = 4. \label{eq7.4.7} \end{equation} [[/math]]
Next, setting [math]x = - 1[/math], we have

[[math]] A \cdot (- 1)(- 2)^2 + B \cdot (- 1)(- 2) - 5 + 16 = 2 - 1 + 2 [[/math]]
or

[[math]] \begin{equation} -4A + 2B = - 8. \label{eq7.4.8} \end{equation} [[/math]]


Subtracting (8) from (7), we get [math]6A = 12[/math] and so [math]A = 2[/math]. It follows that [math]B = 0[/math], and we have therefore found the partial fractions decomposition to be

[[math]] \frac{2x^2 + x + 2}{x(x - 1)^3} = \frac{2}{x - 1} + \frac{5}{(x - 1)^3} - \frac{2}{x}. [[/math]]

Hence

[[math]] \begin{eqnarray*} \int \frac{2x^2 + x + 2}{x(x -1)^3} dx &=& 2\int \frac{dx}{x - 1} + 5 \int \frac{dx}{(x - 1)^3} - 2 \int \frac{dx}{x}\\ &=& 2 \ln | x - 1| - \frac{5}{2} \frac{1}{(x - 1)^2} - 2 \ln |x| + c, \end{eqnarray*} [[/math]]


and this completes the example.

Since any rational function can be written as the sum of a polynomial and a series of partial fractions, the general problem of integrating a rational function reduces to three integration problems: (1) integration of a polynomial; (2) integration of functions of the form [math]\frac{A}{(ax + b)^m}[/math], where [math]m[/math] is a positive integer and [math]a \neq 0[/math]; and (3) integration of functions of the form [math]\frac{Bx + C}{(cx^2 + dx + e)^n}[/math] where [math]n[/math] is a positive integer and [math]d^2 - 4ce \lt 0[/math]. The first, of course, offers no difficulties whatever. The second is also simple, since [math]\frac{A}{(ax + b)^m} = \frac{A}{a} \frac{a} {(ax + b)^m}[/math], and so

[[math]] \int \frac{A}{(ax + b)^m} dx = \frac{A}{a} \int \frac{a}{(ax + b)^m}dx\\ = \left \{ \begin{array}{ll} \frac{A}{a(1 - m)} \frac{1}{(ax + b)^{m-1}} + c &\mbox{if}\; m \neq 1,\\ \frac{A}{a} \ln | ax + b | + c &\mbox{if}\; m = 1. \end{array} \right. [[/math]]

The third problem can be solved, but we have seen in Example 3 that it is complicated even with [math]n = 1[/math]. It is attacked by writing

[[math]] \frac{Bx + C} {(cx^2 + dx + e)^n} = \frac{B}{2c} \frac{2cx + d}{(cx^2 + dx + e)^n} + \Bigl(C - \frac{dB}{2c} \Bigr) \frac{1}{ (cx^2 + dx + e)^n}. [[/math]]
The integra [math]\int \frac{2cx + d}{(cx^2 + dx + e)^n} dx[/math] is easily found, since it is of the form [math]\int \frac{du}{u^n}[/math] with [math]u = cx^2 + dx + e[/math]. The problem therefore reduces to finding [math]\int \frac{dx}{(cx^2 + dx + e)^n}[/math]. However it was demonstrated at the end of Section 3 that this integral can be evaluated by trigonometric substitution, which reduces it to an integral of the form [math]\int \cos^{2n-2} \theta d\theta[/math]. An alternative method is to use the reduction formula given in Problem 9, page 384. Thus all possible partial fractions resulting from the decomposition of a rational function can be integrated. It follows that every rational function can be integrated. The factoring of the denominator into irreducibles may be difficult, and the decomposition into partial fractions and the resulting integrations may be tedious, but the following important result has been established.

Theorem

Every rational function can be integrated by the method of partialfractions.

A table of integrals will show how to integrate many of the functions which are the partial fractions of a rational fraction. For this reason, there is no need to memorize the formulas for integration. However, it is necessary to know the technique of separating a rational function into its partial fractions in order to replace an apparently nonintegrable function by a sum of obviously integrable functions.

\end{exercise}

General references

Doyle, Peter G. (2008). "Crowell and Slesnick's Calculus with Analytic Geometry" (PDF). Retrieved Oct 29, 2024.