##
Section 5.5 Other Options for Finding Algebraic Antiderivatives

¶###### Motivating Questions

How does the method of partial fractions enable any rational function to be antidifferentiated?

What role have integral tables historically played in the study of calculus and how can a table be used to evaluate integrals such as \(\int \sqrt{a^2 + u^2} \, du\text{?}\)

What role can a computer algebra system play in the process of finding antiderivatives?

In the preceding sections, we have learned two very specific antidifferentiation techniques: \(u\)-substitution and integration by parts. The former is used to reverse the chain rule, while the latter to reverse the product rule. But we have seen that each only works in very specialized circumstances. For example, while \(\int xe^{x^2} \, dx\) may be evaluated by \(u\)-substitution and \(\int x e^x \, dx\) by integration by parts, neither method provides a route to evaluate \(\int e^{x^2} \, dx\text{.}\) That fact is not a particular shortcoming of these two antidifferentiation techniques, as it turns out there does not exist an elementary algebraic antiderivative for \(e^{x^2}\text{.}\) Said differently, no matter what antidifferentiation methods we could develop and learn to execute, none of them will be able to provide us with a simple formula that does not involve integrals for a function \(F(x)\) that satisfies \(F'(x) = e^{x^2}\text{.}\)

In this section of the text, our main goals are to better understand some classes of functions that can always be antidifferentiated, as well as to learn some options for so doing. At the same time, we want to recognize that there are many functions for which an algebraic formula for an antiderivative does not exist, and also appreciate the role that computing technology can play in helping us find antiderivatives of other complicated functions. Throughout, it is helpful to remember what we have learned so far: how to reverse the chain rule through \(u\)-substitution, how to reverse the product rule through integration by parts, and that overall, there are subtle and challenging issues to address when trying to find antiderivatives.

######
Preview Activity 5.5.1

For each of the indefinite integrals below, the main question is to decide whether the integral can be evaluated using \(u\)-substitution, integration by parts, a combination of the two, or neither. For integrals for which your answer is affirmative, state the substitution(s) you would use. It is not necessary to actually evaluate any of the integrals completely, unless the integral can be evaluated immediately using a familiar basic antiderivative.

\(\int x^2 \sin(x^3) \, dx\text{,}\) \(\int x^2 \sin(x) \, dx\text{,}\) \(\int \sin(x^3) \, dx\text{,}\) \(\int x^5 \sin(x^3) \, dx\)

\(\int \frac{1}{1+x^2} \, dx\text{,}\) \(\int \frac{x}{1+x^2} \, dx\text{,}\) \(\int \frac{2x+3}{1+x^2} \, dx\text{,}\) \(\int \frac{e^x}{1+(e^x)^2} \, dx\text{,}\)

\(\int x \ln(x) \, dx\text{,}\) \(\int \frac{\ln(x)}{x} \, dx\text{,}\) \(\int \ln(1+x^2) \, dx\text{,}\) \(\int x\ln(1+x^2) \, dx\text{,}\)

\(\int x \sqrt{1-x^2} \, dx\text{,}\) \(\int \frac{1}{\sqrt{1-x^2}} \, dx\text{,}\) \(\int \frac{x}{\sqrt{1-x^2}}\, dx\text{,}\) \(\int \frac{1}{x\sqrt{1-x^2}} \, dx\text{,}\)

###
Subsection 5.5.1 The Method of Partial Fractions

The method of partial fractions is used to integrate rational functions, and essentially involves reversing the process of finding a common denominator. For example, suppoes we have the function \(R(x) = \frac{5x}{x^2 - x - 2}\) and want to evaluate

\begin{equation*}
\int \frac{5x}{x^2-x-2} \, dx\text{.}
\end{equation*}

Thinking algebraically, if we factor the denominator, we can see how \(R\) might come from the sum of two fractions of the form \(\frac{A}{x-2} + \frac{B}{x+1}\text{.}\) In particular, suppose that

\begin{equation*}
\frac{5x}{(x-2)(x+1)} = \frac{A}{x-2} + \frac{B}{x+1}\text{.}
\end{equation*}

Multiplying both sides of this last equation by \((x-2)(x+1)\text{,}\) we find that

\begin{equation*}
5x = A(x+1) + B(x-2)\text{.}
\end{equation*}

Since we want this equation to hold for every value of \(x\text{,}\) we can use insightful choices of specific \(x\)-values to help us find \(A\) and \(B\text{.}\) Taking \(x = -1\text{,}\) we have

\begin{equation*}
5(-1) = A(0) + B(-3)\text{,}
\end{equation*}

and thus \(B = \frac{5}{3}\text{.}\) Choosing \(x = 2\text{,}\) it follows

\begin{equation*}
5(2) = A(3) + B(0)\text{,}
\end{equation*}

so \(A = \frac{10}{3}\text{.}\) Therefore, we now know that

\begin{equation*}
\int \frac{5x}{x^2-x-2} \, dx = \int \frac{10/3}{x-2} + \frac{5/3}{x+1} \, dx\text{.}
\end{equation*}

This equivalent integral expression is straightforward to evaluate, and hence we find that

\begin{equation*}
\int \frac{5x}{x^2-x-2} \, dx = \frac{10}{3} \ln|x-2| + \frac{5}{3}\ln|x+1| + C\text{.}
\end{equation*}

It turns out that for any rational function \(R(x) = \frac{P(x)}{Q(x)}\) where the degree of the polynomial \(P\) is less than^{ 1 }the degree of the polynomial \(Q\text{,}\) the method of partial fractions can be used to rewrite the function as a sum of simpler rational functions of one of the following forms:

\begin{equation*}
\frac{A}{x-c},\ \frac{A}{(x-c)^n},\ \frac{Ax+B}{x^2 + k},\ \text{or }\frac{Ax+B}{\left(x^2 + k\right)^n}
\end{equation*}

where \(A\text{,}\) \(B\text{,}\) and \(c\) are real numbers, and \(k\) is a positive real number. Because each of these basic forms is one we can antidifferentiate, partial fractions enables us to antidifferentiate any rational function. A computer algebra system such as *Maple*, *Mathematica*, or *WolframAlpha* can be used to find the partial fraction decomposition of any rational function. In *WolframAlpha*, entering

`partial fraction 5x/(x^2-x-2)`

results in the output

\begin{equation*}
\frac{5x}{x^2-x-2} = \frac{10}{3(x-2)} + \frac{5}{3(x+1)}\text{.}
\end{equation*}

We will primarily use technology to generate partial fraction decompositions of rational functions, and then work from there to evaluate the integrals of interest using established methods.

######
Activity 5.5.2

For each of the following problems, evaluate the integral by using the partial fraction decomposition provided.

\(\int \frac{1}{x^2 - 2x - 3} \, dx\text{,}\) given that \(\frac{1}{x^2 - 2x - 3} = \frac{1/4}{x-3} - \frac{1/4}{x+1}\)

\(\int \frac{x^2+1}{x^3 - x^2} \, dx\text{,}\) given that \(\frac{x^2+1}{x^3 - x^2} = -\frac{1}{x} - \frac{1}{x^2} + \frac{2}{x-1}\)

\(\int \frac{x-2}{x^4 + x^2}\, dx\text{,}\) given that \(\frac{x-2}{x^4 + x^2} = \frac{1}{x} - \frac{2}{x^2} + \frac{-x+2}{1+x^2}\)

###
Subsection 5.5.2 Using an Integral Table

Calculus has a long history, with key ideas going back as far as Greek mathematicians in 400-300 BC. Its main foundations were first investigated and understood independently by Isaac Newton and Gottfried Wilhelm Leibniz in the late 1600s, making the modern ideas of calculus well over 300 years old. It is instructive to realize that until the late 1980s, the personal computer essentially did not exist, so calculus (and other mathematics) had to be done by hand for roughly 300 years. During the last 30 years, however, computers have revolutionized many aspects of the world we live in, including mathematics. In this section we take a short historical tour to precede the following discussion of the role computer algebra systems can play in evaluating indefinite integrals. In particular, we consider a class of integrals involving certain radical expressions that, until the advent of computer algebra systems, were often evaluated using an integral table.

As seen in the short table of integrals found in Appendix A, there are also many forms of integrals that involve \(\sqrt{a^2 \pm w^2}\) and \(\sqrt{w^2 - a^2}\text{.}\) These integral rules can be developed using a technique known as *trigonometric substitution* that we choose to omit; instead, we will simply accept the results presented in the table. To see how these rules are needed and used, consider the differences among

\begin{equation*}
\int \frac{1}{\sqrt{1-x^2}} \,dx, \ \ \ \int \frac{x}{\sqrt{1-x^2}} \,dx, \ \ \ \text{and} \ \ \ \int \sqrt{1-x^2} \,dx\text{.}
\end{equation*}

The first integral is a familiar basic one, and results in \(\arcsin(x) + C\text{.}\) The second integral can be evaluated using a standard \(u\)-substitution with \(u = 1-x^2\text{.}\) The third, however, is not familiar and does not lend itself to \(u\)-substitution.

In Appendix A, we find the rule

\begin{equation*}
(h) ~ \int \sqrt{a^2 - u^2} \, du = \frac{u}{2}\sqrt{a^2 - u^2} + \frac{a^2}{2} \arcsin \frac{u}{a} + C\text{.}
\end{equation*}

Using the substitutions \(a = 1\) and \(u = x\) (so that \(du = dx\)), it follows that

\begin{equation*}
\int \sqrt{1-x^2} \, dx = \frac{x}{2} \sqrt{1-x^2} - \frac{1}{2} \arcsin x + C\text{.}
\end{equation*}

One important point to note is that whenever we are applying a rule in the table, we are doing a \(u\)-substitution. This is especially key when the situation is more complicated than allowing \(u = x\) as in the last example. For instance, say we wish to evaluate the integral

\begin{equation*}
\int \sqrt{9 + 64x^2} \, dx\text{.}
\end{equation*}

Here, we want to use Rule (c) from the table, and do so with \(a = 3\) and \(u = 8x\text{;}\) we also choose the “\(+\)” option in the rule. With this substitution, it follows that \(du = 8dx\text{,}\) so \(dx = \frac{1}{8} du\text{.}\) Applying this substitution,

\begin{equation*}
\int \sqrt{9 + 64x^2} \, dx = \int \sqrt{9 + u^2} \cdot \frac{1}{8} \, du = \frac{1}{8} \int \sqrt{9+u^2} \, du\text{.}
\end{equation*}

By Rule (c), we now find that

\begin{align*}
\int \sqrt{9 + 64x^2} \, dx =\mathstrut \amp \frac{1}{8} \left( \frac{u}{2}\sqrt{u^2 + 9} + \frac{9}{2}\ln|u + \sqrt{u^2 + 9}| + C \right)\\
=\mathstrut \amp \frac{1}{8} \left( \frac{8x}{2}\sqrt{64x^2 + 9} + \frac{9}{2}\ln|8x + \sqrt{64x^2 + 9}| + C \right)\text{.}
\end{align*}

In problems such as this one, it is essential that we not forget to account for the factor of \(\frac{1}{8}\) that must be present in the evaluation.

######
Activity 5.5.3

For each of the following integrals, evaluate the integral using \(u\)-substitution and/or an entry from the table found in Appendix A.

\(\int \sqrt{x^2 + 4} \, dx\)

\(\int \frac{x}{\sqrt{x^2 +4}} \, dx\)

\(\int \frac{2}{\sqrt{16+25x^2}}\, dx\)

\(\int \frac{1}{x^2 \sqrt{49-36x^2}} \, dx\)

###
Subsection 5.5.3 Using Computer Algebra Systems

A computer algebra system (CAS) is a computer program that is capable of executing symbolic mathematics. For a simple example, if we ask a CAS to solve the equation \(ax^2 + bx + c = 0\) for the variable \(x\text{,}\) where \(a\text{,}\) \(b\text{,}\) and \(c\) are arbitrary constants, the program will return \(x = \frac{-b \pm \sqrt{b^2 - 4ac}}{2a}\text{.}\) While research to develop the first CAS dates to the 1960s, these programs became more common and publicly available in the early 1990s. Two prominent early examples are the programs *Maple* and *Mathematica*, which were among the first computer algebra systems to offer a graphical user interface. Today, *Maple* and *Mathematica* are exceptionally powerful professional software packages that are capable of executing an amazing array of sophisticated mathematical computations. They are also very expensive, as each is a proprietary program. The CAS *SAGE* is an open-source, free alternative to *Maple* and *Mathematica*.

For the purposes of this text, when we need to use a CAS, we are going to turn instead to a similar, but somewhat different computational tool, the web-based “computational knowledge engine” called *WolframAlpha*. There are two features of *WolframAlpha* that make it stand out from the CAS options mentioned above: (1) unlike *Maple* and *Mathematica*, *WolframAlpha* is free (provided we are willing to suffer through some pop-up advertising); and (2) unlike any of the three, the syntax in *WolframAlpha* is flexible. Think of *WolframAlpha* as being a little bit like doing a Google search: the program will interpret what is input, and then provide a summary of options.

If we want to have *WolframAlpha* evaluate an integral for us, we can provide it syntax such as

`integrate x^2 dx`

to which the program responds with

\begin{equation*}
\int x^2 \, dx = \frac{x^3}{3} + \text{constant}\text{.}
\end{equation*}

While there is much to be enthusiastic about regarding CAS programs such as *WolframAlpha*, there are several things we should be cautious about: (1) a CAS only responds to exactly what is input; (2) a CAS can answer using powerful functions from highly advanced mathematics; and (3) there are problems that even a CAS cannot do without additional human insight.

Although (1) likely goes without saying, we have to be careful with our input: if we enter syntax that defines a function other than the problem of interest, the CAS will work with precisely the function we define. For example, if we are interested in evaluating the integral

\begin{equation*}
\int \frac{1}{16-5x^2} \, dx\text{,}
\end{equation*}

and we mistakenly enter

`integrate 1/16 - 5x^2 dx`

a CAS will (correctly) reply with

\begin{equation*}
\frac{1}{16}x - \frac{5}{3} x^3\text{.}
\end{equation*}

It is essential that we are sufficiently well-versed in antidifferentiation to recognize that this function cannot be the one that we seek: integrating a rational function such as \(\frac{1}{16-5x^2}\text{,}\) we expect the logarithm function to be present in the result.

Regarding (2), even for a relatively simple integral such as \(\int \frac{1}{16-5x^2} \, dx\text{,}\) some CASs will invoke advanced functions rather than simple ones. For instance, if we use *Maple* to execute the command

`int(1/(16-5*x^2), x);`

the program responds with

\begin{equation*}
\int \frac{1}{16-5x^2} \, dx = \frac{\sqrt{5}}{20} \arctanh (\frac{\sqrt{5}}{4}x)\text{.}
\end{equation*}

While this is correct (save for the missing arbitrary constant, which *Maple* never reports), the inverse hyperbolic tangent function is not a common nor familiar one; a simpler way to express this function can be found by using the partial fractions method, and happens to be the result reported by *WolframAlpha*:

\begin{equation*}
\int \frac{1}{16-5x^2} \, dx = \frac{1}{8\sqrt{5}} \left(\log(4\sqrt{5}+5\sqrt{x}) - \log(4\sqrt{5}-5\sqrt{x})\right) + \text{constant}\text{.}
\end{equation*}

Using sophisticated functions from more advanced mathematics is sometimes the way a CAS says to the user “I don't know how to do this problem.” For example, if we want to evaluate

\begin{equation*}
\int e^{-x^2} \, dx\text{,}
\end{equation*}

and we ask *WolframAlpha* to do so, the input

`integrate exp(-x^2) dx`

results in the output

\begin{equation*}
\int e^{-x^2} \, dx = \frac{\sqrt{\pi}}{2}\erf (x) + \text{constant}\text{.}
\end{equation*}

The function “erf\((x)\)” is the *error function*, which is actually defined by an integral:

\begin{equation*}
\erf (x) = \frac{2}{\sqrt{\pi}} \int_0^x e^{-t^2} \, dt\text{.}
\end{equation*}

So, in producing output involving an integral, the CAS has basically reported back to us the very question we asked.

Finally, as remarked at (3) above, there are times that a CAS will actually fail without some additional human insight. If we consider the integral

\begin{equation*}
\int (1+x)e^x \sqrt{1+x^2e^{2x}} \, dx
\end{equation*}

and ask *WolframAlpha* to evaluate

`int (1+x) * exp(x) * sqrt(1+x^2 * exp(2x)) dx`

,

the program thinks for a moment and then reports

(*no result found in terms of standard mathematical functions*)

But in fact this integral is not that difficult to evaluate. If we let \(u = xe^{x}\text{,}\) then \(du = (1+x)e^x \, dx\text{,}\) which means that the preceding integral has form

\begin{equation*}
\int (1+x)e^x \sqrt{1+x^2e^{2x}} \, dx = \int \sqrt{1+u^2} \, du\text{,}
\end{equation*}

which is a straightforward one for any CAS to evaluate.

So, the above observations regarding computer algebra systems lead us to proceed with some caution: while any CAS is capable of evaluating a wide range of integrals (both definite and indefinite), there are times when the result can mislead us. We must think carefully about the meaning of the output, whether it is consistent with what we expect, and whether or not it makes sense to proceed.

###
Subsection 5.5.4 Summary

The method of partial fractions enables any rational function to be antidifferentiated, because any polynomial function can be factored into a product of linear and irreducible quadratic terms. This allows any rational function to be written as the sum of a polynomial plus rational terms of the form \(\frac{A}{(x-c)^n}\) (where \(n\) is a natural number) and \(\frac{Bx+C}{x^2 + k}\) (where \(k\) is a positive real number).

Until the development of computing algebra systems, integral tables enabled students of calculus to more easily evaluate integrals such as \(\int \sqrt{a^2 + u^2} \, du\text{,}\) where \(a\) is a positive real number. A short table of integrals may be found in Appendix A.

Computer algebra systems can play an important role in finding antiderivatives, though we must be cautious to use correct input, to watch for unusual or unfamiliar advanced functions that the CAS may cite in its result, and to consider the possibility that a CAS may need further assistance or insight from us in order to answer a particular question.

###
Subsection Exercises

¶######
1 Partial fractions: linear over difference of squares

######
2 Partial fractions: constant over product

######
3 Partial fractions: linear over quadratic

######
4 Partial fractions: cubic over 4th degree

######
5 Partial fractions: quadratic over factored cubic

######
6

For each of the following integrals involving rational functions, (1) use a CAS to find the partial fraction decomposition of the integrand; (2) evaluate the integral of the resulting function without the assistance of technology; (3) use a CAS to evaluate the original integral to test and compare your result in (2).

\(\int \frac{x^3 + x + 1}{x^4 - 1} \, dx\)

\(\int \frac{x^5 + x^2 + 3}{x^3 - 6x^2 + 11x - 6} \, dx\)

\(\int \frac{x^2 - x - 1}{(x-3)^3} \, dx\)

######
7

For each of the following integrals involving radical functions, (1) use an appropriate \(u\)-substitution along with Appendix A to evaluate the integral without the assistance of technology, and (2) use a CAS to evaluate the original integral to test and compare your result in (1).

\(\int \frac{1}{x \sqrt{9x^2 + 25}} \, dx\)

\(\int x \sqrt{1 + x^4} \, dx\)

\(\int e^x \sqrt{4 + e^{2x}} \, dx\)

\(\int \frac{\tan(x)}{\sqrt{9 - \cos^2(x)}} \, dx\)

######
8

Consider the indefinite integral given by

\begin{equation*}
\int \frac{\sqrt{x+\sqrt{1+x^2}}}{x} \, dx\text{.}
\end{equation*}

Explain why \(u\)-substitution does not offer a way to simplify this integral by discussing at least two different options you might try for \(u\text{.}\)

Explain why integration by parts does not seem to be a reasonable way to proceed, either, by considering one option for \(u\) and \(dv\text{.}\)

Is there any line in the integral table in Appendix A that is helpful for this integral?

Evaluate the given integral using *WolframAlpha*. What do you observe?