MIT OpenCourseWare http://ocw.mit.

edu

18.01 Single Variable Calculus

Fall 2006

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.

Lecture 26

18.01 Fall 2006

Lecture 26: Trigonometric Integrals and Substitution

Trigonometric Integrals
How do you integrate an expression like sinnx cosmx dx? (n = 0, 1, 2... and m = 0, 1, 2, . . .)

We already know that: sin x dx = cos x + c

and cos x dx = sin x + c

Method A
Suppose either n or m is odd. Example 1. sin3x cos2x dx. Our strategy is to use sin2x + cos2x = 1 to rewrite our integral in the form: 3 2 sin x cos x dx = f (cosx) sinx dx Indeed, sin3x cos2x dx = sin2x cos2x sin x dx = (1 cos2x) cos2x sin x dx

Next, use the substitution u = cos x and du = sin x dx Then, = (1 cos2 x) cos2 x sin x dx = (1 u2 )u2 (du)

1 1 1 1 (u2 + u4 )du = u3 + u5 + c = cos3 u + cos5 x + c 5 5 3 3

Example 2. cos3 x dx = f (sin x) cos x dx = (1 sin2 x) cos x dx

Again, use a substitution, namely u = sin x and du = cos x dx u3 sin3 x cos3 x dx = (1 u2 )du = u + c = sin x +c 3 3

Lecture 26

18.01 Fall 2006

Method B
This method requires both m and n to be even. It requires double-angle formulae such as cos2 x = 1 + cos 2x 2

(Recall that cos 2x = cos2 x sin2 x = cos2 x (1 sin2 x) = 2 cos2 x 1) Integrating gets us 1 + cos 2x x sin(2x) 2 cos x dx = dx = + +c 2 2 4 We follow a similar process for integrating sin2 x. sin2 x = sin x dx =
2

1 cos(2x) 2

1 cos(2x) x sin(2x) dx = +c 2 2 4

The full strategy for these types of problems is to keep applying Method B until you can apply Method A (when one of m or n is odd). Example 3. sin2 x cos2 x dx.

Applying Method B twice yields 1 cos 2x 1 + cos 2x 1 1 2 cos 2x dx dx = 2 2 4 4 1 1 1 1 (1 + cos 4x) dx = x = sin 4x + c 4 8 8 32 There is a shortcut for Example 3. Because sin 2x = 2 sin x cos x, sin x cos x dx =
2 2

1 sin 2x 2

1 dx = 4

1 cos 4x dx = same as above 2

The next family of trig integrals, which well start today, but will not nish is: secn x tanm x dx where n = 0, 1, 2, . . . and m = 0, 1, 2, . . . Remember that sec2 x = 1 + tan2 x which we double check by writing 1 sin2 x cos2 x + sin2 x = 1 + = cos2 x cos2 x cos3 x sec2 x dx = tan x + c sec x tan x dx = sec x + c

Lecture 26

18.01 Fall 2006

To calculate the integral of tan x, write sin x tan x dx = dx cos x Let u = cos x and du = sin x dx, then sin x du tan x dx = dx = = ln(u) + c cos x u tan x dx = ln(cos x) + c (Well gure out what sec x dx is later.)

Now, lets see what happens when you have an even power of secant. (The case n even.) sec4 x dx = f (tanx) sec2 x dx = (1 + tan2 x) sec2 x dx Make the following substitution: u = tan x and du = sec2 x dx sec4 x dx = (1 + u2 )du = u + u3 tan3 x + c = tan x + +c 3 3

What happens when you have a odd power of tan? (The case m odd.) tan3 x sec x dx = f (sec x) d(sec x) = (sec2 x 1) sec x tan x dx

(Remember that sec2 x 1 = tan2 x.) Use substitution: u = sec x and du = sec x tan x dx Then, tan x secx dx =
3

(u2 1)du =

u3 sec3 x sec x + c u+c= 3 3

We carry out one nal case: n = 1, m = 0 sec x dx = ln (tan x + sec x) + c

Lecture 26

18.01 Fall 2006

We get the answer by advanced guessing, i.e., knowing the answer ahead of time. sec x + tan x sec2 x + sec x tan x sec x dx = sec x dx = dx sec x + tan x tan x + sec x Make the following substitutions: u = tan x + sec x and du = (sec2 x + sec x tan x) dx This gives sec x dx = du = ln(u) + c = ln(tan x + sec x) + c u

Cases like n = 3, m = 0 or more generally n odd and m even are more complicated and will be discussed later.

Trigonometric Substitution
Knowing how to evaluate all of these trigonometric integrals turns out to be useful for evaluating integrals involving square roots. Example 4. y = a2 x2

Figure 1: Graph of the circle x2 + y2 = a2 . We already know that the area of the top half of the disk is a a2 a2 x2 dx = 2 a 4

Lecture 26

18.01 Fall 2006

What if we want to nd this area?

Figure 2: Area to be evaluated is shaded. To do so, you need to evaluate this integral: t=x a2 t2 dt
t=0

Let t = a sin u and dt = a cos u du. (Remember to change the limits of integration when you do a change of variables.) Then, a2 t2 = a2 a2 sin2 u = a2 cos2 u; a2 t2 = a cos u Plugging this into the integral gives us
0 x

a2 t2 dt = (a cos u) a cos u du = a2

u=sin1 (x/a)

cos2 u du

u=0

Heres how we calculated the new limits of integration: t = t =

0 x

0 = a sin u = 0 = u = 0 x = a sin u = x = u = sin1 (x/a)

0 sin1 (x/a)

a2 t2 dt = a2 = a sin
2 1

cos2 u du = a2 + a 4
2

u sin 2u + 2 4

1 sin (x/a)
0

(x/a)

2 sin(sin1 (x/a)) cos(sin1 (x/a))

(Remember, sin 2u = 2 sin u cos u.) Well pick up from here next lecture (Lecture 28 since Lecture 27 is Exam 3). 5

Lecture 28

18.01 Fall 2006

Lecture 28: Integration by Inverse Substitution;

Completing the Square

Trigonometric Substitutions, continued

-a
x

Figure 1: Find area of shaded portion of semicircle. a2 t2 dt

