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

Solution 1 (a).

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

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



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


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

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

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

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

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

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

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

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



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


Form the rational approximation.

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



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


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

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


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

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


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

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


Comparison with the Taylor approximation.  

There were 5 coefficients to determine for the rational approximation, and a Maclaurin polynomial of degree 4 requires 5 coefficients.
Compare with the error in a [Graphics:../Images/RationalApproxMod_gr_31.gif] degree Maclaurin polynomial over the interval  [Graphics:../Images/RationalApproxMod_gr_32.gif].  

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

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

[Graphics:../Images/RationalApproxMod_gr_35.gif]
[Graphics:../Images/RationalApproxMod_gr_36.gif]


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

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


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

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



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


Comparison with the Padé approximation.  

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


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

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



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


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

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


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

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



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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004