Return Pade approximation to a polynomial as the ratio of two polynomials.
Parameters : | an : (N,) array_like
m : int
|
---|---|
Returns : | p, q : Polynomial class
|
Examples
>>> from scipy import misc
>>> e_exp = [1.0, 1.0, 1.0/2.0, 1.0/6.0, 1.0/24.0, 1.0/120.0]
>>> p, q = misc.pade(e_exp, 2)
>>> e_exp.reverse()
>>> e_poly = np.poly1d(e_exp)
Compare e_poly(x) and the pade approximation p(x)/q(x)
>>> e_poly(1)
2.7166666666666668
>>> p(1)/q(1)
2.7179487179487181