dt = a cos u du a2 t2 = a2 a2 sin2 u = a2 cos2 u = a2 t2 = a cos u (No more square root!) Start: x = a u = /2; Finish: x = a u = /2 1 + cos(2u) u sin(2u) a2 t2 dt = a2 cos2 u du = a2 du = a2 + +c 2 2 4 1 + cos(2u) ). 2 We want to express this in terms of x, not u. When t = 0, a sin u = 0, and therefore u = 0. When t = x, a sin u = x, and therefore u = sin1 (x/a). (Recall, cos2 u = sin(2u) 2 sin u cos u 1 = = sin u cos u 4 4 2 x sin u = sin sin1 (x/a) = a

t = a sin u;

Lecture 28

18.01 Fall 2006

How can we nd cos u = cos sin1 (x/a) ? Answer: use a right triangle (Figure 2).

a u a-x
Figure 2: sin u = x/a; cos u = From the diagram, we see cos u = And nally,
0 x

p a2 x2 /a.

a 2 x2 a

a2

t2

dt = a

1 u 1 (x/a) 1 x a2 x2 2 sin + sin u cos u 0 = a + 4 2 2 2 a a a2 x 1 a2 t2 dt = sin1 ( ) + x a2 x2 2 a 2

When the answer is this complicated, the route to getting there has to be rather complicated. Theres no way to avoid the complexity. 1 Lets double-check this answer. The area of the upper shaded sector in Figure 3 is a2 u. The 2 2 2 area of the lower shaded region, which is a triangle of height a2 x2 and base x, is 1 2x a x .

Lecture 28

18.01 Fall 2006

Figure 3: Area divided into a sector and a triangle. Here is a list of integrals that can be computed using a trig substitution and a trig identity. integral dx x2 + 1 dx x2 1 dx 1 x2 substitution x = tan u x = sec u x = sin u trig identity tan2 u + 1 = sec2 u sec2 u 1 = tan2 u 1 sin2 u = cos2 u

Lets extend this further. How can we evaluate an integral like this? dx 2 x + 4x When you have a linear and a quadratic term under the square root, complete the square. x2 + 4x = (something)2 constant In this case, (x + 2)2 = x2 + 4x + 4 = x2 + 4x = (x + 2)2 4 Now, we make a substitution. v =x+2 Plugging these in gives us dx = (x + 2)2 4 dv v2 4 and dv = dx

Now, let v = 2 sec u and dv = 2 sec u tan u dv 2 sec u tan u du = = sec u du 2 tan u v2 4 3

Lecture 28

18.01 Fall 2006

Remember that

sec u du = ln(sec u + tan u) + c

Finally, rewrite everything in terms of x. v = 2 sec u cos u = 2 v

Set up a right triangle as in Figure 4. Express tan u in terms of v .

v u 2

v-4

Figure 4: sec u = v/2 or cos u = 2/v. Just from looking at the triangle, we can read o v2 4 and tan u = 2 v v2 4 2 sec u du = ln + +c 2 2 = ln(v + v 2 4) ln 2 + c v sec u = 2 We can combine those last two terms into another constant, c . dx = ln(x + 2 + x2 + 4x) + c x2 + 4x Heres a teaser for next time. In the next lecture, well integrate all rational functions. By rational functions, we mean functions that are the ratios of polynomials: P (x)
Q(x)
Its easy to evaluate an expression like this: 1 3 + dx = ln |x 1| + 3 ln |x + 2| + c x1 x+2 4

Lecture 28

18.01 Fall 2006

If we write it a bit dierently, however, it becomes much harder to integrate: 1 3 x + 2 + 3(x 1) 4x 1

+ = = 2 x1 x+2 (x 1)(x + 2) x + x 2
4x 1 = ??? x2 + x 2
How can we reorganize what to do starting from (4x 1)/(x2 + x 2)? Next time, well see how. It involves some algebra.

Lecture 29

18.01 Fall 2006

Lecture 29: Partial Fractions

We continue the discussion we started last lecture about integrating rational functions. We dened a rational function as the ratio of two polynomials: P (x)
Q(x)
We looked at the example

1 3 + dx = ln |x 1| + 3 ln |x + 2| + c x1 x+2

That same problem can be disguised: 1 3 (x + 2) + 3(x 1) 4x 1

+ = = 2 x1 x+2 (x 1)(x + 2) x + x 2
which leaves us to integrate this: 4x 1 dx = ??? x2 + x 2 P (x) into simpler pieces. Q(x)

Goal: we want to gure out a systematic way to split First, we factor the denominator Q(x).

4x 1 4x 1 A B
= = + x2 + x 2 (x 1)(x + 2) x 1 x + 2

Theres a slow way to nd A and B . You can clear the denominator by multiplying through by (x 1)(x + 2): (4x 1) = A(x + 2) + B (x 1) From this, you nd 4=A+B and 1 = 2A B You can then solve these simultaneous linear equations for A and B . This approach can take a very long time if youre working with 3, 4, or more variables.

Theres a faster way, which we call the cover-up method. Multiply both sides by (x 1): 4x 1 B =A+ (x 1) x+2 x+2 Set x = 1 to make the B term drop out: 41 =A 1+2 A=1 1

Lecture 29

18.01 Fall 2006

The fastest way is to do this in your head or physically cover up the struck-through terms. For instance, to evaluate B : 4x 1 A B = + (x 1) (x + 2) x1 (x + 2) Implicitly, we are multiplying by (x + 2) and setting x = 2. This gives us 4(2) 1 =B 2 1 = B=3

What weve described so far works when Q(x) factors completely into distinct factors and the degree of P is less than the degree of Q.

If the factors of Q repeat, we use a slightly dierent approach. For example:

A B C
x2 + 2 = + + 2 2 (x 1) (x + 2) x 1 (x 1) x+2 Use the cover-up method on the highest degree term in (x 1). x2 + 1 2 = B + [stu](x 1) x+2 = 12 + 2 =B 1+2 = B=1

Implicitly, we multiplied by (x 1)2 , then took the limit as x 1. C can also be evaluated by the cover-up method. Set x = 2 to get x2 + 2 = C + [stu](x + 2) (x 1) This yields x2 + 2 A 1 2/3 = + + 2 2 (x 1) (x + 2) x 1 (x 1) x+2 Cover-up cant be used to evaluate A. Instead, plug in an easy value of x: x = 0. 2 A 1 1 1 = +1+ = 1 = 1 + A = A = 2 (1) (2) 1 3 3 3 Now we have a complete answer: x2 + 2 1 1 2 = + + (x 1)2 (x + 2) 3(x 1) (x 1)2 3(x + 2)
2

(2)2 + 2 =C (2 1)2

C=

2 3

Not all polynomials factor completely (without resorting to using complex numbers). For exam ple: 1 A1 B1 x + C1 = + 2 (x + 1)(x 1) x1 x2 + 1 We nd A1 , as usual, by the cover-up method. 12 1 = A1 +1 2 = A1 = 1 2

Lecture 29

18.01 Fall 2006

Now, we have 1 1/2 B1 x + C1 = + (x2 + 1)(x 1) x1 x2 + 1 Plug in x = 0. 1 1 C1 = + 1(1) 2 1 = C1 = 1 2

Now, plug in any value other than x = 0, 1. For example, lets use x = 1. 1 1/2 B1 (1) 1/2 B1 1/2 1 = + = 0 = = B1 = 2(2) 2 2 2 2 Alternatively, you can multiply out to clear the denominators (not done here). Lets try to integrate this function, now. dx 1 dx 1 x dx 1 dx = 2 2 2 (x + 1)(x 1) 2 x1 2 x +1 2 x +1 = 1 1 1 ln |x 1| ln | x2 + 1 | tan1 x + c 2 4 2

What if were faced with something that looks like this? dx (x 1)10 This is actually quite simple to integrate: dx 1 = (x 1)9 + c (x 1)10 9

What about this?

dx (x2 + 1)10

Here, we would use trig substitution: x = tan u and and the trig identity tan2 u + 1 = sec2 u to get sec2 u du = (sec2 u)10 cos18 u du dx = sec2 udu

From here, we can evaluate this integral using the methods we introduced two lectures ago.

Lecture 30

18.01 Fall 2006

Lecture 30: Integration by Parts, Reduction Formulae

Integration by Parts
Remember the product rule: (uv ) = u v + uv We can rewrite that as uv = (uv ) u v Integrate this to get the formula for integration by parts: uv dx = uv u v dx

Example 1.

tan1 x dx.

At rst, its not clear how integration by parts helps. Write tan1 x dx = tan1 x(1 dx) = uv dx with u = tan1 x and Therefore, 1 1 + x2 Plug all of these into the formula for integration by parts to get: tan1 x dx = uv dx = (tan1 x)x v = x and u = = x tan1 x v = 1.

1 (x)dx 1 + x2

1 ln |1 + x2 | + c 2

Alternative Approach to Integration by Parts

As above, the product rule: (uv ) = u v + uv can be rewritten as uv = (uv ) u v This time, lets take the denite integral:
a b

uv dx =
a

(uv ) dx
a

u v dx

Lecture 30

18.01 Fall 2006

By the fundamental theorem of calculus, we can say b b uv dx = uv

a a

u v dx

Another notation in the indenite case is u dv = uv v du This is the same because dv = v dx = uv dx = u dv Example 2. (ln x)dx 1 u = ln x; du = dx and dv = dx; v = x x 1 (ln x)dx = x ln x x dx = x ln x dx = x ln x x + c x We can also use advanced guessing to solve this problem. We know that the derivative of something equals ln x: d (??) = ln x dx Lets try d 1 (x ln x) = ln x + x = ln x + 1 dx x Thats almost it, but not quite. Lets repair this guess to get: d (x ln x x) = ln x + 1 1 = ln x dx and du = u dx = u v dx = vu dx = v du

Reduction Formulas (Recurrence Formulas)

Example 3. Lets try: u = (ln x)n = u = n(ln x)n1 v = dx; v = x Plugging these into the formula for integration by parts gives us: (ln x) dx = x(ln x)
n n

(ln x) dx 1 x

n(ln x)

n1

1 1 x dx x

Keep repeating integration by parts to get the full formula: n (n 1) (n 2) (n 3) etc Example 4. xn ex dx Lets try: u = xn = u = nxn1 ; 2 v = ex = v = ex

Lecture 30

18.01 Fall 2006

Putting these into the integration by parts formula gives us: xn ex dx = xn ex nxn1 ex dx Repeat, going from n (n 1) (n 2) etc.

Bad news: If you change the integrals just a little bit, they become impossible to evaluate: 1 2 tan x dx = impossible ex dx = also impossible x

Good news: When you cant evaluate an integral, then 2 x e dx 1 x is an answer, not a question. This is the solution you dont have to integrate it! The most important thing is setting up the integral! (Once youve done that, you can always evaluate it numerically on a computer.) So, why bother to evaluate integrals by hand, then? Because you often get families of related integrals, such as x e F (a) = dx a x 1 where you want to nd how the answer depends on, say, a.

Lecture 30

18.01 Fall 2006

Arc Length
This is very useful to know for 18.02 (multi-variable calculus).

y y=f(x)
ds dx dy

Figure 1: Innitesimal Arc Length ds

ds dx

dy

Figure 2: Zoom in on Figure 1 to see an approximate right triangle. In Figures 1 and 2, s denotes arc length and ds = the innitesmal of arc length. 2 ds = (dx)2 + (dy )2 = 1 + (dy/dx) dx Integrating with respect to ds nds the length of a curve between two points (see Figure 3). To nd the length of the curve between P0 and P1 , evaluate:
P1

ds
P0

Lecture 30

18.01 Fall 2006

P P

Figure 3: Find length of curve between P0 and P1 . We want to integrate with respect to x, not s, so we do the same algebra as above to nd ds in terms of dx. 2 (ds)2 (dx)2 (dy )2 dy = + =1+ (dx)2 (dx)2 (dx)2 dx Therefore,
P1

ds =
a

1+

P0

dy dx

2 dx

Example 5: The Circle. x2 + y 2 = 1 (see Figure 4).

Figure 4: The circle in Example 1.

Lecture 30

18.01 Fall 2006

We want to nd the length of the arc in Figure 5:

a
Figure 5: Arc length to be evaluated.

1 x2 dy 2x 1 x = = 2 dx 2 1x 1 x2 2 x ds = 1 + dx 1 x2 2 x x2 1 x2 + x2 1 1+ =1+ = = 2 2 2 1x 1x 1 x2 1x 1 ds = dx 1 x2 a a dx s= = sin1 x = sin1 a sin1 0 = sin1 a 2 0 1x 0 sin s = a y= This is illustrated in Figure 6.

Lecture 30

18.01 Fall 2006

s a 1 a

Figure 6: s = angle in radians.

Parametric Equations
Example 6. x = a cos t y = a sin t Ask yourself: whats constant? Whats varying? Here, t is variable and a is constant. Is there a relationship between x and y ? Yes: x2 + y 2 = a2 cos2 t + a2 sin2 t = a2 Extra information (besides the circle): At t = 0, x = a cos 0 = a and y = a sin 0 = 0 At t = , 2 x = a cos = 0 and y = a sin = a 2 2 Thus, for 0 t /2, a quarter circle is traced counter-clockwise (Figure 7).

Lecture 30

18.01 Fall 2006

t=/2 (0,a)

(a,0) t=0
Figure 7: Example 6. x = a cos t, y = a sin t; the particle is moving counterclockwise. Example 7: The Ellipse See Figure 8. x = 2 sin t;
2

y = cos t

x + y 2 = 1( = (2 sin t)2 /4 + (cos t)2 = sin2 t + cos2 t = 1) 4

t=0 (0,1)

(2,0) t=/2
Figure 8: Ellipse: x = 2 sin t, y = cos t (traced clockwise). Arclength ds for Example 6. dx = a sin t dt, dy = a cos t dt ds = (dx)2 + (dy )2 = (a sin t dt)2 + (a cos t dt)2 = (a sin t)2 + (a cos t)2 dt = a dt

Lecture 31

18.01 Fall 2006

Lecture 31: Parametric Equations, Arclength, Surface Area

Arclength, continued
Example 1. Consider this parametric equation: x = t2 x3 = (t2 )3 = t6 ; y = t3 for 0 t 1 0x1

y 2 = (t3 )2 = t6

= x3 = y 2 = y = x2/3

ds
ds dy dx

dy

dx

Figure 1: Innitesimal Arclength. (ds)2 = (dx)2 + (dy )2 (2t dt)2 + (3t2 dt)2 = (4t2 + 9t4 )(dt)2 (dx)2 (dy )2 t=1 1 1 2 4 Length = ds = 4t + 9t dt = t 4 + 9t2 dt (ds)2 =
t=0 0 2 3/2 1 0

(4 + 9t ) 1 (133/2 43/2 ) = 27 27 0 Even if you cant evaluate the integral analytically, you can always use numerical methods. =

Lecture 31

18.01 Fall 2006

Surface Area (surfaces of revolution)

y ds y x

a
Figure 2: Calculating surface area

ds (the innitesimal curve length in Figure 2) is revolved a distance 2y . The surface area of the thin strip of width ds is 2y ds. Example 2. Revolve Example 1 (x = t2 , y = t3 , 0 t 1) around the x-axis. Refer to Figure 3.

Figure 3: Curved surface of a trumpet.

Lecture 31

18.01 Fall 2006

Area = 2y ds =

t3 y

1 t 4 + 9t2 dt t4 4 + 9t2 dt = 2 0 ds

Now, we discuss the method used to evaluate t4 (4 + 9t2 )1/2 dt Were going to ignore the factor of 2 . You can reinsert it once youre done evaluating the integral. We use the trigonometric substitution t= 2 tan u; 3 dt = 2 sec2 u du; 3 tan2 u + 1 = sec2 u

Putting all of this together gives us: t4 (4 + 9t2 )1/2 dt = = 4 1/2 2 4 2 tan u 4+9 tan2 u sec2 u du 3 9 3 5 2 tan4 u(2 sec u)(sec2 u du) 3

This is a tan sec integral. Its doable, but it will take a long time for you to work the whole thing out. Were going to stop evaluating it here.

Example 3 Lets use what weve learned to nd the surface area of the unit sphere (see Figure 4).

rotate the curve by 2 radians

. .

Figure 4: Slice of spherical surface (orange peel, only, not the insides).

Lecture 31

18.01 Fall 2006

For the top half of the sphere, y= 1 x2 We want to nd the area of the spherical slice between x = a and x = b. A spherical slice has area A=
x=a x=b

2y ds

From last time, dx ds = 1 x2 Plugging that in yields a remarkably simple formula for A: A=
a b

1 x2

dx = 1 x2

2 dx

= 2 (b a) Special Cases For a whole sphere, a = 1, and b = 1. 2 (1 (1)) = 4 is the surface area of a unit sphere. For a half sphere, a = 0 and b = 1. 2 (1 0) = 2

Lecture 32

18.01 Fall 2006

Lecture 32: Polar Co-ordinates, Area in Polar

Co-ordinates

Polar Coordinates

Figure 1: Polar Co-ordinates. In polar coordinates, we specify an objects position in terms of its distance r from the origin and the angle that the ray from the origin to the point makes with respect to the x-axis. Example 1. What are the polar coordinates for the point specied by (1, 1) in rectangular coordinates?

r
(1,-1)
Figure 2: Rectangular Co-ordinates to Polar Co-ordinates. 2

= =

12 + (1)2 = 4

In most cases, we use the convention that r 0 and 0 2 . But another common convention is to say r 0 and . All values of and even negative values of r can be used.

Lecture 32

18.01 Fall 2006

r
x

Figure 3: Rectangular Co-ordinates to Polar Co-ordinates.

Regardless of whether we allow positive or negative values of r or , what is always true is:
x = r cos and y = r sin

3 For instance, x = 1, y = 1 can be represented by r = 2, = : 4 3 1 = x = 2 cos 4 and 3 1 = y = 2 sin 4

Example 2. Consider a circle of radius a with its center at x = a, y = 0. We want to nd an equation that relates r to .

(a,0)

Figure 4: Circle of radius a with center at x = a, y = 0.

Lecture 32

18.01 Fall 2006

We know the equation for the circle in rectangular coordinates is (x a)2 + y 2 = a2 Start by plugging in: x = r cos This gives us (r cos a)2 + (r sin )2 = a2
r2 cos2 2arcos + a2 + r2 sin2 = a2
r2 2ar cos = 0
r = 2a cos The range of 0 traces out the top half of the circle, while 0 traces out the bottom 2 2 half. Lets graph this. and y = r sin

= /4 r (a,0) =0 x

Figure 5: r = 2a cos , /2 /2.

At = 0, r = 2a = x = 2a, y = 0 At = , r = 2a cos = a 2 4 4 The main issue is nding the range of tracing the circle once. In this case, = = 2 (down) 2 (up) << . 2 2

3 << . When = , r = 2a cos = 2a(1) = 2a. The 2 2 3 radius points backwards. In the range < < , the same circle is traced out a second time. 2 2 Weird range (avoid this one):

Lecture 32

18.01 Fall 2006

r=f()

Figure 6: Using polar co-ordinates to nd area of a generic function.

Area in Polar Coordinates

Since radius is a function of angle (r = f ()), we will integrate with respect to . The question is: what, exactly, should we integrate? 2 ?? d
1

Lets look at a very small slice of this region:

rd

Figure 7: Approximate slice of area in polar coordinates. This innitesimal slice is approximately a right triangle. To nd its area, we take: Area of slice So, Total Area =
1 2

1 1 (base) (height) = r(r d) 2 2 1 2 r d 2

Lecture 32 < < (the circle in Figure 5). 2 2 /2 /2 1 A = area = (2a cos )2 d = 2a2 cos2 d /2 2 /2

18.01 Fall 2006

Example 3. r = 2a cos , and

Because cos2 =

1 1 + cos 2, we can rewrite this as 2 2 /2 A = area = (1 + cos 2) d = a2

/2

/2

d + a2

/2

cos 2 d
/2

/2

= a2 +

/2 0 1 1 sin 2 = a2 + [sin sin( )] 2 2 /2 A = area = a2

Example 4: Circle centered at the Origin.

r=a

Figure 8: Example 4: Circle centered at the origin x = r cos ; y = r sin x + y 2 = r2 cos2 + r2 sin2 = r2
2

The circle is x2 + y 2 = a2 , so r = a and x = a cos ; y = a sin A=

0 2

1 2 1 a d = a2 2 = a2 . 2 2

Lecture 32

18.01 Fall 2006

Example 5: A Ray. In this case, = b.

=b

Figure 9: Example 5: The ray = b, 0 r < .

The range of r is 0 r < ; x = r cos b; y = r sin b.

Example 6: Finding the Polar Formula, based on the Cartesian Formula

1/sin

Figure 10: Example 6: Cartesian Form to Polar Form Consider, in cartesian coordinates, the line y = 1. To nd the polar coordinate equation, plug in y = r sin and x = r cos and solve for r. r sin = 1 = r = 1 sin with 0<<

Lecture 32

18.01 Fall 2006

Example 7: Going back to (x, y ) coordinates from r = f (). Start with 1 r= . 1 1 + 2 sin Hence, r+ Plug in r = x2 + y 2 : Finally, 3y 2 +y =1 4 This is an equation for an ellipse, with the origin at one focus. x2 + Useful conversion formulas: r= x2 + y 2 and = tan1 y x x2 + y 2 = 1 y 2 r sin = 1 2 y =1 2 y 2 y2 x2 + y 2 = 1 =1y+ 2 4

x2 + y 2 +

Example 8: A Rose r = cos(2) The graph looks a bit like a ower:

r<0

/4

r>0 r>0 1

r<0
Figure 11: Example 8: Rose For the rst petal << 4 4

-/4

Note: Next lecture is Lecture 34 as Lecture 33 is Exam 4. 7

Lecture 32: Exam 4 Review

18.01 Fall 2006

Exam 4 Review
1. Trig substitution and trig integrals. 2. Partial fractions. 3. Integration by parts. 4. Arc length and surface area of revolution 5. Polar coordinates 6. Area in polar coordinates.

Questions from the Students

Q: What do we need to know about parametric equations? A: Just keep this formula in mind: ds = Example: Youre given x(t) = t4 and y (t) = 1 + t Find s (length). ds = Then, integrate with respect to t. Q: Can you quickly review how to do partial fractions? A: When nding partial fractions, rst check whether the degree of the numerator is greater than or equal to the degree of the denominator. If so, you rst need to do algebraic longdivision. If not, then you can split into partial fractions. Example. x2 + x + 1
(x 1)2 (x + 2)
We already know the form of the solution: x2 + x + 1 A B C = + + (x 1)2 (x + 2) x 1 (x 1)2 x+2 There are two coecients that are easy to nd: B and C . We can nd these by the cover-up method. 12 + 1 + 1 3 B= = (x 1) 1+2 3 1
(4t3 )2 + (1)2 dt
dx dt 2 + dy dt 2

Lecture 32: Exam 4 Review

18.01 Fall 2006

To nd C , C= (2)2 2 + 1 1 = (2 1)2 3 (x 2)

To nd A, one method is to plug in the easiest value of x other than the ones we already used (x = 1, 2). Usually, we use x = 0. 1 A 1 1/3 = + + (1)2 (2) 1 (1)2 2 and then solve to nd A. The Review Sheet handed out during lecture follows on the next page.

Lecture 32: Exam 4 Review

18.01 Fall 2006

Exam 4 Review Handout

1. Integrate by trigonometric substitution; evaluate the trigonometric integral and work backwards to the original variable by evaluating trig(trig1 ) using a right triangle: a) a2 x2 use x = a sin u, dx = a cos u du. b) a2 + x2 use x = a tan u, dx = a sec2 u du
c) x2 a2 use x = a sec u, dx = a sec u tan u du
2. Integrate rational functions P/Q (ratio of polynomials) by the method of partial fractions: If the degree of P is less than the degree of Q, then factor Q completely into linear and quadratic factors, and write P/Q as a sum of simpler terms. For example, 3x2 + 1 A B1 B2 Cx + D = + + + 2 2 2 2 (x 1)(x + 2) (x + 9) x 1 (x + 2) (x + 2) x +9 Terms such as D/(x2 + 9) can be integrated using the trigonometric substitution x = 3 tan u. This method can be used to evaluate the integral of any rational function. In practice, the hard part turns out to be factoring the denominator! In recitation you encountered two other steps required to cover every case systematically, namely, completing the square1 and long division.2 3. Integration by parts:
a b

