Chapter 7	Techniques of Integration

7.1  Numerical Integration

Students should read Section 7.1 of Rogawski’s Calculus [1] for a detailed discussion of the material presented in this section.

Numerical integration is the process of approximating a definite integral using appropriate sums of function values.  We already saw in Chapter 5 of this text formulas for Right, Left, and Midpoint Rules and their subroutines LRSUM, RRSUM, and MRSUM, respectively.  In this section, we will develop two additional rules: the Trapezoidal Rule and Simpson’s Rule.

7.1.1  Trapezoidal Rule

The Trapezoidal Rule approximates the definite integral (∫_a^b)f(x)dx by using areas of trapezoids and is given by the formula:

(Subscript[T, n])  =  .5 ((b-a)/n )((Subscript[y, 0]) +2(Subscript[y, 1]) + ...  + 2(Subscript[y, n-1]) + (Subscript[y, n]))

where n is the number of trapezoids and (Subscript[y, i])=f(a+i (b-a)/n) .  This formula can be found in your calculus text.  Here is a Mathematica subroutine, called TRAP, for implementing the Trapezoidal Rule:

Clear[f,a,b,n]

TRAP[a_,b_,n_] := (f[a] + 2 Sum[f[a+i*(b-a)/n],{i,1,n-1}] + f[b])(.5 (b-a)/n)

Example 7.1.  Calculate the area under the function f(x)=(x^2) on [0,1] using the Trapezoidal Rule for various values of n.

Solution: The following output gives a table of approximations of (∫_0^1)(x^2)dx based on the Trapezoidal Rule for n=10,20,...,100.

f[x_]:=(x^2)
TableForm[Table[{n,N[TRAP[0,1,n]]},{n,10,100,10}],TableHeadings->{{}, {"n", " (Subscript[T, n])"}}]

It is clear that these values are converging to 1/3, which is the exact value of our definite integral:

(∫_0^1)(x^2)\[DifferentialD]x

7.1.2  Simpson’s Rule

One difference between Simpson’s Rule and all the other rules we have developed so far (TRAP, LRSUM, RRSUM, and MRSUM) is that the number of partition points, n, in this case, must be even.  The other difference is that Simpson’s Rule is a quadratic approximation based on parabolas, whereas the other rules are linear approximations based on lines.  The formula for Simpson’s Rule is given by (refer to your calculus text for details):

(Subscript[S, n])= (1/3) [(Subscript[y, 0])+ 4(Subscript[y, 1])+2(Subscript[y, 2])+4(Subscript[y, 3])+2(Subscript[y, 4])+...+4(Subscript[y, n-3])+2(Subscript[y, n-2])+4(Subscript[y, n-1])+(Subscript[y, n])] (b-a)/n
     =(1/3)[((Subscript[y, 0])+ 4(Subscript[y, 1])+(Subscript[y, 2]))+((Subscript[y, 2])+4(Subscript[y, 3])+(Subscript[y, 4]))+...+((Subscript[y, n-2])+4(Subscript[y, n-1])+(Subscript[y, n]))] (b-a)/n

where (Subscript[y, i])=f(a+i (b-a)/n).  Here is a Mathematica subroutine, called SIMP, for implementing Simpson’s Rule:

Clear[a,b,n]

SIMP[a_,b_,n_] := (1/3)Sum[f[a + (2 i -2)(b-a)/n] + 4 f[a + (2i-1)(b-a)/n] + f[a + 2i (b-a)/n],{i, 1, n/2}] (b-a)/n

Example 7.2.  Calculate the area under the function f(x)=(x^2) on [0,1] using Simpson’s Rule for various values of n.

Solution: We use the same set of values of n as in the previous example.  This will allow us to compare Simpson’s Rule with the Trapezoidal Rule.

f[x_]:=(x^2)
TableForm[Table[{n,N[SIMP[0,1,n]]},{n,10,100,10}],TableHeadings->{{}, {"n", " (Subscript[s, n])"}}]

Notice how fast SIMP converges to the actual value of the integral (1/3) compared to TRAP.

Example 7.3.  Calculate the definite integral of f(x)=sin(25(x^2)) on [0,1] using Simpson’s Rule and approximate it to five decimal places.  What is the minimum number of partition points needed to obtain this level of accuracy?

