pade — SciPy v1.15.3 Manual (original) (raw)

scipy.interpolate.

scipy.interpolate.pade(an, m, n=None)[source]#

Return Pade approximation to a polynomial as the ratio of two polynomials.

Parameters:

an(N,) array_like

Taylor series coefficients.

mint

The order of the returned approximating polynomial q.

nint, optional

The order of the returned approximating polynomial p. By default, the order is len(an)-1-m.

Returns:

p, qPolynomial class

The Pade approximation of the polynomial defined by an isp(x)/q(x).

Examples

import numpy as np from scipy.interpolate import pade e_exp = [1.0, 1.0, 1.0/2.0, 1.0/6.0, 1.0/24.0, 1.0/120.0] p, q = 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