b b uv dx = uv u vdx a
a

This is used when u v is simpler than uv . (This is often the case if u is simpler than u.) 4. Arclength: ds = dx2 + dy 2 . Depending on whether you want to integrate with respect to x, t or y this is written ds = 1 + (dy/dx)2 dx; ds = (dx/dt)2 + (dy/dt)2 dt; ds = (dx/dy )2 + 1 dy 5. Surface area for a surface of revolution: a) around the x-axis: 2yds = 2y 1 + (dy/dx)2 dx (requires a formula for y = y (x)) b) around the y -axis: 2xds = 2x (dx/dy )2 + 1 dy (requires a formula for x = x(y )) 6. Polar coordinates: x = r cos , y = r sin (or, more rarely, r = x2 + y 2 , = tan1 (y/x)) a) Find the polar equation for a curve from its equation in (x, y ) variables by substitution. b) Sketch curves given in polar coordinates and understand the range of the variable (often in preparation for integration). 7. Area in polar coordinates:
2 1

1 2 r d 2

(Pay attention to the range of to be sure that you are not double-counting regions or missing them.)
example, we rewrite the denominator x2 + 4x + 13 = (x + 2)2 + 9 = u2 + a2 with u = x + 2 and a = 3. division is used when the degree of P is greater than or equal to the degree of Q. It expresses P (x)/Q(x) = P1 (x) + R(x)/Q(x) with P1 a quotient polynomial (easy to integrate) and R a remainder. The key point is that the remainder R has degree less than Q, so R/Q can be split into partial fractions.
2 Long 1 For