Solution: We first evaluate SIMP using values for n in increments of 20.

f[x_]:=Sin[25(x^2)]
TableForm[Table[{n,N[SIMP[0,1,n]]},{n,20,200,20}], TableHeadings-> {{},{"n", " (Subscript[s, n])"}}]

Based on the output our approximation, accurate to five decimal places, is 0.10558.  This first occurs between n=100 to n=120.  We evaluate SIMP inside this range to zoom in on the minimum number of partition points needed.

f[x_]:=Sin[25(x^2)]
TableForm[Table[{n,N[SIMP[0,1,n]]},{n,100,120,2}], TableHeadings-> {{},{"n", " (Subscript[s, n])"}}]

Thus, we see that the minimum number of points needed is n=116.  How does this compare with the minimum number of points needed by TRAP to obtain the same level of accuracy?

NOTE: Observe that SIMP does not converge as fast in this example as in the previous example.  This is because the function f(x)=sin(25(x^2)) is oscillatory as the following graph demonstrates:

Plot[f[x],{x,0,1}]

Try increasing the frequency of this function, say to sin(100(x^2)), to see how well SIMP performs.

7.1.3  Midpoint Rule

Since most calculus texts include again the Midpoint Rule in the section on numerical integration, for completeness, we will too.  The Riemann sum using the midpoints of each subinterval is given by the following formula:

Clear[f]
MRSUM[a_,b_,n_]:=Sum[f[a+(i-1/2)*(b-a)/n]*(b-a)/n,{i,1,n}]

Example 7.4.  Calculate the area under the function f(x)=(x^2) on [0,1] using the Midpoint Rule for various values of n.

Solution:

f[x_]:=(x^2)
TableForm[Table[{n, N[MRSUM[0,1,n]]},{n,10,100,10}], TableHeadings-> {{},{"n", "Midpoint rule"}}]

Exercises

1. Consider the definite integral (∫_1^2)ln(x)dx.
a) Using the Trapezoidal Rule, Simpson's Rule, and Midpoint Rule, approximate this integral for n=10,20,...,100.
b) Compare how fast each subroutine (TRAP, SIMP, MRSUM) converges to (∫_a^b)ln(x)dx and decide which of these rules is "best."

2. Repeat Exercise 1 for the following definite integrals:
a)  (∫_0^2)((e^x)/(x+1))dx   		b)  (∫_0^1)cos((x^2))dx		c)  (∫_0^1)(e^(x^2))dx		
Can you make any general conclusions about which rule (Trapezoidal, Simpson’s, Midpoint) is best?

3.  For each of the functions given below, set up a definite integral for the volume of the solid of revolution obtained by revolving the region under f(x) along the given interval and about the given axis.  Then use the subroutines TRAP, SIMP, and MRSUM to approximate the volume of each solid accurate to two decimal places (use various values of n to obtain the desired accuracy).
a) f(x) = cos x;    [0, π/2];    x-axis           b) f(x) = (e^(-(x^2)));    [0,1];   y-axis	c) f(x)=sin x,  [0,π], x-axis

7.2  Techniques of Integration

Students should read Sections 7.2 trhough 7.4 and 7.6 of Rogawski's Calculus [1] for a detailed discussion of the material presented in this section.

All calculus texts have at least a chapter devoted to "Techniques of Integration."  When using Mathematica, these techniques are usually not necessary since Mathematica automatically gives you the answer.

7.2.1 Substitution

On occasion, we do need to use techniques of integration, even when using Mathematica.

Example 7.5.  Evaluate the following integral: ∫(2^x)(Sqrt[((((2^x)))^2)-1])dx.

Solution: We evaluate this integral in Mathematica:

∫(2^x)(Sqrt[((((2^x)))^2)-1])\[DifferentialD]x

To students in a first-year calculus course, this answer makes no sense.  There are many integrals that Mathematica cannot evaluate at all, or cannot evaluate in terms of elementary functions (such as the integral above).  Some of these integrals are doable in terms we should understand, once we first use an appropriate technique of integration. In the above example, all we need to do is first make the following substitution:  u = (2^x)  and d u = (ln 2) (2^x) d x, which transforms the integral to:

(1/(Log[2]))∫(Sqrt[(u^2)-1])\[DifferentialD]u

