Example 5.  Find the rational approximation [Graphics:Images/RationalApproxMod_gr_375.gif]for [Graphics:Images/RationalApproxMod_gr_376.gif]  over the interval [-1,1].  
5 (a).  Use equally spaced interpolation nodes.

Solution 5 (a).

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

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



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


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

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

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

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

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

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

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

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



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


Form the rational approximation.

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



[Graphics:../Images/RationalApproxMod_gr_390.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_391.gif]


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

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


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

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


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

[Graphics:../Images/RationalApproxMod_gr_396.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_397.gif] degree Maclaurin polynomial over the interval  [Graphics:../Images/RationalApproxMod_gr_398.gif].  

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


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

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


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

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

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


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

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



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


Comparison with the Padé approximation.  

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


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

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



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


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

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


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

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



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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004