Lecture 32: Exam 4 Review

18.01 Fall 2006

The following formulas will be printed with Exam 4

sin2 x + cos2 x = 1; sin2 x = 1 1 cos 2x; 2 2 sec2 x = tan2 x + 1

cos2 x = 1 1 + cos 2x
2 2

cos 2x = cos2 x sin2 x;

sin 2x = 2 sin x cos x

d d tan x = sec2 x; sec x = sec x tan x; dx dx tan x dx = ln(cos x) + c;

d 1 d 1
tan1 x = ; sin1 x = dx 1 + x2 dx 1 x2 sec x dx = ln(sec x + tan x) + c

See the next page for a review on integration of rational functions.

Lecture 32: Exam 4 Review

18.01 Fall 2006

Postscript: Systematic integration of rational functions

For a general rational function P/Q, the rst step is to express P/Q as the sum of a polynomial and a ratio in which the numerator has smaller degree than the denominator. For example, x3 3x 2 =x+2+ 2 x2 2x + 1 x 2x + 1 (To carry out this long division, do not factor the denominator Q(x) = x2 2x + 1, just leave it alone.) The quotient x + 2 is a polynomial and is easy to integrate. The remainder term 3x 2 (x 1)2 has a numerator 3x 2 of degree 1 which is less than the degree 2 of the denominator (x 1)2 . Therefore there is a partial fraction decomposition. In fact, 3x 2 (3x 3) + 1 3 1 = = + (x 1)2 (x 1)2 x 1 (x 1)2 In general, if P has degree n and Q has degree m, then long division gives P (x) R(x) = P1 (x) + Q(x) Q(x) in which P1 , the quotient in the long division, has degree n m and R, the remainder in the long division, has degree at most m 1. Evaluation of the simple pieces The integral dx 1 = (x a)1n + c (x a)n n1

