Series 215 problem 05.mws

Integrals can be used to estimate the remainder of a convergent series. For instance, suppose I wish to evaluate the sum of the p-series, . Euler claims that this sums to We compare the series to , noting that f is continuous, positive and decreasing on [ ) and use the remainder estimate to determine a value of n so that our sum is within .0001 of the actual sum.

> f:=x->1/x^2;

> int(f(x),x=10001..infinity);

> evalf(sum(1/n^2,n=1..10000));

> evalf(Pi^2/6);

Notice that the two numbers agree within the error of our estimate.

Submission:

(a) Find the partial sum of the series . Estimate the error in using as an

approximation to the sum of the series.

(b) Use the fact that

and +

with to give an improved estimate of the sum.

(c) Find a value of n so that is within 0.00001 of the true sum.

(d) How many terms of the series would you need to add to find its sum to within 0.01 ?

Submission worksheet:

>

2. Estimating the sum of a positive series using the integral test.

Let be a positive series and suppose that for some function f defined on [ ) which is continuous, positive and decreasing on that interval. Then we know that the series converges precisely when the integral converges. Suppose that the series does converge and that . Let . Then , and we can show using a picture that

<= <=

This formula allows us to determine bounds on the error in using to estimate the value A to which the series converges.

Let us consider the p -series , which the mathematician Euler discovered converges to . We can compare this series to the function .

It is easy to see that this function is positive, continuous and decreasing on [ ). Let us compute the integral

= = =

From this it follows that <= <= , Thus if we want to estimate the sum of the series to within .001, we need to use 1000 terms in the sequence. Let us compare our estimate with Euler's exact answer.

> est:=evalf(sum(1/n^2,n=1..1000));

> ans:=evalf(Pi^2/6);

> ans-est;

We could also improve our estimate slightly by using the fact that we know that our error is at least 1/1001, so we could refine the estimate to be

> est:=evalf(sum(1/n^2,n=1..1000)+1/1001);

> ans-est;

Notice the big improvement in the estimate we just obtained. This suggests that if we use to approximate the sum of the series, we obtain a better estimate of its sum.

Submission:

The Riemann Zeta Function plays an important role in mathematics. The question of where its zeros are located is the subject of the Riemann Hypothesis, an important unsolved conjecture. The value of the Riemann Zeta function is given by when . The objective of this exercise is to estimate the value of the Riemann Zeta function at using the integral test, in other words, to estimate the sum of the p- series for .

a) Define the term for this series in Maple, and then compute the sum of the series using Maple. Finally, get a floating point estimate of this value.

b) Find the function f(x) which we use in the integral test, and show it satisfies the hypotheses of the integral test.

c) Find a number N which is such that if , the error in using the value to approximate the value of the series is no worse than .0001.

d) Determine the error in using to approximate the sum of the series.

Submission worksheet:

>

>

>