This is the correct answer.  All we need to do is substitute (2^x) for u, and add the arbitrary constant of integration, getting:
(1/(2 Log[2]))( (2^x)(Sqrt[-1 + ((((2^x)))^2)]) - Log[(2^x) + (Sqrt[-1 + ((((2^x)))^2)])] ) + C
Note that the Mathematica function Log[x] is equivalent to the standard form ln x.

7.2.2 Trigonometric Substitution

Example 7.6.  Evaluate ∫(1/((x^2)(Sqrt[(x^2)-9])))dx.

Solution: By hand, the integral ∫(1/((x^2)(Sqrt[(x^2)-9])))dx would normally be evaluated with a trigonometric substitution of the form x=3secθ.  But with Mathematica, we can do this directly:

∫(1/((x^2)(Sqrt[(x^2)-9])))\[DifferentialD]x

This, of course, is the correct answer, when we remember that Mathematica does not add an arbitrary constant to indefinite integrals.

7.2.3 Method of Partial Fractions

Integrals of rational expressions often require the Method of Partial Fraction Decomposition to evaluate them (by hand).  For example:

∫((3x-3)/((x^2)+5x+4))dx=∫((5/(x+4))-(2/(x+1)))dx=5ln|x+4|-2ln|x+1|=ln|((((x+4))^5)/(((x+1))^2))|

On the other hand, Mathematica will give us essentially the same answer for this integral, but does its work behind the scenes without revealing its technique:

Simplify[∫((3x-3)/((x^2)+5x+4))\[DifferentialD]x]

If we would like to see the partial fraction decomposition of the integrand, ((3x-3)/((x^2)+5x+4)), Mathematica will also do that for us without strain by using the Apart command:

Apart[((3x-3)/((x^2)+5x+4))]

Example 7.7.  Evaluate ∫((2(x^3)+(x^2)-2x+2)/((((x^2)+1))^2))dx.

Solution:  We simply evaluate this integral using Mathematica:

∫((2 (x^3)+ (x^2) -2 x + 2)/((((x^2)+ 1))^2))\[DifferentialD]x

But again, if we would like to see the partial fraction decomposition of the integrand, ((2(x^3)+(x^2)-2x+2)/((((x^2)+1))^2)), then this is straightforward with Mathematica:

Apart[((2 (x^3)+ (x^2) -2 x + 2)/((((x^2)+ 1))^2))]

Exercises

1.  Evaluate ∫(1+ln(x))(Sqrt[1+(((x ln(x)))^2)])dx with Mathematica. If it doesn't give an understandable answer, use a technique of integration that changes the integral into one that Mathematica will evaluate.

In Exercises 2 through 5, use Mathematica to find the partial fraction decomposition of the given functions and then integrate them:
	2.  (((x^2)+3 x - 44)/((x-3)(x+5)(3x-2)))     		3.  ((3(x^2)-4x+5)/((x-1)((x^2)+1)))     	 4.  (25/(x ((x^2)+2x+5))) 		5.  (10/((x((x^2)+2x+5))^2))

In Exercises 6 through 10, use Mathematica to evaluate the given  integrals.
	6. ∫((x^2)/((((x^2)- 4))^(3/2)))dx		7. ∫(x^3)(Sqrt[9-(x^2)])dx		8. ∫(1/(Sqrt[25+(x^2)]))dx	
	9. ∫(sin^5)xdx			9. ∫(((tan^(-1))t)/(1+(t^2)))dt			10. ∫(sinh^3)x cosh xdx

7.3  Improper Integrals

Students should read Section 7.7 of Rogawski's Calculus [1] for a detailed discussion of the material presented in this section.

Recall that there are two types of improper integrals.

Type I: If we assume that f(x) is integrable over [a,b] for all b≥a, then the improper integral of f(x) over [a,∞) is defined as
	(∫_a^∞)f(x)dx = (lim_(t-> ∞))(∫_a^t)f(x)dx,
provided this limit exists.  Similarly, we define
	(∫_(-∞)^b)f(x)dx = (lim_(t-> -∞))(∫_t^b)f(x)dx,
provided this limit exists.

Type II: If f(x) is continuous on [a,b) but discontinuous at x=b, we define
	(∫_a^b)f(x)dx = (lim_(t-> b-))(∫_a^t)f(x)dx ,