if n = 1 and ln |x a| + c if n = 1. On the other hand the terms xdx dx and (Ax2 + Bx + C )n (Ax2 + Bx + C )n are handled by rst completing the square: B2 Ax + Bx + C = A(x B/2A) + C 4A
2 2

Using the variable u =

A(x B/2A) yields combinations of integrals of the form udu du and 2 2 n 2 (u + k ) (u + k 2 )n

The rst integral is handled by the substitution w = u2 + k 2 , dw = 2udu. The second integral can be worked out using the trigonometric substitution u = k tan du = k sec2 d. This then leads to sec-tan integrals, and the actual computation for large values of n are long. There are also other cases that we will not cover systematically. Examples are below: 1. If Q(x) = (x a)m (x b)n , then the expression is A1 A2 Am B1 B2 Bn + + + + + + + x a (x a)2 x b (x b)2 (x a)m (x b)n 5

Lecture 32: Exam 4 Review

18.01 Fall 2006

2. If there are quadratic factors like (Ax2 + Bx + C )p , one gets terms Ax2 a1 x + b1 a2 x + b2 x ap x + bp + + + 2 2 2 + Bx + C (Ax + Bx + C ) (Ax + Bx + C )p

for each such factor. (To integrate these quadratic pieces complete the square and make a trigonometric substitution.)

Lecture 34

18.01 Fall 2006

Lecture 34: Indeterminate Forms - LHpitals Rule

LHpitals Rule
(Two correct spellings: LHpital and LHospital) Sometimes, we run into indeterminate forms. These are things like 0 0 and For instance, how do you deal with the following? x3 1 0 = ?? 2 x 1 x 1 0 lim

Example 0. One way of dealing with this is to use algebra to simplify things: x3 1 (x 1)(x2 + x + 1) x2 + x + 1 3 = lim = lim = 2 x1 x 1 x1 x1 (x 1)(x + 1) x+1 2 lim In general, when f (a) = g (a) = 0, f (x) f (x) f (a) lim f (x) f (a) xa xa lim = lim x a = = xa g (x) xa g (x) g (x) g (a) g (a) lim xa xa xa This is the easy version of LHpitals rule: f (x) f (a) = xa g (x) g (a) lim Note: this only works when g (a) = 0! In example 0, f (x) = x3 = 1; g (x) = x2 1
f (x) = 3x2 ; g (x) = 2x = f (1) = 3; g (1) = 2
The limit is f (1)/g (1) = 3/2. Now, lets go on to the full LHpital rule.

Lecture 34

18.01 Fall 2006

Example 1. Apply LHpitals rule (a.k.a. LHop) to x15 1 x 1 x3 1 lim to get x15 1 15x14 15 = lim = =5 3 2 x1 x 1 x1 3x 3 Lets compare this with the answer wed get if we used linear approximation techniques, instead of LHpitals rule: x15 1 15(x 1) lim (Here, f (x) = x15 1, a = 1, f (a) = b = 0, m = f (1) = 15, and f (x) m(x a) + b.) Similarly, x3 1 3(x 1) Therefore, x15 1 15(x 1) =5 x3 1 3(x 1)

Example 2. Apply LHop to

x 0

lim

sin 3x x

to get
x 0

lim

3 cos 3x =3 1

This is the same as

d sin(3x) = 3 cos(3x) =3 dx x=0 x=0

Example 3.
x 4

lim

sin x cos x cos x + sin x 1 1 = lim = + = 2 x 4 x 4 1 2 2 f (x) = sin x cos x, f (x) = cos x + sin x f = 2 4 y 0 are always a type of limit. x 0

Remark: Derivatives lim

x0

cos x 1 . x Use LHpitals rule to evaluate the limit: Example 4. lim

x 0

x 0

lim

cos x 1 sin x = lim =0 x0 x x

Lecture 34

18.01 Fall 2006

cos x 1 . x 0 x2 cos x 1 cos x 1 sin x cos x 1 lim = lim = lim = lim = x0 x0 x0 x0 x2 x2 2x 2 2 Just to check, lets compare that answer to the one we would get if we used quadratic approximation techniques. Remember that: 1 cos x 1 x2 (x 0) 2 1 2 1 1 x 1 ( )x2 cos x 1 1 2 2 = = x2 x2 2 x2 Example 5. lim

Example 6. lim

x 0

sin x . x2
x0

sin x cos x = lim By LHpitals rule x0 2x x2 If we apply LHpital again, we get sin x lim =0 x0 2 But this doesnt agree with what we get from taking the linear approximation: lim sin x x 1 2 = as x 0+ x2 x x We can clear up this seeming paradox by noting that lim cos x 1 = 2x 0

x0

0 The limit is not of the form , which means LHpitals rule cannot be used. The point is: look 0 before you LHp!

More interesting cases that work.

, or if x , or x . Lets apply this to rates of growth. Which function goes to faster: x, eax , or ln x? It is also okay to use LHpitals rule on limits of the form

Example 7. For a > 0,

So eax

eax aeax = lim = + x x x 1 grows faster than x (for a > 0). lim

Example 8. eax aeax c2 eax a10 eax = by LHpital = lim = lim = = lim = 10 9 8 x x x 10x x 10 9x x 10! lim

Lecture 34

18.01 Fall 2006

You can apply LHpitals rule ten times. Theres a better way, though: eax x10 1/10 = eax/10 x 10 = 10 =

eax = lim x x10 x lim

eax/10 x

Example 9.
x

lim

ln x 1/x = lim = lim 3x1/3 = 0 x 1/3x2/3 x x1/3 (x , a > 0)

Combining the preceding examples, ln x x1/3 x x10 eax LHpitals rule applies to

