Exercise 3.  Find the rational approximation [Graphics:Images/RationalApproxMod_gr_216.gif]for [Graphics:Images/RationalApproxMod_gr_217.gif]  over the interval [-1,1].  
3 (a).  Use equally spaced interpolation nodes.

Solution 3 (a).

Set up the formula for  [Graphics:../Images/RationalApproxMod_gr_218.gif].

[Graphics:../Images/RationalApproxMod_gr_219.gif]



[Graphics:../Images/RationalApproxMod_gr_220.gif]


Calculate the equally spaced values for the  [Graphics:../Images/RationalApproxMod_gr_221.gif] interpolation nodes.   

[Graphics:../Images/RationalApproxMod_gr_222.gif]

[Graphics:../Images/RationalApproxMod_gr_223.gif]

Form the  [Graphics:../Images/RationalApproxMod_gr_224.gif] ordinates.  

[Graphics:../Images/RationalApproxMod_gr_225.gif]

[Graphics:../Images/RationalApproxMod_gr_226.gif]

Form the set of  [Graphics:../Images/RationalApproxMod_gr_227.gif] equations to solve and find the solution.

[Graphics:../Images/RationalApproxMod_gr_228.gif]



[Graphics:../Images/RationalApproxMod_gr_229.gif]


Form the rational approximation.

[Graphics:../Images/RationalApproxMod_gr_230.gif]



[Graphics:../Images/RationalApproxMod_gr_231.gif]


Plot graphs of the function and its rational approximation over the interval  [-1,1].  But we will draw the graphs over [-2,2].

[Graphics:../Images/RationalApproxMod_gr_232.gif]


[Graphics:../Images/RationalApproxMod_gr_233.gif]

[Graphics:../Images/RationalApproxMod_gr_234.gif]


Find the error  over the interval  [-1,1].  

[Graphics:../Images/RationalApproxMod_gr_235.gif]


[Graphics:../Images/RationalApproxMod_gr_236.gif]

[Graphics:../Images/RationalApproxMod_gr_237.gif]


Comparison with the Taylor approximation.  

There were 9 coefficients to determine for the rational approximation, and a Maclaurin polynomial of degree 8 requires 9 coefficients.
Compare with the error in a [Graphics:../Images/RationalApproxMod_gr_238.gif] degree Maclaurin polynomial over the interval  [Graphics:../Images/RationalApproxMod_gr_239.gif].  

[Graphics:../Images/RationalApproxMod_gr_240.gif]


[Graphics:../Images/RationalApproxMod_gr_241.gif]

[Graphics:../Images/RationalApproxMod_gr_242.gif]


[Graphics:../Images/RationalApproxMod_gr_243.gif]

[Graphics:../Images/RationalApproxMod_gr_244.gif]

[Graphics:../Images/RationalApproxMod_gr_245.gif]

We can determine how much smaller the error is for the rational approximation.

[Graphics:../Images/RationalApproxMod_gr_246.gif]



[Graphics:../Images/RationalApproxMod_gr_247.gif]


Comparison with the Padé approximation.  

[Graphics:../Images/RationalApproxMod_gr_248.gif]


[Graphics:../Images/RationalApproxMod_gr_249.gif]

[Graphics:../Images/RationalApproxMod_gr_250.gif]



[Graphics:../Images/RationalApproxMod_gr_251.gif]


[Graphics:../Images/RationalApproxMod_gr_252.gif]

[Graphics:../Images/RationalApproxMod_gr_253.gif]


We can determine how much smaller the error is for the rational approximation.

[Graphics:../Images/RationalApproxMod_gr_254.gif]



[Graphics:../Images/RationalApproxMod_gr_255.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004