provided this limit exists.  Similarly, if f(x) is continuous on (a,b] but discontinuous at x=a,
	(∫_a^b)f(x)dx = (lim_(t-> a+))(∫_t^b)f(x)dx,
provided this limit exists.  Finally, if f(x) is continuous for all x on [a,b] except at x=c, where a<c<b, we define
	(∫_a^b)f(x)dx=(lim_(t->c-))(∫_a^t)f(x)dx+(lim_(t->c+))(∫_t^b)f(x)dx,
provided both of these limits exist.

By using the Limit command in Mathematica along with Integrate, Mathematica eliminates the drudgery of having to evaluate these integrals by hand.

Example 7.8.  Evaluate the following improper integrals:
a)  (∫_20^∞)(1/y)dy
b) (∫_2^∞)(e^(-2x))dx
c) (∫_0^1)x ln xdx
d) (∫_(-∞)^∞)(1/(1+(x^2)))dx

Solution:

a)  We evaluate

(∫_20^∞)(1/y)\[DifferentialD]y

Thus, evaluating this integral directly using Mathematica tells us it does not exist.  Alternatively, we could have used the limit definition:

Limit[(∫_20^t)(1/y)\[DifferentialD]y,t-> ∞]

Observe the difference in the two outputs above.  Both correctly express the answer as divergent; however, the second answer is better since it reveals the nature of the divergence (infinity), which is the answer we would expect if solving this problem by hand.

b)  We evaluate

(∫_2^∞)(E^(-2x))\[DifferentialD]x

Again, we obtain the same answer using the limit definition (as it should):

Limit[(∫_2^t)(E^(-2x))\[DifferentialD]x , t-> ∞]

Mathematica will similarly handle discontinuities. In the following example, the function has a discontinuity at x = 0.

c) We evaluate

(∫_0^1)x Log[x]\[DifferentialD]x

Limit[(∫_t^1)x Log[x]\[DifferentialD]x, t-> 0, Direction->-1]

d)  We evaluate

(∫_(-∞)^∞)(1/(1+(x^2)))\[DifferentialD]x

Note that Mathematica does not require us to break the integral up into two integrals, which would be required according to its definition, if evaluated by hand.  On the other hand, there is nothing wrong with dividing this integral into two in Mathematica:

(∫_(-∞)^0)(1/(1+(x^2)))\[DifferentialD]x + (∫_0^∞)(1/(1+(x^2)))\[DifferentialD]x

NOTE: Observe that it does not matter where we divide the integral.  It is valid to express (∫_(-∞)^a)(1/(1+(x^2)))dx+(∫_a^∞)(1/(1+(x^2)))dx for the integral  (∫_(-∞)^∞)(1/(1+(x^2)))dx for any real value a as long as they are convergent.  However, evaluating this sum in Mathematica yields different expressions for the answer, which depend on the sign of a and whether it is real or complex.  This is shown in the following output:

Clear[a]
(∫_(-∞)^a)(1/(1+(x^2)))\[DifferentialD]x + (∫_a^∞)(1/(1+(x^2)))\[DifferentialD]x

If instead, a is given a fixed value, then Mathematica will give us our answer of π:

a=1
(∫_(-∞)^a)(1/(1+(x^2)))\[DifferentialD]x + (∫_a^∞)(1/(1+(x^2)))\[DifferentialD]x

Exercises

In Excercises 1 through 8, evaluate the given improper integrals:
1.  (∫_(-∞)^4)(e^(.01 t))dt 		2.  (∫_(-3)^∞)(1/(((x+4))^(3/2)))dx 		3.  (∫_(-2)^4)(1/(((x+2))^(1/3)))dx 		4.  (∫_(-∞)^∞)x (e^(-(x^2)))dx		
5. (∫_0^3)(1/(x-1))dx			6. (∫_(-∞)^∞)(1/((e^x)+(e^(-x))))dx		7. (∫_1^∞)(1/(x^(.999)))dx			8.  (∫_1^∞)(1/(x^1.003))dx

11.  Find the volume of the solid obtained by rotating the region below the graph of  y=(e^(-x)) about the x-axis for  0≤x<∞.