0 and . But, we sometimes face other indeterminate limits, such 0 0 as 1 , 0 , and 0 . Use algebra, exponentials, and logarithms to put these in LHpital form. Example 10. lim xx for x > 0.
x 0

Because the exponent is a variable, use base e:

x 0

lim xx = lim ex ln x
x 0

First, we need to evaluate the limit of the exponent

x0

lim x ln x 0 or . 0

This limit has the form 0 . We want to put it in the form Lets try to put it into the 0 form: 0

x 1/ ln x

We dont know how to nd lim

1 , though, so that approach isnt helpful. ln x Instead, lets try to put it into the form:

x0

ln x 1/x Using LHpitals rule, we nd lim x ln x = lim ln x 1/x = lim = lim (x) = 0 x0 1/x x0 1/x2 x0

x 0

Therefore,
x 0

lim (x ln x) lim xx = lim ex ln x = e x0 = e0 = 1

x 0

Lecture 35

18.01 Fall 2006

Lecture 35: Improper Integrals

Denition.
An improper integral, dened by
a

f (x)dx = lim
M a

f (x)dx

is said to converge if the limit exists (diverges if the limit does not exist). Example 1.
0

ekx dx = 1/k
0 M

(k > 0)

M ekx dx = (1/k )ekx = (1/k )(1 ekM )

0

Taking the limit as M , we nd ekM 0 and ekx dx = 1/k

0

We rewrite this calculation more informally as follows, k ekx dx = (1/k )ekx ) = 1/k = (1/k )(1 e
0 0

(since k > 0)

Note that the integral over the innite interval ekx dx = 1/k has an easier formula than the 0 M corresponding nite integral ekx dx = (1/k )(1 ekM ). As a practical matter, for large M , the term ekM is negligible, so even the simpler formula 1/k serves as a good approximation to the nite
integral. Innite integrals are often easier than nite ones, just as innitesimals and derivatives are
easier than dierence quotients.
Application: Replace x by t = time in seconds in Example 1.
R = rate of decay = number of atoms that decay per second at time 0.
At later times t > 0 the decay rate is Rekt (smaller by an exponential factor ekt )
Eventually (over time 0 t < ) every atom decays. So the total number of atoms N is calculated using the formula we found in Example 1, N= Rekt dt = R/k
0 0

The half life H of a radioactive element is the time H at which the decay rate is half what it was at the start. Thus ekH = 1/2 = kH = ln(1/2) 1 = k = (ln 2)/H

Lecture 35

18.01 Fall 2006

Hence R = N k = N (ln 2)/H Let us illustrate with Polonium 210, which has been in the news lately. The half life is 138 days or H = (138days)(24hr/day)(602 sec/hr) = (138)(24)(60)2 seconds Using this value of H , we nd that one gram of Polonium 210 emits (1 gram)(6 1023 /210 atoms/gram)(ln2)/H = 1.661014 decays/sec 4500 curies At 5.3 MeV per decay, Polonium gives o 140 watts of radioactive energy per gram (white hot). Polonium emits alpha rays, which are blocked by skin but when ingested are 20 times more dangerous than gamma and X-rays. The lethal dose, when ingested, is about 107 grams. Example 2. We calculate,
0

dx/(1 + x2 ) = /2.
0 M

M dx 1 = tan x = tan1 M /2 1 + x2
0

as M . (If = tan

M then /2 as M . See Figures 1 and 2.)

y = tan(x)
x = -/2

.
x

x = /2

Figure 1: Graph of the tangent function, M = tan .

Lecture 35

18.01 Fall 2006

y = arctan(x)

y= /2

Figure 2: Graph of the arctangent function, = tan1 M . Example 3.

0

.
M x = tan(y)
y = -/2

ex dx =

/2

Recall that we already computed this improper integral (by computing a volume in two ways, slices and the method of shells). This shows vividly that a nite integral can be harder to understand than its innite counterpart: M 2 ex dx
0

can only evaluated numerically. It has no elementary formula. By contrast, we found an explicit formula when M = . Example 4.
1

dx/x
1 M

M dx/x = ln x = ln M ln 1 = ln M
1

as M . This improper integral is innite (called divergent or not convergent). Example 5.

1

dx/xp
M

(p > 1)

M dx/xp = (1/(1 p))x1p = (1/(1 p))(M 1p 1) 1/(p 1)

1

as M because 1 p < 0. Thus, this integral is convergent. Example 6.

1

dx/xp

(0 < p < 1)

This is very similar to the previous example, but diverges M M p 1p dx/x = (1/(1 p))x = (1/(1 p))(M 1p 1) 1
1

as M because 1 p > 0. 3

Lecture 35

18.01 Fall 2006

Determining Divergence and Convergence

To decide whether an integral converges or diverges, dont need to evaluate. Instead one can compare it to a simpler integral that can be evaluated. dx The General Story for powers: xp 1 From Examples 4, 5 and 6 we know that this diverges (is innite) for 0 < p 1 and converges (is nite) for p > 1. The comparison of integrals says that a larger function has a larger integral. If we restrict ourselves to nonnegative functions, then even when the region is unbounded, as in the case of an improper integral, the area under the graph of the larger function is more than the area under the graph of the smaller one. Consider 0 f (x) g (x) (as in Figure 3)
y

f(x)

g(x) x x =a

Figure 3: The area under f (x) is less than the area under g(x) for a x < . If
a

g (x) dx converges, then so does

a

f (x) dx. (In other words, if the area under g is nite,

then the area under f , being smaller, must also be nite.) If f (x) dx diverges, then so does g (x) dx. (In other words, if the area under f is innite,
a a

then the area under g , being larger, must also be innite.) The way comparison is used is by replacing functions by simpler ones whose integrals we can calculate. You will have to decide whether you want to trap the function from above or below. This will depend on whether you are demonstrating that the integral is nite or innite.

Lecture 35

18.01 Fall 2006

Example 7.
0

dx It is natural to try the comparison x3 + 1 1 1 3/2 3 x x +1

But the area under x3/2 on the interval 0 < x < , dx 3/2 x 0 turns out to be innite because of the innite behavior as x 0. We can rescue this comparison by excluding an interval near 0. 1 dx dx dx = + 3 3 x +1 x +1 x3 + 1 0 0 1 The integral on 0 < x < 1 is a nite integral and the second integral now works well with comparison, dx dx < 3 x3/2 x +1 1 1 because 3/2 > 1. Example 8.
0

ex dx
1 x3

For x 1, x3 x, so

dx
1

ex dx = 1 <
3

Thus the full integral from 0 x < of ex converges as well. We can ignore the interval 3 0 x 1 because it has nite length and ex does not tend to innity there. Limit comparison: Suppose that 0 f (x) and lim f (x)/g (x) 1. Then f (x) 2g (x) for x a (some large a).
x
Hence f (x) dx 2 g (x) dx.
a a

(x + 10) dx x2 + 1 0 The limiting behavior as x is Example 9. 1 (x + 10)dx x 2 = x2 + 1 x x Since

1

dx = , the integral x

(x + 10) dx also diverges. x2 + 1

Lecture 35

18.01 Fall 2006

Example 10 (from PS8).

0

xn ex dx

This converges. To carry out a convenient comparison requires some experience with growth rates of functions. xn << ex not enough. Instead use xn /ex/2 0 (true by LHop). It follows that xn << ex/2 = xn ex << ex/2 ex = ex/2 Now by limit comparison, since ex/2 dx converges, so does our integral. You will deal with this integral on the problem set.
0

Improper Integrals of the Second Type

0 1

dx x

1 We know that as x 0. x
0 1

dx = lim x a0+

x1/2 dx

1 x1/2 dx = 2x1/2 = 2 2a1/2

a

As a 0, 2a1/2 0. So, Similarly,

0 1 1

x1/2 dx = 2

xp dx =

1 p + 1

for all p < 1. 1 For p = , 2

1 =2 1 +1 2

However, for p 1, the integral diverges.

Lecture 36

18.01 Fall 2006

Lecture 36: Innite Series and Convergence Tests

Innite Series
Geometric Series
A geometric series looks like 1 + a + a2 + a3 + ... = S Theres a trick to evaluate this: multiply both sides by a: a + a2 + a3 + ... = aS Subtracting, (1 + a + a2 + a3 + ) (a + a2 + a3 + ) = S aS In other words, 1 = S aS = 1 = (1 a)S = S = This only works when |a| < 1, i.e. 1 < a < 1. a = 1 cant work: 1 + 1 + 1 + ... = a = 1 cant work, either: 1 1 + 1 1 + ... = 1 1 = 2 1 (1) 1 1a

Notation
Here is some notation thats useful for dealing with series or sums. An innite sum is written:
k=0

ak = a0 + a1 + a2 + ...

The nite sum Sn =

n k=0

ak = a0 + ... + an

is called the nth partial sum of the innite series.

Lecture 36

18.01 Fall 2006

Denition
k=0

ak = s

means the same thing as

n

lim Sn = s, where Sn =

n k=0

ak

We say the series converges to s, if the limit exists and is nite. The importance of convergence is illustrated here by the example of the geometric series. If a = 1, S = 1 + 1 + 1 + ... = . But S aS = 1 or does not make sense and is not usable! =1

Another type of series:

1 p n n=1

We can use integrals to decide if this type of series converges. First, turn the sum into an integral: 1 dx p n xp 1 n=1 If that improper integral evaluates to a nite number, the series converges. Note : This approach only tells us whether or not a series converges. It does not tell us what number the series converges to. That is a much harder problem. For example, it takes a lot of work to determine 1 2 = n2 6 n=1 Mathematicians have only recently been able to determine that
1 3 n n=1

converges to an irrational number! Harmonic Series 1 dx n x 1 n=1 We can evaluate the improper integral via Riemann sums. Well use the upper Riemann sum (see Figure 1) to get an upper bound on the value of the integral.

Lecture 36

18.01 Fall 2006

y=x 1

1 2 3

Figure 1: Upper Riemann Sum.

dx 1 1 1 + + ... + = sN 1 sN x 2 N 1
N

We know that

As N , ln N , so sN

dx = ln N x 1 as well. In other words,

1 n n=1

diverges. Actually, sN approaches rather slowly. Lets take the lower Riemann sum (see Figure 2).

y=x

1 2

3 4

Figure 2: Lower Riemann Sum. sN = 1 + Therefore, ln N < sN < 1 + ln N 3 N N 1 dx 1 1 1+ = 1 + ln N + ... + =1+ n x 2 N 1 n=2

Lecture 36

18.01 Fall 2006

Integral Comparison
Consider a positive, decreasing function f (x) > 0. (For example, f (x) = f (n) f (x)dx < f (1) 1
n=1

1 ) xp

So, either both of the terms converge, or they both diverge. This is what we mean when we say 1 dx p n xp 1 n=1 Therefore,
1 diverges for p 1 and converges for p > 1. p n n=1

Lots of fudge room: in comparison. 1

2 n + 10
n=1 diverges, because 1 1 1 2 1/2 = 2 n (n ) n + 10 Limit comparison : If f (x) g (x) as x , then f (n) and g (n) either both converge or both diverge. What, exactly, does f (x) g (x) mean? It means that lim f (x) =c g (x)

where 0 < c < . Lets check: does the following series converge? n 5 n 10 n=1 n n 1 5/2 = n3/2 = 3/2 5 n n n 10 Since 3 > 1, this series does converge. 2

Lecture 36

18.01 Fall 2006

Playing with blocks

At this point in the lecture, the professor brings out several long, identical building blocks. Do you think its possible to stack the blocks like this?

Top block is farther out than the bottom block.

Figure 3: Collective center of mass of upper blocks is always over the base block. In order for this to work, you want the collective center of mass of the upper blocks always to be over the base block. The professor successfully builds the stack. Is it possible to extend this stack clear across the room? The best strategy is to build from the top block down. Let C0 be the left end of the rst (top) block. Let C1 = the center of mass of the rst block (top block). Put the second block as far to the right as possible, namely, so that its left end is at C1 (Figure 4). Let C2 = the center of mass of the top two blocks. Strategy : put the left end of the next block underneath the center of mass of all the previous ones combined. (See Figure 5).

Lecture 36

18.01 Fall 2006

1/2

C0

C1

C2

Figure 4: Stack of 2 Blocks.

2

1 2 3
1 1/2 1/3

C0

C1

C2 C3

Figure 5: Stack of 3 Blocks. Left end of block 3 is C2 = center of mass of blocks 1 and 2.

C0 = 0
C1 = 1
C2 = 1 + 1
2

Cn+1 =

nCn + 1(Cn + 1) (n + 1)Cn + 1 1

= = Cn + n+1 n+1 n + 1
1 1
C3 = 1 + + 2 3
1 1 1
C4 = 1 + + + 2 3 4
1 1 1 1
C5 = 1 + + + + > 2
2 3 4 5
6

Lecture 36

18.01 Fall 2006

}
center of mass of the rst n blocks

