Introduction to Integral Calculus

Table of contents

Note: it is highly recommended to navigate by clicking links in the table of contents! It means you can use the back button in your browser to go back to any section you were reading, so you can jump back and forth between sections!

These are notes taken during RPI's MATH 1010 course, on the topic of integrals and their applications in calculus.

Antiderivatives

To introduce the topic of integral calculus, we must begin with understanding the notion of an antiderivative. An antiderivative (also called an indefinite integral, and we'll see why later) is the opposite of a derivative. The antiderivative of a particular function f(x)f(x) is the function whose derivative is f(x)f(x). For example, x2x^2 is the antiderivative of 2x2x, because 2x2x is the derivative of x2x^2.

Antiderivatives are represented by the "curly S" integral sign (again for reasons that will make sense later), with a dxdx representing the variable to take the antiderivative with respect to - in this case xx:

f(x)dx=F(x)+C \int f(x) dx = F(x) + C

Since the derivative of a constant is zero, the antiderivative of f(x)+1f(x) + 1, f(x)+3f(x) + 3, or f(x)+100f(x) + 100 is the same. Thus we add the CC known as the constant of integration to remind us that multiple antiderivatives differing by a constant have the same derivative.

Rules for antiderivatives (indefinite integral formulas)

The rules of antiderivatives (indefinite integrals) are the inverse of the rules for derivatives:

k dx=kx+Cxn dx=xn+1n+1+C1x dx=lnx+Cex dx=ex+Cbx=bxlnb+Ccos(x) dx=sin(x)+Csin(x) dx=cos(x)+Cdx1x2=arcsin(x)+Cdx1+x2=arctan(x)+Cdxxx21=arcsec(x)+C \begin{align*} \int k~dx &= kx + C \\ \int x^n~dx &= \frac{x^{n + 1}}{n + 1} + C \\ \int \frac{1}{x}~dx &= \ln |x| + C \\ \int e^x~dx &= e^x + C \\ \int b^x &= \frac{b^x}{\ln b} + C \\ \int \cos(x)~dx &= \sin(x) + C \\ \int \sin(x)~dx &= -\cos(x) + C \\ \int \frac{dx}{\sqrt{1 - x^2}} &= \operatorname{arcsin}(x) + C \\ \int \frac{dx}{1 + x^2} &= \operatorname{arctan}(x) + C \\ \int \frac{dx}{|x|\sqrt{x^2 - 1}} &= \operatorname{arcsec}(x) + C \end{align*}

Note that the last few formulas for inverse trigonometric functions can be generalized to the following:

dua2u2=arcsin(ua)+Cdua2+u2=1aarctan(ua)+Cduuu2a2=1aarcsec(ua)+C \begin{align*} \int \frac{du}{\sqrt{a^2 - u^2}} &= \operatorname{arcsin} \left(\frac{u}{a}\right) + C \\ \int \frac{du}{a^2 + u^2} &= \frac{1}{a} \operatorname{arctan} \left(\frac{u}{a}\right) + C \\ \int \frac{du}{|u|\sqrt{u^2 - a^2}} &= \frac{1}{a} \operatorname{arcsec} \left(\frac{u}{a}\right) + C \end{align*}

Where u=u(x)u = u(x) is some function of xx, and du=u(x)dxdu = u'(x) dx. For more formulas of integrals, a good resource to reference is https://www.integral-table.com/.

Initial value problems

To find the value of the CC when doing derivatives, we can use an initial value problem. For instance, say:

f(0)=1, f(x)=2x dx f(0) = 1,~f(x) = \int 2x~dx

Thus, f(x)=x2+Cf(x) = x^2 + C. Then, plugging in the initial value, f(0)=1f(0) = 1, so (02)+C=1(0^2) + C = 1, so C=1C = 1. Therefore:

f(x)=x2+1 f(x) = x^2 + 1

Sums

We will now jump to a topic that might initially appear to be completely unrelated - sums and areas. A sum is just what it means - adding a group of numbers (or variables) together. Sums are denoted by the summation sign \displaystyle \sum. For instance, one sum could be:

i=0nai \sum_{i = 0}^n a_i

This means to sum from a0a_0 to ana_n - that is:

i=0nai=a0+a1+a2+a3++an1+an \sum_{i = 0}^n a_i = a_0 + a_1 + a_2 + a_3 + \dots + a_{n - 1} + a_n

Rules of sums

If we are given a sum, we might be interested in what the sum evaluates to (or we might not, but pretend that we are for the time being). Evaluating sums is actually not possible in many cases. However, in cases where a sum can be evaluated, it often helps to use the below identities:

i=1ncai=ci=1naii=1nc=cni=1nai±bi=i=1nai±i=1nbi \begin{align*} \sum_{i = 1}^n c a_i &= c \sum_{i = 1}^n a_i \\ \sum_{i = 1}^n c &= cn \\ \sum_{i = 1}^n a_i \pm b_i &= \sum_{i = 1}^n a_i \pm \sum_{i = 1}^n b_i \end{align*}

There also exist special sums that have well-known solutions - these are extremely helpful for evaluating more complicated sums, and they are given below:

i=1ni=n(n+1)2i=1ni2=n(n+1)(2n+1)6i=1ni3=(n(n+1)2)2 \begin{align*} \sum_{i = 1}^n i &= \frac{n(n + 1)}{2} \\ \sum_{i = 1}^n i^2 &= \frac{n(n + 1)(2n + 1)}{6} \\ \sum_{i = 1}^n i^3 &= \left(\frac{n(n + 1)}{2}\right)^2 \end{align*}

Finding and estimating areas

We will now turn to a very useful application of sums - finding areas. The idea is this: if we had some shape modelled by a function f(x)f(x), then we can approximate the area under the curve by "slicing" the area under the curve into small rectanges, each of width Δx\Delta x (which is some small number) and of height f(xi)f(x_i). This is called the left-hand approximation of an area under the curve, and is given by:

AL=i=1nf(xi)Δx A_L = \sum_{i = 1}^{n} f(x_i) \Delta x

Given that:

Δx=ban,xi=a+(i1)Δx.xi+1=a+iΔx \begin{matrix*} \Delta x = \dfrac{b - a}{n}, &x_i = a + (i - 1)\Delta x. &x_{i + 1} = a + i \Delta x \end{matrix*}

The right-hand approximation of an area under the curve is given by:

AR=i=1nf(xi+1)Δx A_R = \sum_{i = 1}^n f(x_{i + 1}) \Delta x

The midpoint approximation to the area is given by:

AM=i=1nf(xi+xi+1n)Δx A_M = \sum_{i = 1}^{n} f\left(\frac{x_i + x_{i + 1}}{n}\right) \Delta x

Now, these are approximations for finite sums, but in the limit as we sum over increasingly many and increasingly tiny rectangles, we obtain a formula for the exact area under a curve:

A=limni=1nf(xi)Δx A = \lim_{n \to \infty} \sum_{i = 1}^n f(x_i) \Delta x

For example, consider the area under the curve f(x)=4x2f(x) = 4 - x^2 from -2 to 2. We would then have:

A=limni=1n(4(2+iΔx)2)Δx A = \lim_{n \to \infty} \sum_{i = 1}^n (4 - (-2 + i \Delta x)^2) \Delta x

Since Δx=ban\Delta x = \frac{b - a}{n}, where b=2b = 2 and a=2a = -2, we would have:

A=limni=1n(4(2+4in)2)4n A = \lim_{n \to \infty} \sum_{i = 1}^n \left(4 - \left(-2 + \frac{4i}{n}\right)^2\right) \frac{4}{n}

This simplifies to:

A=limni=1n64in264i2n3 A = \lim_{n \to \infty} \sum_{i = 1}^n \frac{64i}{n^2} - \frac{64i^2}{n^3}

Using the sum rules, we get:

A=limni=1n64n2n(n+1)264n3n(n+1)(2n+1)6 A = \lim_{n \to \infty} \sum_{i = 1}^n \frac{64}{n^2} \frac{n(n + 1)}{2} - \frac{64}{n^3} \frac{n(n + 1)(2n + 1)}{6}

Which simplifies to:

A=limni=1n64(n+1)2n64(n+1)(2n+1)6n2 A = \lim_{n \to \infty} \sum_{i = 1}^n \frac{64(n + 1)}{2n} - \frac{64 (n + 1)(2n + 1)}{6n^2}

Expanding this and simplifying, we get:

A=limni=1n323323n2 A = \lim_{n \to \infty} \sum_{i = 1}^n \frac{32}{3} - \frac{32}{3n^2}

Now, taking the limit as nn \to \infty, we get:

A=323 A = \frac{32}{3}

Riemann sums and the definite integral

From our previous discussion of the limit of sums, we can now formalize the idea of computing the area under a curve. A Riemann sum is given by:

k=1Nf(xk)Δxk \sum_{k = 1}^N f(x_k) \Delta x_k

The definite integral is given by the limit of the Riemann sum, and represents the area under a curve, where AA is the area and y=f(x)y = f(x) is the equation of the curve:

A=abf(x)dx=limNk=1Nf(xk)Δxk A = \int_a^b f(x)dx = \lim_{N \to \infty} \sum_{k = 1}^N f(x_k) \Delta x_k

The definite integral also has some convenient properties that will come in handy later:

abf(x)dx=baf(x)dxaaf(x)dx=0abcf(x)dx=cabf(x)dxacf(x)dx=abf(x)dx+bcf(x)dx \begin{gather*} \int_a^b f(x) dx = -\int_b^a f(x) dx \\ \int_a^a f(x) dx = 0 \\ \int_a^b c f(x) dx = c\int_a^b f(x) dx \\ \int_a^c f(x) dx = \int_a^b f(x) dx + \int_b^c f(x) dx \end{gather*}

Evaluating the definite integral is called integration. The definite integral can be evaluated in simple cases from just geometry, especially when the curve forms a semicircle, triangle, rectangle, or trapezoid with the x-axis. However, in most cases, integration is not an easy task - some definite integrals are almost impossible to solve!

Applications of integration

We described integration from the point of view of finding the area under a curve, but integration is far more general than an area-finding method. Indeed, integration is ubiquitous in the mathematical sciences, with myriad applications:

The Fundamental Theorem of Calculus

There are two parts to the Fundamental Theorem of Calculus (FTC). Part 1 of the FTC states that:

F(x)=f(x)dxabf(x)dx=F(b)F(a) F(x) = \int f(x) dx \Rightarrow \int_a^b f(x) dx = F(b) - F(a)

That is, if F(x)F(x) is the indefinite integral of f(x)f(x), then the definite integral of f(x)f(x) can be found by evaluating F(x)F(x) at 2 points.

Part 2 of the FTC states that:

ddxaxf(t)dt=f(x) \frac{d}{dx} \int_a^x f(t) dt = f(x)

That is, the derivative of the indefinite integral of a function from any number to xx is equal to the same function evaluated at xx. Or more simply, integrals are the inverse of derivatives.

Integration by substitution

Consider a composite function of the form f(g(x))g(x)f(g(x)) g'(x). If we make the substitution u=g(x)u = g(x), then we can rewrite this as f(u)uf(u) u'. Then the chain rule applies in reverse:

f(u(x))u(x)dx=f(u)du=F(u)+C \int f(u(x)) u'(x) dx = \int f(u) du = F(u) + C

This is called integration by u-substitution. The trick is to find an "inner function" uu in a composite function whose derivative also appears in that same composite function.

As an example, we want to evaluate:

2t4dt \int \frac{2}{t- 4} dt

We want to find an inner function uu whose derivative also appears. That is, we want to rewrite the function in the form:

f(u)udt \int f(u) u' dt

Here, we can choose u=t4u = t - 4, which means f(u)=2uf(u) = \frac{2}{u}. So this means we have:

f(u)udt=2uudt \int f(u) u' dt = \int \frac{2}{u} u' dt

But we still need to find uu'. Taking the derivative of uu, we have:

dudt=1 \frac{du}{dt} = 1

Then we can apply the substitution rule:

2uudt=2u(1)du=2lnu+C=2lnt4+C \int \frac{2}{u} u' dt = \int \frac{2}{u} (1) du = 2\ln |u| + C = 2\ln | t- 4| + C

For a more complicated example, we may perhaps want to integrate:

cos(lnx)xdx \int \frac{\cos (\ln x)}{x} dx

Here, we let u=lnxu = \ln x. We can employ a notational shorthand for expediting the substitution rule that does not require rewriting the function in a form f(u)uf(u) u'. Instead, we can say that:

u=lnxdudx=1x u = \ln x \Rightarrow \frac{du}{dx} = \frac{1}{x}

xdu=dx x\, du = dx

Then, we can substitute in x,dux, du for dxdx:

cos(lnx)xdx=cosuxxdu=cosu du=sin(u)+C=sin(lnx)+C \int \frac{\cos (\ln x)}{x} dx = \int \frac{\cos u}{x} x\, du = \int \cos u~du = \sin(u) + C = \sin(\ln x) + C

To evaluate definite integrals with substitution, we have:

abf(u)udx=u(a)u(b)f(u)du \int_a^b f(u) u' dx = \int_{u(a)}^{u(b)} f(u) du

Integration of functions using the natural log trick

A direct consequence of the substitution method for integrals is the "natural log trick":

uu=lnu+C \int \frac{u'}{u} = \ln | u | + C

As an example to illustrate why this formula is the case, let us take the example of evaluating the following integral:

tan(5x)dx \int \tan(5x) dx

We first rewrite tangent in terms of sine and cosine:

sinxcosxdx \int \frac{\sin x}{\cos x} dx

If we use the substitution u=cosxu = \cos x, we have::

dudx=sinxdx=dusinx \begin{gather*} \frac{du}{dx} = -\sin x dx = \frac{du}{-\sin x} \end{gather*}

This means that:

sinxcosxdx=sinxsinx1udu=1udu=lncos(u)+C \begin{align*} \int \frac{\sin x}{\cos x} dx &= \int \frac{\sin x}{-\sin x} \frac{1}{u} du \\ &= -\int \frac{1}{u} du \\ &= -\ln | \cos(u) | + C \end{align*}

Differential equations with the natural log trick

We may use the natural log trick to solve not only integrals, but also differential equations. Differential equations specify a derivative of an unknown function as a function of xx, yy, or some other variables; the goal is to find the exact form of the unknown function. Consider, for instance, the differential equation:

dydx=2xx29 \frac{dy}{dx} = \frac{2x}{x^2 - 9}

To solve this differential equation, we can multiply by dxdx on both sides. However, this must be done in secret because this operation is technically "illegal" in mathematical terms (derivatives are not fractions; it just happens that sometimes they can be treated as such). One must then leave the scene of the mathematical crime and erase all the evidence used to find the answer. So we shall now commit the crime in question, and multiply by dxdx on both sides, which results in:

dy=2xx29dx dy = \frac{2x}{x^2 - 9} dx

If we integrate both sides, we then have:

dy=y=2xx29dx \int dy = y = \int \frac{2x}{x^2 - 9} dx

This is where we can use the natural log trick. Let u=x29u = x^2 - 9. Then:

dudx=2xdu2x=dxy=2x2x1udu=1udu=lnx29+C \begin{align*} \frac{du}{dx} &= 2x \\ \frac{du}{2x} &= dx \\ y &= \int \frac{2x}{2x} \frac{1}{u} du \\ &= \int \frac{1}{u} du \\ &= \ln|x^2 - 9| + C \end{align*}

Which is the general solution to our differential equation.

Average value of a function

We may also use integration to find the average value of a function, which is very useful for continuously-varying data (like velocity or position). The average value of a function is given by:

1baf(x)dx \frac{1}{b-a} \int f(x) dx

The mean value theorem for integrals states that if ff is continuous on [a,b][a, b], there is at least one cc value that satisfies:

f(c)=1baf(x)dx f(c) = \frac{1}{b - a} \int f(x) dx

Remember, as with the MVT for derivatives, to ignore any values not in the interval!

Area between curves

A region is vertically simple if it is bounded on the left and right by vertical lines, and bounded on the top and bottom by functions of xx. A region is horizontally simple if it is bounded on the top and bottom by horizontal lines, and bounded on the left and right by functions of yy.

The area between curves f(x)f(x) and g(x)g(x), assuming f(x)g(x)f(x) \geq g(x) for all x[a,b]x \in [a, b], is given by:

abf(x)g(x)dx \int_a^b f(x) - g(x) dx

Other specific cases: For curves f(y)f(y) and g(y)g(y), just replace every xx with yy. If f(x)g(x)f(x) \ngeq g(x) for all x[a,b]x \in [a, b], the area must be split into regions that are separately evaluated to make the condition true.

Solids of revolution

Another application of integrals generalizes the area-finding process to finding the surface area and volume of a particular shape. Such shapes must have a specific requirement: they have to be rotationally-symmetric, and thus we call them solids of revolution. By nature of their rotational symmetry, they are obtained by the rotation of a region in the xy-plane about an axis.

Note: Unlike many other uses of integration, this particular application of integrals is not very widely used in physics or applied math areas, but they are a good teaching/learning tool to practice integration.

Disk method

Using the disk method of finding the volume of a solid of revolution, the solid of revolution created from the region underneath the curve f(x)f(x) has the volume:

V=abπ[f(x)]2dx V = \int_a^b \pi [f(x)]^2 dx

And for a function f(y)f(y), the volume is given by:

V=abπ[f(y)]2dy V = \int_a^b \pi [f(y)]^2 dy

For example, say we wanted to find the volume of a sphere. A sphere is simply the solid of revolution created by rotating the curve f(x)=r2x2f(x) = \sqrt{r^2 - x^2} about the x-axis. So, we would have:

V=πrr(r2x2)2 dx V = \pi \int_{-r}^r (\sqrt{r^2 - x^2})^2~dx

Which becomes:

V=πrrr2x2dx V = \pi \int_{-r}^r r^2 - x^2 dx

Since rr is a constant, the integral becomes:

V=43πr3 V = \frac{4}{3} \pi r^3

Which is the volume of a sphere!

Washer method

Using the washer method of finding the volume of a solid of revolution, the solid of revolution created by the region between two curves f(x)f(x) and g(x)g(x) is given by:

V=abπ[f(x)]2π[g(x)]2dx V = \int_a^b \pi [f(x)]^2 - \pi [g(x)]^2 dx

And similarly for yy:

V=abπ[f(y)]2π[g(y)]2dy V = \int_a^b \pi [f(y)]^2 - \pi [g(y)]^2 dy

Shells method

Using the shells method of finding the volume of a solid of revolution, the solid of revolution created by the region underneath the curve f(x)f(x) is given by:

V=ab2πrf(x)dx V = \int_a^b 2\pi rf(x) dx

And similarly for yy:

V=ab2πrf(y)dy V = \int_a^b 2\pi rf(y) dy

If the region is instead between two curves f(x)f(x) and g(x)g(x), then the volume is instead given by:

V=ab2πr[f(x)g(x)]dx V = \int_a^b 2\pi r[f(x) - g(x)] dx

Or for yy:

V=ab2πr[f(y)g(y)]dy V = \int_a^b 2\pi r[f(y) - g(y)] dy

For instance, suppose we were to find the volume of the solid of revolution created by revolving the region bound by f(x)=x3,x[0,2]f(x) = x^3, x \in [0, 2] and the x-axis about the line x=3x = 3. A visualization of this region can be found at https://www.desmos.com/calculator/ymuwqk1xk7.

To compute the area, we'll first find the radii of the shells, and set that. In this case, the radius would be 3x3 - x, as that is the distance from the rotational axis to the edge of the shell. The bounds of integration would be from 0 to 2, as we want to create shells that fill up the space of the region, and the region is bounded by x=0x = 0 and x=2x = 2. Therefore, the integral is:

2π02(3x)x3dx=56π5 2\pi \int_0^2 (3 - x) x^3 dx = \frac{56\pi}{5}

Integration by Parts

Integration by parts is analagous to the product rule for derivatives. It is given by:

udv=uvvdu \int udv = uv - \int vdu

where:

v=dv v = \int dv

Typically, we want uu to be a function that has a simple derivative (such as a polynomial or lnx\ln x), and vv to be a function that has a simple integral (such as cosx\cos x or sec2x\sec^2 x or exe^x). It's typically much harder to integrate than differentiate, so pick dvdv first, and uu second.

In addition, integration by parts can be done multiple times over and over to repeatedly expand out an integral until it can be solved.

As an example of integration by parts, consider:

xcos(2x)dx \int x \cos (2x) dx

Here, we want our uu to be xx (because that's easy to differentiate), and dv=cos2xdv = \cos 2x. This means that du=dxdu = dx. To find vv, we simply do:

v=dv=cos(2x)=12sin(2x) v = \int dv = \int \cos (2x) = \frac{1}{2} \sin(2x)

After substitution, we obtain:

udv=12xsin(2x)12sin(2x)dx \int udv = \frac{1}{2}x \sin(2x) - \int \frac{1}{2} \sin(2x) dx

Which evaluates to:

12xsin(2x)+14cos(2x)+C \frac{1}{2} x \sin(2x) + \frac{1}{4} \cos(2x) + C

For a faster method of integration by parts, we can use the tabular method, which works so long as uu is a polynomial. For instance, consider:

x3exdx \int x^3 e^x dx

We construct a table of uu and all its derivatives (D) on the left, and dvdv and all its integrals (I) on the right, and alternating signs in front, until the derivative is zero:

SignDI
+x3x^3exe^x
-3x23x^2exe^x
+6x6xexe^x
-66exe^x
+00exe^x

Now, we start from the first "D" term, and match terms crosswise down diagonally like this:

SignDI
+1
-21
+32
-43
+54
.........

In front of each derivative term we put the sign associated with it. Therefore, our result is:

(+)x3ex()3x2ex(+)6xex()6ex+C (+) x^3 e^x (-) 3x^2 e^x (+)6x e^x (-) 6e^x + C

Or simplified:

x3ex3x2ex+6xex6ex+C x^3 e^x - 3x^2 e^x + 6xe^x - 6e^x + C

A third type of integration-by-parts integrals utilize a specialized technique. Consider:

excosxdx \int e^x \cos x dx

Doing a first integration by parts with u=exu = e^x and dv=cosxdv = \cos x, we get:

excosxdx=exsinxexsinxdx \int e^x \cos x dx = e^x \sin x - \int e^x \sin x dx

Now repeating the integration by parts on the second integral with u=exu = e^x and dv=sinxdv = \sin x, we get:

excosxdx=exsinx(excosx(excosxdx)) \int e^x \cos x dx = e^x \sin x - \left(-e^x \cos x - \int(-e^x\cos x dx)\right)

Or simplifying:

excosxdx=exsinx+excosxexcosxdx \int e^x \cos x dx = e^x \sin x + e^x \cos x - \int e^x \cos x dx

Notice how the integral appears twice in the expression. If we say that:

I=excosxdx I = \int e^x \cos x dx

We can rewrite the integral as:

I=exsinx+excosxI I = e^x \sin x + e^x \cos x- I

Solving for II, we get:

I=exsinx+excosx2+C I = \frac{e^x \sin x + e^x \cos x}{2} + C

which is the solution to the integral.

Trigonometric integrals

Trigonometric integrals are integrals that involve combination of trigonometric functions.

The first general category of trigonometric integrals are in the form:

sinmxcosnxdx \int \sin^m x \cos^n x dx

There are several common approaches for such integrals:

sin2x=1cos2x2cos2x=1+cos2x2sinxcosx=sin2x2 \begin{align*} \sin^2x &= \frac{1 - \cos 2x}{2} \\ \cos^2 x &= \frac{1 + \cos 2x}{2} \\ \sin x \cos x &= \frac{\sin 2x}{2} \end{align*}

Meanwhile, for trigonometric integrals in the general form:

tanmxsecnxdx \int \tan^m x \sec^n x dx

For instance, consider:

sin5xcos4xdx \int \sin^5 x \cos^4 x dx

We can break this into:

sin4xsinxcos4xdx \int \sin^4 x \sin x \cos^4 xdx

Which we can simplify into:

(sin2x)2sinxcos4xdx \int (\sin^2 x)^2 \sin x \cos^4 x dx

And again, into:

(1cos2x)2sinxcos4xdx \int (1 - \cos^2 x)^2 \sin x \cos^4 x dx

Now, we can use u=cosxu = \cos x, thus the integral becomes:

(1u2)2u4du -\int (1 - u^2)^2 u^4 du

Or expanded and with the negative sign distributed:

2u6u4u8du \int 2u^6 - u^4 - u^8 du

This becomes:

2u77u55u99+C \frac{2u^7}{7} - \frac{u^5}{5} - \frac{u^9}{9} + C

And substituting back, we have:

2cos7x7cos5x5cos9x9+C \frac{2 \cos^7 x}{7} - \frac{\cos^5 x}{5} - \frac{\cos^9 x}{9} + C

Trigonometric substitutions

For integrals in the form:

(a2±x2n)mdx \int (a^2 \pm x^{2n})^m dx

or generally, any integrals that have a2u2\sqrt{a^2 - u^2}, a2+u2\sqrt{a^2 + u^2}, or a2+u2a^2 + u^2, a trigonometric substitution is often helpful.

The substitutions make use of the 3 Pythagorean identities:

{sin2θ+cos2θ=11+tan2θ=sec2θ1+cot2θ=csc2θ \begin{cases} \sin^2 \theta + \cos^2 \theta = 1 \\ 1 + \tan^2 \theta = \sec^2 \theta \\ 1 + \cot^2 \theta = \csc^2 \theta \end{cases}

Performing a trigonometric substitution involves fitting an integrand into one of these forms. For instance, consider the integral:

4x2x2dx \int \frac{\sqrt{4 - x^2}}{x^2} dx

Here, we see the square root, which suggests the use of a trigonometric substitution. We first simplify the integral with:

4(1x24)x2dx \int \frac{\sqrt{4(1 - \frac{x^2}{4})}}{x^2} dx

We now have an expression in the form 1u21 - u^2. It makes sense to use the first Pythagorean identity. We say that:

x24=sin2θ \frac{x^2}{4} = \sin^2 \theta

Which means:

1x24=1sin2θ=cos2θ 1 - \frac{x^2}{4} = 1 - \sin^2 \theta = \cos^2 \theta

If we solve for xx, we get:

x=2sinθ x = 2\sin \theta

So:

dx=2cosθdθ dx = 2\cos \theta d\theta

Plugging in 1x24=cos2θ1 - \frac{x^2}{4} = \cos^2 \theta and dx=2cosθdθdx = 2\cos \theta d\theta into the integral, we get:

4cos2θ4sin2θ2cosθdθ \int \frac{\sqrt{4\cos^2 \theta}}{4\sin^2 \theta} 2 \cos \theta d\theta

Which simplifies to:

cot2θdθ \int \cot^2 \theta d\theta

We use the identity cot2θ+1=csc2θ\cot^2 \theta + 1 = \csc^2 \theta to get:

csc2θ1dθ \int \csc^2\theta - 1 d\theta

Whose solution is:

cotθθ+C -\cot \theta - \theta + C

We now need to convert the answer in terms of θ\theta back into an answer in terms of xx. To do so, remember that:

x=2sinθ x = 2 \sin \theta

Therefore:

sinθ=x2=OppositeHypotenuse \sin \theta = \frac{x}{2} = \frac{\text{Opposite}}{\text{Hypotenuse}}

We also know that:

cotθ=1tanθ=AdjacentOpposite \cot \theta = \frac{1}{\tan \theta} = \frac{\text{Adjacent}}{\text{Opposite}}

And by Pythagoras's theorem, we have:

Adjacent2=Hypotenuse2Opposite2=22x2=4x2 \begin{align*} \text{Adjacent}^2 &= \text{Hypotenuse}^2 - \text{Opposite}^2 \\ &= 2^2 - x^2 \\ &= 4 - x^2 \end{align*}

So:

Adjacent=4x2 \text{Adjacent} = \sqrt{4 - x^2}

And:

cotθ=4x2xθ=arcsin(x2) \begin{align*} \cot \theta &= \frac{\sqrt{4 - x^2}}{x} \\ \theta &= \operatorname{arcsin} \left(\frac{x}{2}\right) \end{align*}

We can now rewrite the solution in terms of xx:

4x2xarcsin(x2)+C -\frac{\sqrt{4 - x^2}}{x} - \operatorname{arcsin} \left(\frac{x}{2}\right) + C

Similarly, we can solve the following integral by trigonometric substitution:

dt(9t2+4)2 \int \frac{dt}{(9t^2 + 4)^2}

We first factor:

dt(4(94t2+1))2 \int \frac{dt}{(4(\frac{9}{4}t^2 + 1))^2}

We set:

1+94t2=sec2θ94t2=tan2θ \begin{align*} 1 + \frac{9}{4} t^2 &= \sec^2 \theta \\ \frac{9}{4} t^2 &= \tan^2 \theta \end{align*}

Therefore:

t=23tanθdt=23sec2θdθ \begin{align*} t &= \frac{2}{3} \tan \theta \\ dt &= \frac{2}{3} \sec^2 \theta d \theta \end{align*}

Subtituting, we get:

23dt(4sec2θ)2)sec2θdθ \frac{2}{3} \int \frac{dt}{(4\sec^2 \theta)^2)} \sec^2 \theta d \theta

Which becomes:

148cos2θdθ \frac{1}{48} \int \cos^2 \theta d \theta

Trigonometric substitutions are incredibly easy to mess up, so be careful with them, and make sure to avoid these errors:

Integration by long division and partial fractions

Consider an integral in the form:

P(x)Q(x)dx \int \frac{P(x)}{Q(x)} dx

If the degree of P(x)P(x) is higher than the degree of Q(x)Q(x), then the integrand can be simplified by dividing P(x)P(x) by Q(x)Q(x) using polynomial long division. The process is analogous to arithmetic long division.

When doing long division, be very careful of signs while subtracting. It is easy to get the wrong result just by forgetting a negative sign when doing a subtraction. Also, don't forget to carry over the remainder. Whenever possible, split up the resulting fraction when you integrate, so for instance, something like:

2xdx+5x+3x1dx \int 2x dx + \frac{5x + 3}{x - 1} dx

becomes:

2xdx+5xx1dx+3x1dx \int 2x dx + \int \frac{5x}{x - 1} dx + \frac{3}{x - 1} dx

In cases where long division cannot be done, a partial fraction decomposition can be performed instead. The goal of partial fraction decomposition is to write the integrand as a sum of simpler fractions. That is:

P(x)Q(x)=AF1(x)+BF2(x)+CF3(x)++ZFn(x) \frac{P(x)}{Q(x)} = \frac{A}{F_1(x)} + \frac{B}{F_2(x)} + \frac{C}{F_3(x)} + \dots + \frac{Z}{F_n(x)}

The first step of a partial fraction decomposition is to factor Q(x)Q(x) as completely as possible. That is, the integrand should take the following form:

P(x)Q(x)=P(x)F1(x)F2(x)F3(x)Fn(x) \frac{P(x)}{Q(x)} = \frac{P(x)}{F_1(x)F_2(x)F_3(x)\dots F_n(x)}

Where F1,F2,F3,,FnF_1, F_2, F_3, \dots, F_n are the factors of Q(x)Q(x). If a factor of Q(x)Q(x) is in the form (xa)n(x - a)^n, where aa is not a complex number or radical, it is called a linear factor. Such a factor contributes a partial fraction decomposition of the form:

1(xa)n=Axa+B(xa)2+C(xa)3++Z(xa)n \frac{1}{(x - a)^n} = \frac{A}{x - a} + \frac{B}{(x - a)^2} + \frac{C}{(x - a)^3} + \dots + \frac{Z}{(x - a)^n}

Note that a factor in the form xnx^n (such as x2x^2 or xx) can technically be written in the form (x0)n(x - 0)^n, and so also qualifies as a linear factor.

Meanwhile, if the factor is in the form (x2+ax+b)n(x^2 + ax + b)^n, and cannot be simplified further without using radicals or complex numbers, it is called an irreducible quadratic factor. Such a factor contributes a partial fraction decomposition of the form:

1(x2+ax+b)n=Ax+C1x2+ax+b+Bx+C2(x2+ax+b)2++Zx+Cn(x2+ax+b)n \frac{1}{(x^2 + ax + b)^n} = \frac{Ax + C_1}{x^2 + ax + b} + \frac{Bx + C_2}{(x^2 + ax + b)^2} + \dots + \frac{Zx + C_n}{(x^2 + ax + b)^n}

As an example, consider the integral:

x+5x2+x2dx \int \frac{x + 5}{x^2 + x - 2} dx

We first factor:

x+5(x1)(x+2)dx \int \frac{x + 5}{(x - 1)(x + 2)} dx

Now we perform a partial fraction composition of the integrand. Since we see that the factors are linear, we propose that the integrand take the form:

x+5(x1)(x+2)=Ax1+Bx+2 \frac{x + 5}{(x - 1)(x + 2)} = \frac{A}{x - 1} + \frac{B}{x + 2}

If we distribute by multiplying every term by (x1)(x+2)(x - 1)(x + 2), we get:

x+5=A(x+2)+B(x1) x + 5 = A(x + 2) + B(x - 1)

We can now substitute in convenient values to find AA and BB. For instance, if we evaluate with x=1x = 1, then:

(1+5)=A(1+2)+B(11) (1 + 5) = A(1 + 2) + B(1 - 1)

which simplifies to:

6=3A 6 = 3A

Therefore A=2A = 2. Similarly, if we evaluate with x=2x = -2, then:

(2+5)=A(2+2)+B(21) (-2 + 5) = A(-2 + 2) + B(-2 - 1)

which simplifies to:

3=3B 3 = -3B

Therefore A=2,B=1A = 2, B = -1. Placing the partial fraction decomposition back into the integrand, we have:

2x11x+2 dx \int \frac{2}{x - 1} - \frac{1}{x + 2}~dx

which becomes:

2lnx1lnx+2+C 2 \ln | x - 1 | - \ln | x + 2 | + C

When doing partial fraction decompositions, make sure that if the fraction is improper, that you first do long division, and then (if the result isn't simple enough to directly integrate) decompose using partial fractions. In addition, be careful about identifying factors - sometimes factors that might look linear factorable might be irreducible quadratics, and factors that look like quadratics are linear factorable.

Finally, when convenient cancelling isn't possible, use variations of x-values (e.g. evaluate when x=0x = 0, then x=1x = 1, then x=2x = 2) for getting partial fraction decompositions when there are multiple repeated factors. This allows you to setup a system of equations that you can solve to find the coefficients of the partial fraction decomposition.

After finding the partial fraction decomposition of a rational function, it is always best to plug in a value of xx (such as x=0x = 0 or x=1x = 1) and check that the partial fraction outputs the same value as the original integrand. This means the partial fraction decomposition is correct.

Improper integrals

A improper integral is an integral defined either on an unbounded domain, or a domain that is not integratable, then:

akf(x)dx=limbkabf(x)dx \int_a^k f(x) dx = \lim_{b \to k} \int_a^b f(x) dx

An improper integral satisfies the identities:

0f(x)dx=limb0bf(x)dx0f(x)dx=limbb0f(x)dxf(x)dx=limaa0f(x)dx+limb0bf(x)dx \begin{align*} \int_0^\infty f(x)dx &= \lim_{b \to \infty} \int_0^b f(x) dx \\ \int_{-\infty}^0 f(x)dx &= \lim_{b \to -\infty} \int_b^0 f(x) dx \\ \int_{-\infty}^\infty f(x) dx &= \lim_{a \to -\infty} \int_a^0 f(x) dx + \lim_{b \to \infty} \int_0^b f(x) dx \end{align*}

If the limit exists (evaluates to a finite number), then the integral is said to converge; otherwise, the integral is said to diverge.

For instance, consider the integral:

11x2dx \int_1^\infty \frac{1}{x^2} dx

Using the approach outlined previously, we rewrite as a limit:

limb1b1x2dx=limb1x1b=limb11b=1 \lim_{b \to \infty} \int_1^b \frac{1}{x^2} dx = \lim_{b \to \infty} -\frac{1}{x} \bigg |_1^b = \lim_{b \to \infty} 1 - \frac{1}{b} = 1

The p-integral test

For an improper integral of the function xpx^{-p}, there are three possible cases.

Value of ppConvergent integralDivergent integral
p>1p > 1adxxp\displaystyle \int_a^\infty \frac{dx}{x^p}0adxxp\displaystyle \int_0^a \frac{dx}{x^p}
p<1p < 10adxxp\displaystyle \int_0^a \frac{dx}{x^p}adxxp\displaystyle \int_a^\infty \frac{dx}{x^p}
p=1p = 1NoneAll integrals diverge

Arc length

We will cover one last application of integration: finding the total length (called the arc length) of a curve. For a curve y=f(x)y = f(x), an infinitesimal segment can be found by the Pythagorean theorem:

ds2=dx2+dy2 ds^2 = dx^2 + dy^2

The right-hand side can be expressed purely in terms of dxdx by factoring (here note that dy2=(dy)(dy)dy^2 = (dy)(dy) and similarly for dxdx):

ds2=dx2(1+dy2dx2) ds^2 = dx^2 \left(1 + \frac{dy^2}{dx^2}\right)

Square rooting both sides, we have:

ds=1+(dydx)2dx ds = \sqrt{1 + \left(\frac{dy}{dx}\right)^2} dx

Now taking the integral of this infinitesimal arc length, we find the length SS of the complete curve:

S=ds=ab1+(dydx)2dx S = \int ds = \int_a^b \sqrt{1 + \left(\frac{dy}{dx}\right)^2} dx

Finding the surface area of a surface of revolution from arc length

We may combine our technique for finding the arc length and the methods we used for finding volumes of solids of revolution to find their surface areas instead. To start, we can write out the area of an infinitesimal portion of a rotationally-symmetric surface by:

dA=2πrds dA = 2\pi r ds

where dsds is the infinitesimal arc length, and r=f(x)r = f(x). The complete surface area is given by:

A=dA A = \int dA

Therefore, the surface area of a surface of revolution is given by:

A=ab2πrds=ab2πf(x)1+(dydx)2dx A = \int_a^b 2\pi r ds = \int_a^b 2\pi f(x) \sqrt{1 + \left(\frac{dy}{dx}\right)^2} dx

Show table of contents Back to all notes