12.  Determine how large the number b has to be in order that   (∫_b^∞)(1/((x^2)+1))dx <  .0001.

13. Evaluate the improper integral (∫_(-1)^1)(1/(|(x^(1/3))|))dx.

14.  Determine how large the number b should be so that  (∫_b^∞)(1/((x^2)+1))dx   <   .0001.

15. Consider the function defined by

G(x)=(∫_0^∞)(t^(x-1))(e^(-t))dt

a)  Evaluate G(n)  for n=0, 2,,3, 4, ...., 10.  Make  a conjecture about these values.  Verify your conjecture.
b) Evaluate G((2n-1)/2), for n=1, 2, 3, ... 10.  Make a conjecture about these values.  Verify your conjecture.
c) Plot the graph of G(x) on the interval [0, 5].

NOTE: The function G is called the gamma function and is denoted by Γ[x]. In Mathematica it is denoted by Gamma[x].   The gamma function was first introduced by Euler as a generalization of the factorial function.

7.4  Hyperbolic and Inverse Hyperbolic Functions

Students should read Section 7.5 of Rogawski's Calculus [1] for a detailed discussion of the material presented in this section.

7.4.1. Hyperbolic Functions

The hyperbolic functions are defined in terms of the exponential functions.  They have a direct connection to engineering mathematics, including bridge construction.  For example, cables from suspension bridges typically form a curve called a catenary (derived from the Latin word catena, which means chain) that is described by these functions.

The six hyperbolic functions are denoted and defined as follows:

sinh x =(((e^x)- (e^(-x)))/2), 	 cosh x=(((e^x)+ (e^(-x)))/2),	tanh x=(((e^x)- (e^(-x)))/((e^x)+ (e^(-x))))

coth x =(((e^x)+ (e^(-x)))/((e^x)- (e^(-x)))), 	 sech x=(2/((e^x)+ (e^(-x)))),	csch x=(2/((e^x)- (e^(-x))))

The reason these functions are called hyperbolic functions is due to their connection with the equilateral hyperbola (x^2)-(y^2)=1.  Here, one defines x=cosh t and y=sinh t.  Hence, one obtains the basic hyperbolic identity (cosh^2)t-(sinh^2)t=1, much the same manner as the corresponding trigonometric identity (cos^2)t+(sin^2)t=1, when one considers the unit circle (x^2)+(y^2)=1 with x=cos t and y=sin t.

In Mathematica, we use the same notation with the obvious convention that the first letter of each function is capitalized and square brackets must be used in place of parentheses.  Thus, sinh x will be entered as Sinh[x].

Example 7.9.  Consider the hyperbolic sine function f(x)=sinh x.
a) Plot the graph of f. 
b) From the graphs deduce the domain and range of the function. 
c) Is f bounded? 
d) Does f attain an absolute minimum? Maximum?
e) Repeat a) through d) for the hyperbolic function g(x)=cosh x  
f) Repeat a) through d) for the hyperbolic function h(x)=tanh x.

Solution: We begin by defining f in Mathematica:

Clear[f,x]
f[x_]=Sinh[x]

a) We next plot its graph on the interval [-3,3].

Plot[f[x], {x, -3,3}]

b) The preceding graph indicates that the domain and range of sinh x is (-∞,∞). To convince yourself, you should plot the graph over wider intervals.  We should also expect this from the definition of sinh x itself.  Can you explain why?

c) The function sinh x is not bounded. The graph earlier should not be used as a proof of this.  However, we can evaluate its limit at -∞ and ∞ to see that this is indeed true.

Limit[f[x],x->-∞]
Limit[f[x],x->∞]

d) The limits just computed show that sinh x has no absolute maximum or minimum since it is unbounded.

e) Next, we consider the hyperbolic cosine function denoted by cosh x.

Clear[g,x]
g[x_]=Cosh[x]

Plot[g[x],{x,-3,3}]

The preceding graph indicates that the domain of cosh x is (-∞,∞). The range appears to be  [1,∞). Can you prove this?

The hyperbolic cosine function, cosh x, is not bounded from above. This can be seen from the following limits:

Limit[Cosh[x],x->-∞]
Limit[Cosh[x],x->∞]

Again, since cosh x is not bounded from above, it follows that cosh x has no absolute maximum. As we have observed in part b) of this example, cosh x has absolute minimum value 1, attained at x=0.