n
n+1 block

Figure 6: Stack of n + 1 Blocks. So yes, you can extend this stack as far (horizontally) as you want provided that you have enough blocks. Another way of looking at this problem is to say
N 1 = SN n n=1

Recall the Riemann Sum estimation from the beginning of this lecture: ln N < SN < (ln N ) + 1 as N , SN . How high would this stack of blocks be if we extended it across the two lab tables here at the front of the lecture hall? The blocks are 30 cm by 3 cm (see Figure 7). One lab table is 6.5 blocks, or 13 units, long. Two tables are 26 units long. There will be 26 2 = 24 units of overhang in the stack.

3 cm 30 cm

Figure 7: Side view of one block. If ln N = 24, then N = e24 . Height = 3 cm e24 8 108 m That height is roughly twice the distance to the moon. If you want the stack to span this room ( 30 ft.), it would have to be 1026 meters high. Thats about the diameter of the observable universe. 7

Lecture 37

18.01 Fall 2006

Lecture 37: Taylor Series

General Power Series

What is cos x anyway? Recall: geometric series 1 + a + a2 + = General power series is an innite sum: f (x) = a0 + a1 x + a2 x2 + a3 x3 + represents f when |x| < R where R = radius of convergence. This means that for |x| < R, |an xn | 0 as n (geometrically). On the other hand, if |x| > R, then |an xn | does not tend to 0. For 1 1 example, in the case of the geometric series, if |a| = , then |an | = n . Since the higher-order terms 2 2 get increasingly small if |a| < 1, the tail of the series is negligible. Example 1. If a = 1, |an | = 1 does not tend to 0. 1 1 + 1 1 + The sum bounces back and forth between 0 and 1. Therefore it does not approach 0. Outside the interval 1 < a < 1, the series diverges. 1 1a for | a |< 1

Basic Tools
Rules of polynomials apply to series within the radius of convergence. Substitution/Algebra 1 = 1 + x + x2 + 1x

Example 2. x = -u. 1 = 1 u + u2 u3 + 1+u

Example 3. x = v 2 . 1 = 1 v2 + v4 v6 + 1 + v2

Lecture 37

18.01 Fall 2006

Example 4. 1 1x 1 1x = (1 + x + x2 + )(1 + x + x2 + )

Term-by-term multiplication gives: 1 + 2x + 3x2 + 1 Remember, here x is some number like . As you take higher and higher powers of x, the result 2 gets smaller and smaller. Dierentiation (term by term) d 1 d = 1 + x + x2 + x3 + dx 1 x dx 1 = 0 + 1 + 2x + 3x2 + (1 x)2 where 1 is a0 , 2 is a1 and 3 is a2

Same answer as Example 4, but using a new method. Integration (term by term) where f (x) = a0 + a1 x + a2 x2 + a1 a2 f (x) dx = c + a0 + x2 + x3 + 2 3

Example 5.

du 1+u 1 = 1 u + u2 u3 + 1+u

