_Mh#@ddlmZmZmZmZmZmZddlmZdgZ ddZ dS))zerosasarrayeyepoly1dhstackr_)linalgpadeNcJt|}|*t|dz |z }|dkrtd|dkrtd||z}|t|dz krtd|d|dz}t|dz|dz|j}t |dz|f|j}t d|dzD]}|d|ddd ||d|f<t |dz|dzD] }|||z |ddd ||ddf<!t||f}tj ||}|d|dz} td ||dzdf} t| dddt| dddfS) a 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 W/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/scipy/interpolate/_pade.pyr r sN By GGaK!O q55MNN N1uuABBB AA3r7719}}KLLL DQqSDB ac1Q3bh ' ' 'C 1a ) ) )CQ!}}**TcT(DDbD))CH QqS!A#,,#a%)}ddd++CE SzA a  B 4AaC4A 31Q344=A !DDbD'??F1TTrT7OO ++)N) numpyrrrrrrscipyr __all__r rrr$sl9999999999999999 (<,<,<,<,<,<,r