(ph#@SSKJrJrJrJrJrJr SSKJr S/r SSjr g))zerosasarrayeyepoly1dhstackr_)linalgpadeNc[U5nUc"[U5S- U- nUS:a [S5eUS:a [S5eX-nU[U5S- :a [S5eUSUS-n[US-US-URS9n[ US-U4URS9n[ SUS-5HnUSUSSS2*XVSU24'M [ US-US-5HnXU- USSS2*XVSS24'M [XE45n[R"Xp5nUSUS-n [S XS-S4n [U SSS25[U SSS254$) a5 Return Pade approximation to a polynomial as the ratio of two polynomials. Parameters ---------- an : (N,) array_like Taylor series coefficients. m : int The order of the returned approximating polynomial `q`. n : int, optional The order of the returned approximating polynomial `p`. By default, the order is ``len(an)-1-m``. Returns ------- p, q : Polynomial class The Pade approximation of the polynomial defined by `an` is ``p(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 Nrz.Order of q must be smaller than len(an)-1.z&Order of p must be greater than 0.z0Order of q+p must be smaller than len(an).)dtypeg?) rlen ValueErrorrr rrangerr solverr) anmnNAkjBkjrowCpqpqs J/var/www/html/venv/lib/python3.13/site-packages/scipy/interpolate/_pade.pyr r s{N By GaK!O q5MN N1uABB A3r719}KLL DQqSB ac1Q3bhh 'C 1a )CQ!}Tc(DbD))H QqS!A#a%}dd++E zA a B 4AaCA 3Q34=A !DbD'?F1TrT7O ++)N) numpyrrrrrrscipyr __all__r rrr$s99 (<,r