f) Finally, we consider the hyperbolic tangent function, tanh x:

Clear[h,x]
h[x_]=Tanh[x]

Plot[h[x],{x,-3,3}]

Again, the preceding graph indicates that the domain of tanh x is (-∞,∞). The range appears to be  (-1,1). This can be seen from the following limits:

Limit[Tanh[x],x->-∞]
Limit[Tanh[x],x->∞]

The graph of tanh x also indicates that it is strictly increasing on its domain.  This can be proven by showing that its derivative, which we will calculate later, is strictly positive.  It is clear that tanh x has no absolute extrema.

NOTE: The reader will notice some similarities between the hyperbolic functions and the associated trigonometric functions. Moreover, if one studies the theory of functions of a complex variable, the relationship between these classes of transcendental functions becomes even more transparent; for numerous identities exist between the classes of functions.

7.4.2 Identities Involving Hyperbolic Functions

It is immediate that the ratio and reciprocal identities for the hyperbolic functions coincide with their trigonometric counterparts. In fact, for each trigonometric identity, there is a corresponding (not necessarily the same) hyperbolic identity.  Following are some examples.

Example 7.10.  Show that the following identities hold true.
a) 1-(tanh^2)x=(sech^2)x	b)  cosh(x+y)=cosh x cosh y+sinh x sinh y

Solution:

a) We use the definitions for tanh x and sech x to express each side of the identity in terms of exponentials:

Simplify[(1-((Tanh[x])^2))/.Tanh[x]->(E^x-E^(-x))/(E^x+E^(-x))]

Simplify[((Sech[x])^2)/.Sech[x]->2/(E^x+E^(-x))]

We leave it for the reader to verify that both of these outputs agree, that is, ((4(e^(2x)))/(((1+(e^(2x))))^2))=(4/((((e^(-x))+(e^x)))^2)) (cross-multiply and then simplify).

The identity can also be confirmed in Mathematica by evaluating the difference between its left- and right-hand sides, which should equal zero:

Simplify[1-((Tanh[x])^2)-((Sech[x])^2)]

NOTE: We can also confirm the identity graphically by plotting the graphs of each side of the identity, which should coincide.

Plot[{1-Tanh[x]^2,Sech[x]^2},{x,-2,2}]

b) We again evaluate the difference between the left- and right-hand sides of the identity:

Simplify[Cosh[x+y]-(Cosh[x] Cosh[y]+Sinh[x] Sinh[y])]

7.4.3 Derivatives of Hyperbolic Functions

We next contrast the formulas for the derivatives of the trigonometric functions versus the formulas for the derivatives of the companion hyperbolic functions.

Example 7.11.  Compare the derivatives of the given pair of functions.
a) sinh x and sin x	b)  cosh x and cos x	 c) tanh x and tan x

Solution: We use the derivative command, D, to evaluate derivatives of each pair.

a)

D[Sinh[x],x]
D[Sin[x],x]

b)

D[Cosh[x],x]
D[Cos[x],x]

b)

D[Tanh[x],x]
D[Tan[x],x]

It is clear that derivatives of hyperbolic and trigonometric functions are quite similar.

7.4.4 Inverse Hyperbolic Functions

In light of the fact that hyperbolic functions are defined in terms of the exponential functions, it is readily apparent that the inverse hyperbolic functions are defined in terms of the natural logarithmic function.  The inverses of the hyperbolic functions have notation similar to those of inverse trigonometric functions.  Thus, the inverse of sinh x is denoted by arcsinh x or (sinh^(-1)) x.  In Mathematica, the notation is (sinh^(-1))x is ArcSinh[x].

Example 7.12.  Plot the graphs of (sinh^(-1))x and sinh x on the same axis.

Solution: Recall that the graph of a function and the graph of its inverse are reflections of each other across the line y=x.  This is confirmed by the following plot of (sinh^(-1))x (in blue) and sin x (in red).

Plot[{Sinh[x], x, ArcSinh[x]},{x, -3,3}, PlotStyle->{Blue,Green,Red},AspectRatio->Automatic,PlotRange->{-3,3}]

Example 7.13.  Show that (tanh^(-1))x=(1/2)ln(((1+x)/(1-x))) for -1<x<1.