du u2 u3 u4 =c+u + + 1+u 2 3 4 x du x2 x3 x4 ln(1 + x) = =x + + 3 4 2 0 1+u So now we know the series expansion of ln(1 + x).

Example 6. Integrate Example 3. 1 = 1 v2 + v4 v6 + 1 + v2 dv v3 v5 v7 =c+ v + + 1 + v2 3 5 7 x dv x3 x5 x7 tan1 x = =x + + 2 5 7 3 0 1+v

Lecture 37

18.01 Fall 2006

Taylors Series and Taylors Formula

If f (x) = a0 + a1 x + a2 x2 + , we want to gure out what all these coecients are. Dierentiating, f (x) = a1 + 2a2 x + 3a3 x2 + f (x) = (2)(1)a2 + (3)(2)a3 x + (4)(3)a4 x2 + f (x) = (3)(2)(1)a3 + (4)(3)(2)a4 x + Lets plug in x = 0 to all of these equations. f (0) = a0 ; f (0) = a1 ; f (0) = 2a2 ; f (0) = (3!)a3 Taylors Formula tells us what the coecients are: f (n) (0) = (n!)an Remember, n! = n(n 1)(n 2) (2)(1) and 0! = 1. Coecients an are given by: an = 1 n! f (n) (0)

Example 7. f (x) = ex . f (x) = ex f (x) = ex f (n) (x) = ex f (n) (0) = e0 = 1 Therefore, by Taylors Formula an = ex = Or in compact form, ex = Now, we can calculate e to any accuracy: e=1+1+ 1 1 1 1 + + + + 2 3! 4! 5!
xn n! n=0

1 and n! 1 1 1 1 + x + x2 + x3 + 0! 1! 2! 3!

Example 7. f (x) = cos x. f (x) = sin x f (x) = cos x 3

Lecture 37

18.01 Fall 2006

f (x) = sin x
f (4) (x) = cos x
f (0) = cos(0) = 1
f (0) = sin(0) = 0
f (0) = cos(0) = 1
f (0) = sin(0) = 0
Only even coecients are non-zero, and their signs alternate. Therefore, 1 1 1 1 cos x = 1 x2 + x4 x6 + x8 + 2 4! 8! 6! Note: cos(x) is an even function. So is this power series as it contains only even powers of x. There are two ways of nding the Taylor Series for sin x. Take derivative of cos x, or use Taylors formula. We will take the derivative: d 1 4 6 8 sin x = cos x = 0 2 x + x3 x5 + x7 + dx 2 4! 6! 8! = x + sin(x) = x x3 x5 x7 + + 3! 5! 7!

x3 x5 x7 + + 3! 5! 7!

Compare with quadratic approximation from earlier in the term: 1 cos x 1 x2 2 We can also write: cos x =
x2k x0 x2 1 (1)k = (1)0 + (1)2 + = 1 x2 + (2k )! 0! 2! 2 x2k+1 (1)k n = 2k + 1 (2k + 1)!

sin x x

k=0

sin x =

k=0

Example 8: Binomial Expansion. f (x) = (1 + x)a a a(a 1) 2 a(a 1)(a 2) 3 (1 + x)a = 1 + x + x + x + 1 2! 3!

Lecture 37

18.01 Fall 2006

Taylor Series with Another Base Point

A Taylor series with its base point at a (instead of at 0) looks like: f (3) (b) f (b) (x b)2 + (x b)3 + ... 2 3! Taylor series for x. Its a bad idea to expand using b = 0 because x is not dierentiable at x = 0.
Instead use b = 1.
1 1
1
1 2 2 x1/2 = 1 + (x 1) + (x 1)2 + 2 2! f (x) = f (b) + f (b)(x b) +

Lecture 38

18.01 Fall 2006

Lecture 38: Final Review

Review: Dierentiating and Integrating Series.
If f (x) =
n=0 n=1

an xn , then and f (x)dx = C +

an xn+1 n+1 n=0

f (x) =

nan x

n1

Example 1: Normal (or Gaussian) Distribution. x x 2 (t2 )2 (t2 )3 et dt = 1 t2 + + + dt 2! 3! 0 0 x t4 t6 t8 2 = 1 t + + ... dt 2! 3! 4! 0 =x Even though
0 x
2

x3 1 x5 1 x7 + + ... 3 2! 5 3! 7

et dt isnt an elementary function, we can still compute it. Elementary functions

are still a little bit better, though. For example: sin x = x x3 x5 (/2)3 (/2)5 + = sin = + 3! 5! 2 2 3! 5!

But to compute sin(/2) numerically is a waste of time. We know that the sum if something very simple, namely, sin = 1 2 Its not obvious from the series expansion that sin x deals with angles. Series are sometimes com plicated and unintuitive. Nevertheless, we can read this formula backwards to nd a formula for . Start with sin = 1. 2 2 Then, 1 1 dx = sin1 x = sin1 1 sin1 0 = 0 = 2 2 2 0 1x 0 We want to nd the series expansion for (1 x2 )1/2 , but lets tackle a simpler case rst: 1 1 1 1 1 1 2 1 1 2 2 2 2 2 (1 + u)1/2 = 1 + u+ u2 + u3 + 12 123 2 1 1 3 2 1 3 5
3 =1 u+ u u +
24 2 246 Notice the pattern: odd numbers go on the top, even numbers go on the bottom, and the signs alternate. 1

Lecture 38

18.01 Fall 2006

Now, let u = x2 . 1 13 4 135 6 x + x + (1 x2 )1/2 = 1 + x2 + 2 24 246 1 x3 1 3 x5 1 3 5 x7 2 1/2 (1 x ) dx = C + x + + + + 2 3 24 5 246 7 1 1 1 13 1 135 1 = (1 x2 )1/2 dx = 1 + + + + 2 2 3 24 5 246 7 0 Heres a hard (optional) extra credit problem: why does this series converge? LHpitals rule to nd out how quickly the terms decrease. Hint: use

The Final Exam

Heres another attempt to clarify the concept of weighted averages.

Weighted Average
A weighted average of some function, f , is dened as: b Average(f ) = Here,
a b a

w(x)f (x) dx b w(x) dx a

w(x) dx is the total, and w(x) is the weighting function.

Example: taken from a past problem set. You get $t if a certain particle decays in t seconds. How much should you pay to play? You were given that the likelihood that the particle has not decayed (the weighting function) is: w(x) = ekt Remember,
0

ekt dt =

1 k

The payo is f (t) = t The expected (or average) payo is kt f (t)w(t) dt te dt 0 = 0 kt w(t) dt e dt 0 0 =k
0

tekt dt =

(kt)ekt dt

Do the change of variable: u = kt and du = k dt

Lecture 38

18.01 Fall 2006

Average =
0

ueu

du k

On a previous problem set, you evaluated this using integration by parts:

0

ueu du = 1.

Average =
0

ueu

du 1 = k k

On the problem set, we calculated the half-life (H ) for Polonium120 was (131)(24)(60)2 seconds. We also found that ln 2 k= H Therefore, the expected payo is 1 H = k ln 2 where H is the half-life of the particle in seconds. Now, youre all probably wondering: who on earth bets on particle decays? In truth, no one does. There is, however, a very similar problem that is useful in the real world. There is something called an annuity, which is basically a retirement pension. You can buy an annuity, and then get paid a certain amount every month once you retire. Once you die, the annuity payments stop. You (and the people paying you) naturally care about how much money you can expect to get over the course of your retirement. In this case, f (t) = t represents how much money you end up with, and w(t) = ekt represents how likely your are to be alive after t years. What if you want a 2-life annuity? Then, you need multiple integrals, which you will learn about in multivariable calculus (18.02). Our rst goal in this class was to be able to dierentiate anything. In multivariable calculus, you will learn about another chain rule. That chain rule will unify the (single-variable) chain rule, the product rule, the quotient rule, and implicit dierentiation. You might say the multivariable chain rule is One One One And thing to rule them all thing to nd them thing to bring them all in a matrix bind them.

(with apologies to JRR Tolkien).