Using Basic Integration Formula
Using Basic Integration Formula
Formula
ˆ ˆ
. xn dx xn+1 C
2 =n + 1 + (n /= −1 13. cot(x) dx = ln | sin(x)| + C
ˆ
1
3. dx = ln |x| ˆ
+ Cx 4.
sec(x) dx = ln | sec(x) + tan(x)| +
ˆ C
4. ex dx = ex + C
15. ˆ
ˆ x
x a csc(x) dx = − ln | csc(x) + cot(x)|
5 = +
. a dx C +C
ln(a) (a > 0, a /= 1) 16.
ˆ ˆ
6. sin(x) dx = − cos(x) + C sinh(x) dx = cosh(x) + C
17.
ˆ ˆ
7. cos(x) dx = sin(x) + C 18. cosh(x) dx = sinh(x) + C (a > 0)
ˆ ˆ 1 x2
1
8. 2
sec (x) dx = tan(x) + C √ dx = sin−1
19. 2
+ Ca −x 2 (a > 0)
a
ˆ ˆ 1 x2
1 1
9. 2
csc (x) dx = − cot(x) + C 20. dx = tan−1 (a > 0)
+ C a 2 + x2 a
a
ˆ ˆ
1 1 x
10. sec(x) tan(x) dx = sec(x) + 21. √ dx = sec−1 - - (a > 0)
C
+ x x −a 2 2 a a
C ˆ 1 x2
1
√ dx = sinh−1
ˆ 22. + C a 2 + x2 a (x > a > 0
11. csc(x) cot(x) dx = − csc(x) ˆ 1 x2
1
+C √ dx = cosh−1
+ C x2 − a2 a
Example 1 - Substitution: Evaluate the integral
ˆ 5
2x − 3
√ 3 x 2−3 x +1 dx .
ˆ 5
2x − 3 ˆ 11
1
u = x2 − 3x + 1 √ dx = √ du
3 x2 − 3x + 1 1 u
du = 2x − 3 dx
=
u−1/2 du
11 ˆ
11
2
u = (3) − 3(3) + 1
=1 = 2u1/2
u = (5)2 − 3(5) + 1
= 11 -1
√
√
=22 1√
11 − 22 1
11 − 1
=
1
ˆ 8x − x2 dx.
√
8x − x2 = −(x2 − 8x) ˆ ˆ
1 1
√ dx √
= −((x − 4)2 − 42 8x − x 2
4 − (x − 4)2
2
) = ˆ dx
1
= 42 − (x − 4)2 =
2 du
2
u=x−4
=
du = dx
=
Example 3 - Trig Identities: Calculate
ˆ
cos(x) sin(2x) + sin(x) cos(2x) dx.
ˆ ˆ
cos(x) sin(2x) + sin(x) cos(2x) dx = sin (x + 2x) dx
ˆ
u = 3x = sin (3x) dx
du = 3 dx
ˆ
1 1
= sin (u) du
3
du = dx 1
3 = − cos (u) + C
3
1
= − cos (3x) + C
3
Example 4 - Trig Identities:
Find
ˆ π
4 1
0 dx.
1 − sin(x)
ˆ π ˆ π
1 1 1 + sin(x)
0 4
dx = 4
· dx
1 − sin(x) 0 1 − sin(x) 1 + sin(x)
ˆ π
4
1 + sin(x)
= dx
0 1 − sin2(x)
ˆ π
1 1 sin(x)
= 4
+ dx
0 cos2(x) cos(x) cos(x)
ˆ π
- π4
= x) + x-
tan( sec( ) -0
1π 1
= tan 2
2 + sec − (tan(0) + sec(0))
π 4
4√
=1+ 2 − (0 + 1)
=
√2
√
u=1+ x ˆ 2 u3 2
= − du
1 u2 u3
ˆ
du = √ dx = 2u−2 − 2u−3 du
2 x
√
2 xdu = dx
= −2u−1 + u−2 + C
2(u − 1) du = dx 2 1
= Properties of Trig Integrals
Exampl
e6-
− + +C
ˆ π
2
u u2
x3 cos(x) dx.
= 2 1
−π
− √ + √ 2 +C
1 + x (1 + x)
2
f (x) = x3 =⇒ f (−x) = (−x)3 = −x3 = −f (x) g(x) = cos(x) =⇒ f (−x) = cos(−x) = cos(x) = f (x)
2 2
Putting
(by these8,two
Theorem facts together we see that x3 cos(x) is an odd function and is symmetric over the
Section
5.6) # $
interval − π , π . Thus
ˆ π
2
0
x3 cos(x) dx =
−π
2
Section 8.2: Techniques of Integration
It is useful when one of the functions (f (x) or g(x)) can be differentiated repeatedly and the other
function can be integrated repeatedly without difficulty. The following are two such integrals:
ˆ ˆ
x cos(x) dx and x2ex dx.
Notice f (x) = x or f (x) = x2 can be differentiated repeatedly (they are even eventually zero) and g(x) =
cos(x) and
g(x) = ex can be integrated repeatedly without difficulty.
An Application of the Product Rule: If f (x) and g(x) are differentiable functions of x, the product rule says that
d
[f (x)g(x)] = f'(x)g(x) + f (x)g'(x). dx
Integrating both sides and rearranging gives us the Integration by Parts formula!
ˆ d ˆ
[f (x)g(x)] dx = f'(x)g(x) dx + f (x)g'(x) dx
ˆ
dx
ˆ ˆ ˆ
d
=⇒ '
f (x)g (x) dx = [f (x)g(x)] dx − f'(x)g(x) dx
ˆ dx ˆ
f (x)g'(x) dx = f (x)g(x) − f'(x)g(x) dx
=⇒
ˆ
x cos(x) dx.
We have to decide what to assign to u and what to assign to dv. Our goal is to make the integral easier.
One thing to bear in mind is that whichever term we let equal u we need to differentiate - so if
differentiating makes a part of the integrand simpler that’s probably what we want! In this cases
differentiating cos(x) gives − sin(x), which is no easier to deal with. But differentiating x gives 1 which is
simpler. So we have,
ˆ ˆ
u=x dv = cos(x) dx
x cos(x) dx = x sin(x) − sin(x) dx
du = dx v = sin(x)
=
x sin(x) + cos(x) + C
Example 2: Evaluate
ˆ
x2ex dx.
Here we go through the same thought process. If u = ex then du = ex dx, which doesn’t make the problem
1 3
any easier (though it doesn’t make it any harder either). But in this
3 case dv = x2 would give v = x
which arguably is not simpler that x2. So,
u = x2 dv = ex dx ˆ ˆ
du = 2x dx v=e x 2 x 2 x
x e dx = x e − 2
xex dx.
It’s at this point we see that we still cannot integrate the integral on the write easily. This is okay. Sometimes we may have to
apply the integration by parts formula more than once!
ˆ ˆ
2 x 2 x
x e dx = x e − 2 xex dx
5 ˆ 6
u=x x
dv = e dx 2 x
= x e − 2 xe − x x
e dx
= ! "
x2 − 2x + 2 ex + C
ˆ
The previous technique works for any integral of the form xnemx dx, where n is any positive integer and m
is any integer. What if n was negative? Then this case we would set u = ex.
Example 3 - Integration by Parts for Definite Integrals: Find the area of the region bounded by the curve y = xe−x
ˆ 4
A= xe−x dx
0
4
u=x dv = e−x dx
du = dx v = −e−x
ˆ 4
-4 ˆ 4
xe−x dx = xe−x −e −x dx
− - −
0
0 -0
4 4
ˆ
= −xe−x-
-0 + e−x dx
!
0
−4 4
" -
−x
= −4e −0 e
− -
0
= −4e−4 −
! "
e−4 −
= −4e−4 − !e1−4 − 1
"
= −5e−4 + 1
Example 4 - Tabular Method: In Example 2 we have to apply the Integration by Parts Formula
multiple times. There is a convenient way to “book-keep” our work. This is done by creating a table.
Let’s see how by examining Example 2 again.
Evaluate
ˆ
x2ex dx.
x2 + ex
2x ex
−
2
ex
Then the integral is,
ˆ
x2ex dx = +x2 · ex − 2x · ex + 2 · ex + C
! "
x2 − 2x + 2 ex + C
=
We have actually used the integration by parts formula, but we have just made our lives easier by
condensing the work into a neat table. This method is extremely useful when Integration by Parts
needs to be used over and over again.
Example 5: Find the integral
ˆ
1 x3 cos (nx) dx,
π
π 0
x3 cos(nx)
+
3x2
6x − 1
sin(nx)
n
6 1
+ − n2
cos(nx)
1
− − n3
sin(nx)
1
0 cos(nx)
n4
-
53 cos(nx) − 3 sin(nx)
43 − 4 cos(nx) 46
1 3π2 n2 −
nπ n−6 −n 60
= 0 π+ cos( ) 0+ 0+0 0
5 2 6
1 n2 n 6(−1)n
3π (−1) n64 −
= − + n4
π n2 n4 n4
3 π2n2(−1)n − 2(−1)n + 2
π
= n4
cos(x):
π 2π 3π 4π
−1
1
sin(x):
π 2π 3π 4π
−1
Example 6 - Recurring Integrals: Find the integral
ˆ
ex sin(x) dx.
(2) ˆ
Notice that now the integral we are interested in, ex sin(x) dx, appears on
u = ex dv = cos(x) both the left and right hand side of the equation. So, if we add this
integral to both sides we get
du = ex dx v = sin(x) ˆ
ˆ
=⇒ ex sin(x) dx = ex e(− cos(x) + sin(x))
=⇒ ex sin(x) dx =
x
(sin(x) 2− cos(x))
This “trick” comes up often when we are dealing with the product of two functions with “non-
terminating” derivatives. By this we mean that you can keep differentiating functions like ex and trig
functions indefinitely and never reach 0. Polynomials on the other hand will eventually “terminate” and
their nth derivative (where n is the degree of the polynomial) is identically 0.
Section 8.3: Trigonometric Integrals -
Worksheet
Goal: By using trig identities combined with u-substitution, we’d like to find antiderivatives of the form
ˆ
sinm(x) cosn(x) dx
(for integer values of m and n). The goal of this worksheet1 is for you to work together in groups of 2-3 to
discover the techniques that work for these anti-derivatives.
ˆ
cos4(x) sin(x) dx.
u = cos(x) ˆ ˆ
4
cos (x) sin(x) dx = − u4 du
du = − sin(x) dx
5
u
= − + C
5 5
cos (x)
= − +C
5
Example 2: Find
ˆ
sin3(x) dx.
u = cos(x) ˆ ˆ
3
! "
du = − sin(x) dx sin (x) dx = 1ˆ − cos2(x) sin(x) dx
! "
= − 1 − u2 du
u
= −u + 3 C
+
3
cos3(x)
= − cos(x) + +C
3
1
Worksheet adapted from BOALA, math.colorado.edu/activecalc
Example 3: Find
ˆ
sin5(x) cos2(x) dx.
! "2
(Hint: Write sin5(x) as sin2(x)
sin(x).)
ˆ ˆ
5 2
! "2
sin (x) cos (x) dx = sin2(x) cos2(x) sin(x) dx
ˆ
! "2
u = cos(x) = 1 − cos2(x) cos2(x) sin(x) dx
du = − sin(x) dx ˆ
! "2
=− 1 − u2 u2 du
ˆ
= − !1 − 2u2 + u4 du
ˆ "
=− u2 − 2u4 + u6 du
3
u 2u5 u7
=−
3 − +C
+ cos 3 5
(x) 7 5(x) cos7(x)
2 cos
− + − +C
= 3 7
5
Example 4: Find
ˆ
sin7(x) cos5(x) dx.
(The algebra here is long. Only set up the substitution - you do not need to fully evaluate.)
ˆ ˆ
7 5
! "3
sin (x) cos (x) dx = sin2(x) cos5(x) sin(x) dx
ˆ
! "3
= 1 − cos2(x) cos5(x) sin(x) dx
u = cos(x)
du = − sin(x) dx ˆ
! "3
= − 1 − u2 u5 du
u = cos(x)
du = − sin(x) dx ˆ
! "(m−1)/2
= − 1 − u2 un du
Example 6: Note that the same kind of trick works when the power on cos(x) is odd. To check that
you understand, what trig identity and what u-substitution would you use to integrate
ˆ
cos3(x) sin2(x) dx?
ˆ ˆ
sin2(x) + cos2(x) 3 2
= 1 cos (x) sin (x) dx = cos2(x) sin2(x) cos(x) dx
ˆ
cos2(x) = 1 − sin2(x) ! "
= 1− sin2(x ) sin2x( ) cos(
x dx)
ˆ
u = sin(x) ! "
= 1 − u2 u2 du
du = cos(x) dx
Example 7: Now what if the power on cos(x) and sin(x) are both even? Find
ˆ
sin2(x) dx,
~~~~~~~~~~
ˆ
u = sin(x) dv = sin(x) dx = − sin(x) cos(x) − − cos(x) cos(x) dx
ˆ
du = cos(x) dx v = − cos(x) = − sin(x) cos(x) + cos2(x) dx
ˆ
= − sin(x) cos(x) + 1 − sin2(x) dx
ˆ
= − sin(x) cos(x) + x − sin2(x) dx
~~~~~~~~~~
ˆ
=⇒ 2 sin2(x) dx = − sin(x) cos(x) + x + C
ˆ
=⇒ sin2(x) dx = − sin(x) cos(x)
x +C
2
(c) Show that your answers to parts (a) and (b) above are the same by giving a suitable trig identity.
1 1
sin(x) cos(x) = 2 sin(x) cos(x) = sin(2x).
2 2
(d) How would you evaluate the integral
ˆ
sin2(x) cos2(x) dx?
ˆ ˆ
2 2 1 1
sin (x)cos (x) dx = (1 − cos(2x)) · (1 + cos(2x)) dx
~~~~~~ ~~~~~~~~~~~~~
2 2
ˆ
1
= 1ˆ− cos2(x) dx
14 1
= x− cos2(2x) dx
4 ˆ
4
11 1
= x− (1 + cos(4x)) dx
4 4 2ˆ
1 1 1
= x− x− cos(4x) dx
4 8 8
= 1 x − 1 sin(4x) + C
8 32
Example 8: Evaluate the integral in problem (2) above, again, but this time by parts using u = sin2(x) and
dv − sin(x) dx. (After this, you’ll probably need to do a substitution.)
ˆ ˆ
3
sin (x) dx = sin2(x) sin(x) dx
ˆ
2
u = sin2(x) dv = sin(x) dx = − sin (x) cos(x) − − cos(x) · 2 sin(x) cos(x) dx
ˆ
du = 2 sin(x) cos(x) dx v = − cos(x) = − sin2(x) cos(x) + 2 cos2(x) sin(x) dx
ˆ
2
= − sin (x) cos(x) − 2 w2 dw
w = cos(x) 2u3
= − sin2(x ) cos(
x −) +C
dw = − sin(x) dx 3
2 2 cos3(x)
= − sin ( x) cos(x) − +C
3
Example 9 - For fun: Can you show your answers to problem (2) and (8) above are the same? It’s
another great trigonometric identity.
2 cos3(x) 2 2 cos3(x)
2 ! "
sin (x) cos(x) − = − 1 − cos2(x) cos(x) − cos3(x) = − cos(x) + cos3(x) − cos3(x) = − cos(x) +
− 3 3 3
3
Example 10 - Further investigations: (especially for mathematics, physics and engineering majors) We also would like to be able
to solve integrals of the form
ˆ
tanm(x) secn(x) dx.
These two functions play well with each other, since the derivative of tan(x) is sec2(x), the derivative of
sec(x) is sec(x)|tan(x) and since there is a Pythagorean identity relating them. It sometimes works to use
u = tan(x) and it sometimes works to use u = sec(x). Based on the values of m and n, which substitution
should you use? Are there cases for which neither substitution works? (See page 472 of the text.)
Section 8.4: Trigonometric
Substitution
Motivation: If we want to find the area of a circle or ellipse, we have an integral of the form
ˆ √
a2 − x2 dx
u = a2 − x2
√ ñ ñ ñ ! " √
a2 − x2 = a2 − (a sin (θ))2 = 12 − a2 sin2 (θ) = a2 1 − sin2 (θ) = a2 cos2 (θ) = a |cos (θ)| .
Common Trig Substitutions: The following is a summary of when to use each trig substitution.
If you are worried about remembering the identities, then don’t! They can all be derived easily,
assuming you know three basic ones (which by now you should):
sin2 1 sin (θ)
(θ) + cos2 (θ) = 1, sec (θ) = , tan (θ) =
√
9 x2
ˆ −
dx.
x2
√ 2 2
ˆ 2 ˆ √ 2
x = 3 sin (θ) 9− 3 − 3 sin
dx = · 3 cos (θ) dθ
x x 2
(θ) 32 sin2 (θ)
è π π ˆ √
θ é = 3/ 1 − sin2 (θ)
, · /3 cos (θ)
∈ −2 2 2
2 /3/ sin (θ)
ˆ
√
=
dx = 3 cos (θ) dθ
cos2 (θ)
·
sin2 (θ) cos (θ) dθ
cos2 (θ)
ˆ dθ
= sin2 (θ)
ˆ
= cot2 (θ) dθ
ˆ
= csc2 (θ) − 1 dθ
= − cot (θ) − θ + C
√
3 2 − x2
= − − arcsin (θ) + C
How did we recover x? x
x A2 + x 2 = 3 2
x = 3 sin (θ) =⇒ = sin
(θ) 3
hyp. 3 A2 = 3 2 − x 2
opp. =⇒ x
√
θ A= 3 2 − x2
θ
A
cot ( adj.
1 √
adj
. 3 2 − x2
θ) = = =
tan (θ) opp. x
This is a common process in trig substitution. When you substitute back for your original variable, in this
case x, you will always be able to find the correct substitutions by drawing out and labelling a right triangle
correctly.
Example 2: Find
ˆ 1
√
x2 + 4 dx.
x2
ˆ 1 2 sec2 (θ)
x = 2 tan (θ) √ dx √ dθ
x2 x2 + 4 ˆ
=ˆ 22 tan2 (θ) 22 tan2 (θ) + 22
2/ sec2 (θ)
1 π π2 =
√ dθ
θ∈ − , 22 tan2 (θ) 2/ tan 2 (θ) + 1
2 2
dθ ˆ
1 cos (θ)
=4 sin2 (θ) dθ
u = sin (θ) ˆ
1 1
= du
4 u2
du = cos (θ) 1 1
= − + C
4u
1
= − +
4 sin
C
(θ)
1
= − csc (θ) + C
4√
x2 + 4
= − +C
How did we recover 4x
x?
x = 2 tan (θ) =⇒
x = hyp. H
2 opp. =⇒ x
tan (θ) H2 = x2 + 22
θ
√
adj. 2 H = x2 + 4
hyp √
. x2 + 4
csc
( 1
θ) = = =
sin (θ) opp. x
Example 3: Evaluate
ˆ x2
√ dx.
9 − x2
ˆ ˆ 3 cos ( )
x2
32 sin2 (θ)
√ dx = √ · θ dθ
9 − x2 32 − 32 sin2 (θ)
32 sin2 (θ)
√ ·
ˆ 3/ 1 − sin2 (θ) 3/ cos (θ) dθ
=
x = 3 sin (θ) 32 sin2 (θ)
ˆ cos (θ) dθ
= √ ·
1 π π2 2
θ ˆ cos (θ)
, =
∈ −2 9 sin2 (θ)
2
dx = 3 cos (θ) dθ ˆ
9
= 1 − cos (2θ) dθ
2 3 4
9 1
= θ − sin(2θ) + C
2 2
9
= (θ − sin(θ) cos(θ)) + C
2 A √ B
1 x2 C
9 −1 x 9 − x2
= sin − +
2 3 3
3
x
x = 3 sin (θ) =⇒ = hyp. 3
3 opp. =⇒ x
sin (θ)
3 2 = x 2 + A2
θ θ √
A = 9 − x2
adj. A
cos ( √
9 − x2
adj
.
θ) = =
hyp. x
Section 8.5: Integration by Partial
Fractions
Our next technique: We can integrate some rational functions using u-substitution or trigonometric
substitution, but these methods do not always work. Our next method of integration allows us to
express any rational function as a sum of functions that can be integrated using methods with which we
are already familiar. That is, we cannot integrate
1
x2 − x
but it is equivalent to
as-is, 1 1
− ,
x x−1
each term of which we can
integrate.
ˆ
x−7
dx.
(
x + 1) (x − 3)
ˆ
1
(a) dx = ln |x + 1| + C
x+1
(b) 2 1 2 (x − 3) − (x + 1) 2x − 6 − x − 1 x−7
− = = =
x+1 x−3 (x + 1) (x − 3) (x + 1) (x − 3) (x + 1) (x − 3)
ˆ ˆ
x−7 2 1
(c) dx = − dx = 2 ln |x + 1| − ln |x − 3| + C
(
x + 1) (x − 3) x+1 x−3
ˆ
2: Compute
Example ( 10x − 31 dx.
x − 1) (x − 4)
(a) 7 3 7 (x − 4) + 3 (x − 1) 10x − 31
+ = =
x−1 x−4 (x − 1) (x − 4) (x − 1) (x − 4)
ˆ ˆ
10x − 31 7 3
(b) dx = + dx = 7 ln |x − 1| + 3 ln |x − 4| + C
(
x − 1) (x − 4) x−1 x−4
The previous two examples were nice since we were given a different expression of our integrand
before hand. But what about when we don’t? It is clear that the key step is decomposing our integrand
into simple pieces, so how do we do it? The next example outlines the method.
ˆ
x + 14
Example 3: Goal: Compute
( + 5) (x + 2)
x dx.
first step is to
decompose x + 14
Our ( as
x + 5) (x + 2)
x + 14
? ?
= + .
( x + 5) (x + 2)
x+5 x+2
There is no indicator of what the numerators should be, so there is work to be done to find them. If we
let the numerators be variables, we can use algebra to solve. That is, we want to find constants A and B
that make the equation below true for all x /= −5, −2.
x + 14
A B
= + .
( x + 5) (x + 2)
x+5 x+2
We solve for A and B by cross multiplying and equating the numerators.
x + 14
( A B A (x + 2) + B (x + 5)
= + = =⇒ x + 14 = A (x + 2) + B (x + 5)
x + 5) (x + 2) x+5 x+2 (x + 5) (x + 2)
= Ax + 2A + Bx + 5B
= (A + B) x + 2A + 5B
1 = A + B =⇒ B = 1 − A
14 = 2A + 5B
ˆ ˆ
=2 + 5 (1 ) x + 14 −3 4
dx = + dx
A −A (x + 5) (x + 2) x+5 x+2
= 2A + 5 − 5A
= −3 ln |x + 5| + 4 ln |x + 2| + C
= 5 − 3A
=⇒ 9 = −3A
=⇒ −3 = A
=⇒ B = 1 − (−3)
=4
Example 4: Find ˆ
dx.
x + 15
(3
x − 4) (x + 1)
x + 15
A B A (x + 1) + B (3x − 4)
(3 = + = =⇒ x + 15 = A (x + 1) + B (3x − 4)
x − 4) (x + 1) 3x − 4 x+1 (3x − 4) (x + 1)
= Ax + A + 3Bx − 4B
= (A + 3B) x + A − 4B
1 = A + 3B =⇒ A = 1 − 3B
15 = A − 4B
= (1 − 3B) − 4B 3 − dx
ˆ x−4 x+1
= 1 − 7B ˆ x + 15 7 2
7
=⇒ 14 = −7B dx = ln | 3x + 5| − 2 ln |x + 1| + C
(3x − 4) (x + 1)
=⇒ −2 = B 3
=
=⇒ A = 1 − 3 (−2) = 7
ˆ
Example 4 - An alternative approach: Find
(3
x + 15
dx.
x − 4) (x + 1)
x + 15
A B A (x + 1) + B (3x − 4)
(3 = + = =⇒ x + 15 = A (x + 1) + B (3x − 4)
x − 4) (x + 1) 3x − 4 x+1 (3x − 4) (x + 1)
Instead of expanding everything, comparing coefficients and solving a system of linear equations,
sometimes it may be helpful to plug in strategic values of x to solve. Good values to choose are those
that are roots of the polynomials that appear on the denominators of the fraction. Observe,
4
! 4
" !! 4
" " ! ! 4
"
x = −1 : (−1) + 15 = A ((−1) + 1) + B (3(−1) − 4) x= 3
: 3
+ 15 = A
3
+ 1 +3 B 3
"
−4
49
=⇒ 14 = 0 − 7B =⇒ 3 = 73A + 0
=⇒ ˆ −2 = B =⇒ 7=A
ˆ
7 2
(3 x + 15 dx = − dx =
7
ln | 3x + 5| − 2 ln |x + 1| + C
x − 4) (x + 1) 3x − x+1 3
4
ˆ
5x − 2
Example 5: Goal: Find dx.
(x + 3)2
Here, there are not two different linear factors in the denominator. This CANNOT be expressed in the form
5x − 2 5x − 2 A B A+B
= /= + = .
(x + 3)2 (x + 3)(x + 3) x+3 x+3 x+3
5x − 2 A B A(x + 3) + B
= + = =⇒ 5x − 2 = A(x + 3) + B
ˆ ˆ
5x − 2 5 17 17
dx = − dx = ln |x + 3| + +C
(x + 3)2 5 x+3
x+3 (x + 3)2
Example 6: What if the denominator is an irreducible quadratic of the form x2 + px + q? That is, it can
! n
not be factored (does not have any real roots). "In this case, suppose that x2 + px + q is the highest
power of this factor that divides the denominator. Then, to this factor, assign the sum of the n partial
fractions:
B 1x + C 1 B 1x + C 1 B 1x + C 1 B 1x + C 1
+
2 + 3 + ···+ .
(x2 + px + q) (x2 + px + q) (x2 + px + q)
2
(x + px + q) n
Compute
ˆ 2x + 4
2
− dx.
( 2
x + 1)(x − 1)
There are four unknowns here, A, B, C and D. In this case we’re going to want to minimise the amount of work we do here. In general it
is going to be beneficial to solve for as many coefficients as we can by plugging in numbers, and then expand everything to compare coefficient
after reducing the workload.
=⇒ 2 = 0 + 0 + 2D
=⇒ 1=D
So we got one coefficient this way. That’s better than nothing! Now if we use this new information and then rearrange a little we end up
with less solving to do. This does require you however to be comfortable with algebra.
Now we have already seen what happens when x = 1, so we can go right ahead and divide by the (x−1) term that appears on both sides.
=⇒ −x − 3 = (Ax + B) (x − 1) + C( x2 + 1)
ficients. When there are lots of coefficients it is a good idea of coming
= Ax2 + Bx − Ax − B + 2
+C up with a way to book-keep your algebra - it can get very messy if you
Cx
don’t. Below is just one way you can do it.
2
= (A + C) + (B − A)x + C − B
x
(1) 0 = A +C (1) 0 = A +C (1) −4 = 2C
=
(2)+(3)
=
(1)+(2)
− ⇒ − ⇒ −
(2) 1 = −A +B (2) 4 = −A +C (2) 4 = −A +C
(3) −3 = −B +C (3) −3 = −B +C (3) −3 = −B +C
ˆ ˆ
x2 +−2x
1)(x+
−41)2 2x + 1 2 1
( dx = x2 + 1 − − 1)2
x − 1 (x+ dx
=⇒ −4 = 2C =⇒ −2 = C
=⇒ −4 = −A − =⇒ 2=A ˆ
x22x+ 1 x21 + 1 x − 2 1 (x − 1 1)2
2 = + − + dx
=⇒ −3 = −B − =⇒ 1 = B, So . . . = 1
ln( 2 −1
+C
2
x + 1) + tan (x) − 2 ln |x − 1| −
x− 1
Summary: Method of Partial Fractions when
f (x)
(
g x) is proper (deg f (x) < deg g(x))
1. Let x − r be a linear factor of g(x). Suppose that (x − r)m is the highest power of x − r that divides
g(x). Then, to this factor, assign the sum of the m partial fractions:
A1 A2 A3 + · · · + Am
+ .
( 2
+ 3 m
x − r) (x − r) (x − r) (x − r)
B 1x + C 1 B 1x + C 1 B1x + C1
+ B1 x + C 1 + · · · + .
2
+
(x2 + px + q) (x2 + px + q) (x2 + px + q) (x2 + px + q)n
3
3. Continue with this process with all irreducible factors, and all powers. The key things to remember are
(i) One fraction for each power of the irreducible factor that appears
(ii) The degree of the numerator should be one less than the degree of the denominator
f (x)
4. Set the original fraction equal to the sum of all these partial fractions. Clear the resulting
equation of fractions ( )
gx
5. Solved for the undetermined coefficients by either strategically plugging in values or comparing
coefficients of powers of x.
Section 8.7: Numerical Integration
What to do when there’s no nice antiderivative? The antiderivatives of some functions, like sin(x2), 1/ ln(x) and
√
1 + x4 have no elementary formulas/ When we cannot find a workable antiderivative for a function f (x) that
we have to
integrate, we can partition the interval of integration, replace f (x) by a closely fitting polynomial on
each subinterval, integrate the poynomials and add the results to approximate the definite integral of f
(x). This is an example of numerical integration. There are many methods of numerical integration but
we will study only two: the Trapezium Rule and Simpson’s Rule.
Trapezoidal Approximations: As the name implies, the Trapezium Rule for the value of a definite
integral is based on approximating the region between a curve and the x-axis with trapeziums instead
of rectangles - which, if you recall, we studied when we looked at Riemann integration in Calculus I.
f (x) f (x)
y0 y1 y2 y3 y4 ...
∆x ∆x ∆x ∆x
b−a
Assume the length of each subinterval is ∆x = . Then the area of the trapezium that lies above the x-axis
in the i th
n
δx
subinterval is Ti = (yi−1 + yi) where yi−1 = f (xi−1) and yi = f (xi). Then the area of the under the curve
and above 2
the x-axis is approximated by the sum of the trapeziums:
∆x ∆x ∆x
T = (y0 + y1) + (y1 + y2) + · · · + (
2 2 y +y )
n−1 n
2
∆x
= (y0 + y1 + y1 + y2 + · · · + yn−2 + yn−1 + yn−1 + yn)
2
∆x
= (y0 + 2y1 + 2y2 + · · · + 2yn−1 + yn)
2
A n−1
∆x B
= 0 + Σ
y yn + 2 yi
2
A B
∆x n−1
Σ
= ( x0 ( i)
f ) +f x
( n) + 2 f x
2
i=1
ˆ b
The Trapezium Rule: To approximate ( ) , use
f x dx
a
∆x
T= (y0 + 2y1 + 2y2 + · · · + 2y 1 +y )
n−
2 A n B
n−1
∆x ( )
= ( x0 Σ f xi ,
f f x ) + 2
2 )+ ( n
i=1
b−a
and ∆x = .
ˆ 2
Example 1: Use the Trapezium Rule with n = 4 to estimate x2
1 dx. Compare the estimate with the exact value.
2−
1
∆ x0 = a x1 = a + ∆x x2 = a + 2∆x x3 = a + 3∆x x4 = a + 4∆x
1
= (16 + 50 + 72 + 98 1
12 = (8 − 1)
+ 64) 3
8 7
1 =
= (300) 3
12
8
= 75
32
75 7
225 224 1
− = − = .
32 3 96 96 96
1
So the approximation overestimated the actual area by , which is pretty good considering we only
9
6
used 4 trapeziums. Just like when we looked at Riemann sums, using more trapeziums results in a
better approximation.
Parabolic Approximations: Instead of using the straight-line segments that produced the trapeziums,
we can use parabolas to approximate the definite integral of a continuous function. We partition the interval
[a, b] into n subintervals
b−a
of equal length ∆x = but this time we require n to be an even number. On each consecutive pair
of intervals
n
we approximate the curv e y = f (x) ≥ 0 by a parabola. A typical parabola passed through three
consecutive points:
(xi−1, yi−1), (xi, yi) and (xi+1, yi+1) on the curve.
n=4:
f (x)
a := x0 x1 x2 x3 x4 =: b
n=8:
f (x)
a := x0 x1 x2 x3 x4 x5 x6 x7 x8 =: b
So how do we compute the area under each parabola y = Ax2 + Bx + C? By translating we can
assume that the centre point of our parabola is at xi = 0
The area under the parabola and above the x-axis is given by
(0, yi)
ˆ
∆x
Si = Ax2 (−∆x, yi−1)
−∆x
+ Bx + C dx
∆x
Ax3 Bx2 -
= + Cx (∆x, yi+1)
3 -
2
+ -
A(∆x)3 B(∆x)2 −∆x
5 6
= + + (−∆x)3 B(−∆x)2
(∆ ) A + + ( ∆ )
3 2 C x −
∆ 3
2A∆x3
= +2
C x
3
C − x
2
∆x
= (2A∆x + 6C)
3 −∆x 0 ∆x
yi−1 = A∆x2 − B∆x + C
! " ! "
yi =C
=⇒ yi−1 + 4yi + yi+1 = A∆x2 − B∆x + C + 4C + A∆x2 + B∆x + C = 2A∆x + 6C
yi+1 = A∆x2 + B∆x + C
∆x
=⇒ S i yi−1 + 4y i + y i+1 )
= 3 (
So if we sum up the areas under all of the parabolas, we obtain our approximation.
ˆ b
Simpson’s Rule: To approximate ( ) , use
f x dx
a
∆x
S= (y0 + 4y1 + 2y2 + 4y3 + · · · + 2yn−2 + 4yn−1 + yn)
3
n−1
∆x Σ 2
= f (x0) + f (x ) + 2 f (x2 1) + 2f (x2 ) ,
n i− i
3 i=1
b−a
and ∆x = with n an even number.
n
ˆ 2
4
Example 2: Use the Simpson’s Rule with = 4 to approximate 5 . Compare
the estimate with the exact value.
n x dx
0
2−
0
∆ x1 = a + ∆x x2 = a + 2∆x x3 = a + 3∆x x4 = a + 4∆x
x0 = a
33 2 5x4 dx = x5
1 0 2
1
16 2 81 2 4 2
128
2
= 5 + 4 ·5 + 2 ·5 + 4 ·5 0 -0
6 16 16 16 16 -
+5
16 = 25 − 05
5
= (0 + 4 + 32 + 324
+ 256)
96 = 32 − 0
5
= (616) = 32
9
385
= 12
1
So the approximation overestimated the actual area by , which is pretty good considering we only used 2
1
parabolas.
2
Just like Riemann sums and the Trapezium rule, using more parabolas results in a better
approximation. In fact, of the three rules Simpson’s Rule gives the best approximation. This can be
seen by looking at the error estimates.
Error Estimates in the Trapezium and Simpson’s Rules If f''(x) is continuous and M is any upper bound for the values of |f''(x)|
on [a, b], then the error ET in the Trapezium Rule for approximating the definite integral of f (x) over the interval [a, b] using n
trapeziums satisfies the inequality
3
|ET | ≤ M (b − a)
.
12n2
If f (4)(x) is continuous and M is any upper bound for the values of |f (4)(x)| on [a, b], then the error ES in
n
Simpson’s Rule for approximating the definite integral of f (x) over the interval [a, b] using parabolas
satisfies the inequality 2
5
|ES| ≤ M (b − a)
.
180n4
First we differentiate f (x) 4 times and check that it is continuous on the interval [0, 2].
f (x) = 5x4
f'(x) = 20x3
f''(x) = 60x2
f'''(x) = 120x
f (4)(x) = 120
Thus M = 120 works as a bound. So, with n = 4, the error is bounded by:
M (b − a)5 120(2 − 0)5 120 · 25 1 1
|ES | ≤ = = = = .
180n4 180(4)4 180 · 28 3 · 22 12
To achieve an approximation with |ES| ≤ 0.001, we again find a bound for M but this time we solve the
inequality for n.
26 1 4
=
⇒ ≤ n
3 1000
26 · 1000 4
=
⇒ ≤n
28 · 23·53 4
=
⇒ ò ≤n
4
32 53
=⇒ 4 · ≤n
3
ñ
2 3
·54
So setting n ≥ 4 3 ≈ 12.086 would ensure an approximation of the desired accuracy.
Section 8.8: Improper Integrals
Switching up the Limits of Integration: Up until now, we have required two properties of definite integral:
We will now see what happens if we allow the domain or range to be infinite!
Infinite Limits of Integration: Let’s consider the infinite region (unbounded on the right) that lies under the curve
f (x) = e−x/2
First, we examine what the area looks like over finite intervals. That is, we integrate over [0, b].
b
b è é=2 2
A(b) := ˆ = −2e−b/2 − −0/2
1 1 − e−b/2 .
e−x/2 dx = − 2e−x/2 -
0 -0 −2e
Now we have an expression for the area over a finite integral, we can let b −→ ∞ by calculating the limit of this
expression.
1 2
A = lim A(b) = lim 2 1 − e −b/2
= 2 (1 − 0) = 2.
b→∞ b→∞
So
,
ˆ ∞ b
e−x/2 dx = lim ˆ
e−x/2 dx = 2.
0 b→∞
0
So this is how we deal with infinite limits of integration - with a limit! Remember those?
Definition: Integrals with infinite limits of integration are called improper integrals of Type I.
b
ˆ ˆ f (x) dx.
f (x) dx = lim
2. If f (x) is continuous on (−∞, b],
then a b→∞ a
ˆ b ˆ b f (x) dx.
f (x) dx = lim
−∞ a→∞ −a
In each case, if the limit is finite we sat that the improper integral converges and that the
limit is the
value of the improper integral. If the limit fails to exist, the improper integral
diverges Any of the integrals in the above definition can be interpreted as an area if f
Example 1: Evaluate
ˆ ∞
ln(x)
1 dx.
x2
-b ˆ b
ˆ ln(x) 1
b
ln(x) dx = − −
- − 2
x2 x dx x
1
1 b
1 1
-
u = ln(x) dv = dx
x21 x x)
ln(
=− − -
du = 1dx v=− 1 x 5 6
− −
x x ln(b) 11 1 − 1
ln(1)
= − −
b b
1
Now we take a
limit, ln(b) 1
=− +1 −
5 b b 5
ˆ ∞
ln(b) +6 5 6 1/b
1− − =
lim 6 b) − 0 + L ’ H
ln( −
dx = lim ˆ b
− = +1=0+1=
2
1 1 lim 1
2
dx = lim
1 x b→∞ 1 x b→∞ b b b→∞ b b→∞ 1
L’Hôpital’s Rule Suppose that f (a) = g(a) = 0, that f (x) and g(x) are differentiable on an open interval I containing a
f ( x) f'(x)
lim = lim
,
x→a g(x) x→a g'(x)
assuming that the limit on the left and right both exist.
Example 2: Evaluate
ˆ ∞
1
dx.
−∞ 1 + x2
According to part 3 of our definition, we can choose any real number c and split this integral into two
integrals and then apply parts 1 and 2 to each piece. Let’s choose c = 0 and write
∞
ˆ 1 0
ˆ 1 ∞
1
dx ˆ
−∞ 2
1+x = −∞ dx +
1 + x2 dx.
0 1 + x2
Now we will evaluate each piece separately.
0
ˆ 0
1 b
1 ˆ
∞ 1
dx = lim ˆ 1 dx = lim
ˆ
2
2
dx 2 2
dx
-0 -b
= −1
lim = lim tan−1(x)
tan (x) -
a→−∞ - b→∞ -
1 0
ˆ ∞ ˆ 0 ˆ ∞
1 1 1 π π
= + = + =
2 dx dx π
−∞ 1+x −∞ 1 + x2 0 1 + x2 2 2
Since 1/(1 + x2) > 0 on R, the improper integral can be interpreted as the (finite) area between the curve
and the x-axis.
y
y=
1 + x2 Area = π
x
0
A Special Example: For what values of p does the integral
ˆ ∞
1
dx
1 xp
converge? When the integral does converge, what is its
value?
If p /= 1:
If p = 1:
ˆ ∞ ˆ b
1
= lim x−p dx
dx ˆ b
ˆ 1
1 b→∞ 1
∞
= lim
1
xp
b dx dx
= lim +1 - x x
x−p 1 b→∞ 1
b→∞ −p + 1 - x-
-b
1 = lim )
-b ln(b→∞ -1
1 1
=
· = lim [ln(b)
lim 1 − p xp−1 -
5 1
−
46
= lim 1
1
1 1 − 1 =
3 p− , p>1 = lim ln(b) =∞
b→∞
1 p−
1
Combining ˆ ∞ 1
1 , p>1
dx =
1
∞,
xp p≤1
Integrands with Vertical Asymptotes: Another type of improper integral that can arise is when the
integrand has a vertical asymptote (infinite discontinuity) at a limit of integration or at a point on the
interval of integration. We apply a similar technique as in the previous examples of integrating over an
altered interval before obtaining the integral we want by taking limits.
ˆ 1 1
1 √ =2− √
1
dx = 2x 1/2 =2 x √1 a = 2(1 − a).
√ dx = ˆ
a x −1/2
- 2
x -a a
1 √ = 2. !
lim ˆ 1 a
a→0+ 1 dx = lim 2 − "
√
Therefore, a x a→0+
1
1 2
ˆ 1
1
√ ˆ
dx = lim √ dx =
+
0 x a→0 a x
Definition: Integrals of functions that become infinite at a point within the interval of integration are called improper integrals of Type II.
ˆ b
a
ˆ
f (x) dx = lim f (x) dx.
a c→a+
c
a c→b− a
3. If f (x) is discontinuous at c, where a < c < b, and continuous on [a, c) ∪ (c, b], then
ˆ b
c b
ˆ ˆ
f (x) dx = f (x) dx + f (x) dx.
a
a c
In each case, if the limit is finite we sat that the improper integral converges and that the
limit is the
value of the improper integral. If the limit fails to exist, the improper integral
diverges
0 1 −x
ˆ 1 ˆ b
1 1
dx = lim dx
1 − 0 1−x
b→1−
0 x ˆ b
= 1
lim
− dx
b→1− 0 x− 1
-b
b→1−
= lim
b→1− -
= lim
b→1−
= lim
− ln |x − 1| − ln(x − 1)
0
0
− ln(1 − b)
b
= − (−∞) -
=
∞
Tests for Convergence: When we cannot evaluate an improper integral directly, we try to determine
whether it converges of diverges. If the integral diverges, we are done. If it converges we can use
numerical methods to approximate its value. The principal tests for convergence or divergence are the
Direct Comparison Test and the Limit Comparison Test.
Direct Comparison Test for Integrals: If 0 ≤ f (x) ≤ g(x) on the interval (a, ∞], where a ∈ R, then,
ˆ ∞ ˆ∞
1. If ( ) converges, then so does ( ) .
g x dx f x dx
a a
ˆ ∞ ˆ ∞
2. If ( ) diverges, then so does ( ) .
f x dx g x dx
a a
2. If the area under the curve of f (x) is infinite and g(x) is bounded below by f (x), then the area under
the curve of g(x) must be “less than or equal to” the area under the curve of g(x). Since there is no finite
number “greater than” infinity, the area under g(x) must also be infinite.
dx.
2 x2
We want to find a
function cos2(x)
,cos(2(x)
)= ( )=
( ) such that for
some ( )
or ( ) for all .
gx a∈R f x ≤g x f x ≥g x x≥a
x2 x2
2
One way we can do this is by finding bounds for f (x). Since 0 ≤ cos (x) ≤ 1 for all x,
cos2(x) 1
≤ .
x2 x2
1
So then we can use g(x) :=
2
. So,
x
ˆ ∞ ˆ ∞ ˆ b 3 3 44
cos2(x) 1 1 1 1 1
0≤ = lim = lim =
dx ≤ dx b→∞ dx b→∞ − .
2 x2 2 x 2 2 x2 −b − 2
ˆ ∞ 2 2
cos (x)
So converges.
dx
2 x2
Example 7: Determine if the following integral is convergent of divergent.
ˆ ∞
1
dx.
3 x − e−x
1 1
Since x ≥ x − e−x, f (x) := ≤ =: g(x) for all x ≥ 3. So,
x x − e−1
ˆ ∞ ˆ ∞
0 ( ) ( )
≤
g x dx.
f x dx ≤
3 3
ˆ ∞ ˆ ∞
1 1
By the Direct Comparison Test then, diverges since diverges.
dx dx
3 x − e−x 3 x
Limit Comparison Test for Integrals: If the positive functions f (x) and g(x) are continuous on [a, ∞), and if
f (x)
lim = L, 0 < L < ∞,
x→∞ g(x)
then
ˆ ∞ ∞
ˆ
f (x) dx and g(x) dx
both converge or
diverge. a a
Why does this make sense? The convergence is really only dependent on the “tail” of the integral. That
is, the convergence is dictated by what happens “at infinity.” If for sufficiently large values of x, f (x) ≈
Lg(x) and one of the two integrals converges, then the other one should also converge, since it is only
off by “about a scalar multiple.” The same goes for diverging, if one diverges, then multiplying it by a
positive number won’t suddenly make it converge, so the other one should also diverge.
ˆ ∞
1
dx
1 1 + x2
converges.
1
Let
1 f (x) := and g(x) :=
. Then,
1 + x2 x2
2
+ x2 − 1 4
lim f (x) x 1 = lim 1
= lim = lim − =1
x→∞ g(x) x→∞ 1 + 3 1+ .
x2 x→∞ 1 x→∞ x2
ˆ 1 + x2
So,
∞ by the Limit Comparison Test, the integral
1
converges.