Solution: We plot the graphs of y=(tanh^(-1))x and y=(1/2)ln(((1+x)/(1-x))) on the same axes. Note that Mathematica's notation of (tanh^(-1))x is ArcTanh [x] and ln y is entered as Log[y]:

Plot[{ArcTanh[x],(1/2)Log[((1+x)/(1-x))]},{x,-2,2}]

The fact that there is only one graph indicates that the functions are the same. We prove this by letting y=(tanh^(-1))x and solving for y as follows.  From y=(tanh^(-1))x we get x=tanh y=(((e^y)-(e^(-y)))/((e^y)+(e^(-y)))).  Now solving this last equation for y in Mathematica yields:

Solve[x==(E^y-E^(-y))/(E^y+E^(-y)),y]

The first solution in the preceding output is imaginary, which we ignore, and consider only the second solution.  Hence,

(tanh^(-1)) x =y=ln((Sqrt[-1-x])/(Sqrt[-1+x]))=ln(Sqrt[((1+x)/(1-x))])=(1/2)ln(((1+x)/(1-x))).

NOTE: The message in the previous output refers to the fact that when solving equations involving inverse functions, not all solutions are necessarily found by Mathematica since there may be infinitely many of them or they depend on the domain of definition.  For example, the equation sin x=1 has infinitely many solutions, in particular all values of the form x=π/2+2π n, where n is any integer.  On the other hand, solving this equation in Mathematica yields only the solution in its principal domain, that is, x=π/2:

Solve[Sin[x]==1,x]

Exercises

In Exercises 1 through 5, verify the given hyperbolic identities using the Simplify command. Also state the corresponding trigonometric identity.

1.  sinh(x+y)=sinh x cosh y+cosh x sinh y 			2.  cosh 2x = (cosh^2)x+(sinh^2)x		3.   tanh 2x=((2 tanh x)/(1+(tanh^2)x)) 
4.  cosh(x+y)=cosh x cosh y + sonh x sinh y		5.  tanh(x+y)=((tanh x+tanh y)/(1+tanh x tanh y))

6.  Determine the first few positive integral powers of cosh x + sinh x. Can you form a general conjecture for the nth case, namely (((cosh x+ sinh x))^n), where n is any natural number?  Then justify your conclusion via mathematical induction.

In Exercises 7 through 12, determine the derivatives of thegiven functions and simplify your answers where possible. Compare your solution via paper and pencil methods with the one generated by Mathematica.

7.  f(x)=tanh (1+(x^2)) 			8.  f(x)=x sinh x - cosh x					9. f(x)=(Sqrt[((1+tanh x)/(1- tanh x))])
10. f(x)=(x^2)(sinh^(-1))(2x)	 		11.  f(x)=x (tanh^(-1))x+ ln ((Sqrt[1-(x^2)]))			12. f(x)=x coth x -sech x

13. The Gateway Arch in St. Louis was designed by Eero Saarinen and was constructed using the equation

y = 211.49 - 20.96 cosh (0.03291765 x)

for the central curve of the arch, where x and y are measured in meters and |x|≤ 91.20.

a) Plot the graph of the central curve.
b) What is the height of the arch at its center?
c) At what points is the arch 100 meters in height?
d) What is the slope of the arch at the points in part (c)?

14.  A flexible cable always hangs in the shape of a catenary y = c + a cosh (x / a), where c and a are constants and a > 0.  Plot several members of the family of functions y = a cosh (x / a) for various values of a.  How does the graph change as a varies?

In Exercises 15 through 17,  evaluate each of the given integrals:

15.  ∫sinh x (cosh^n)xdx			16. ∫((cosh x)/((cosh^2)x-1))dx			17.  ∫(((sech^2)x)/(2+tanh x))dx

18.   Let t=ln(((1+(Sqrt[5]))/2)) and define

f(n)={({{(2/(Sqrt[5]))cosh(t n), if n is odd}, {(2/(Sqrt[5]))sinh(t n), if n is even}})

Evaluate f(n) for n=1, 2, 3, ...,20.  Do these values seem familiar? If not, we highly recommend the interesting article by Thomas Osler, Vieta-like products of nested radicals with Fibonacci and Lucas numbers, to appear in the journal Fibonacci Quarterly.