Math 240, Power Series Lab

Radius of Convergence Lab Instructions

In this lab, we will look at Taylor polynomials that approximate the solution of a differential equation. In this case, we won't know the exact solution, so we won't be able to compare our approximations to the actual function value in most cases. However, it usually isn't hard to see where things seem to work, and this will build some intuition about where series converge that we can justify with theorems in Wednesday's lecture. We will also look at a couple of dangers in using Taylor polynomial approximations.

The first problem we will consider is (x²+1)y" + xy' + 2y = 0, y(0)=0, y'(0)=1. Launch the lab (see the link at the top of the page) and change the coefficients to match this equation. Now look at the 5, 10, 25, 50, and 1000 degree Taylor polynomials. What do you think is the radius of convergence of the series solution for this problem?
Now change the coefficients of y and y', along with the initial values. How does changing the lower order coefficients seem to affect the radius of convergence of the power series? (You may find it easiest to check the radius of convergence if you leave the degree set at 1000).
Next consider the problem (.2x²+.4x+1)y" + (x + 1)y' + 3y = 0, y(0)=1, y'(0)=0. Where do you think the series solution for this problem converges? How does changing the lower order terms affect the radius of convergence?.
Based on your work so far, make a conjecture about the relationship between radius of convergence and the coefficients of the differential equation. Test your conjecture with at least two additional initial value problems.
The next two problems point out places where looking at Taylor approximations may be misleading.
A series solution will be valid up to the radius of convergence. The approximation provided by a finite Taylor polynomial (as opposed to the infinite Taylor series) won't work over quite as large a region. In fact, sometimes it can be difficult to impossible to discover features of the solution near the radius of convergence from examining Taylor polynomials. Consider the equation (x²+2x+1)y" + (x+1)y' + y = 0, y(0)=1, y'(0)=0. The exact solution to this initial value problem is y = cos(log(x+1)) (we will see how to find this solution next week). The radius of convergence of this equation 1, so the Taylor series represents the function in the region -1 < x < 1. Look at the behaviour of the Taylor polynomials near x = -1. Sketch how you would guess the solution to the initial value problem behaves near x = -1 based on the behaviour of the Taylor polynomials up to degree 1000. Show that the solution actually oscillates between -1 and 1 infinitely many times near x = -1. Explain why no polynomial can represent infinitely many oscillations (and in this case, it doesn't even come close).
Other problems may arise with the computation of the Taylor polynomials even away from the radius of convergence. Consider the equation y" + y = 0, y(0)=1, y'(0)=0. The true solution to this problem is y = cos(x) of course, and the radius of convergence is infinity. But try to get a Taylor polynomial approximation that works well past about x = 35. Sketch the behaviour of the Taylor polynomial approximations for 35 < x <40. The difficulty you observe is caused by a build up of round-off error in the computation of the Taylor polynomial. Note that the applet uses double precision arithmetic and is coded to minimize the number of multiplications (and hence the amount of accumulated round-off error). Sometimes the computation just doesn't want to behave (and then you try to come up with some means other than blindly computing the polynomial to find the answer).

Write up a lab report with your answers to the questions in bold-face above. Be sure to use complete sentences in your explanations.

Please report any problems with this page to bennett@math.ksu.edu