+
    nDj	                     v    ^ RI t ^ RIt^ RIHtHtHtHtHtHt ^ RI	H
t
 ^ RIHt R.t]! RR7      R	R l4       tR# )
    N)zerosasarrayeyepoly1dhstackr_)linalg)xp_capabilitiespadeT)out_of_scopec                   \         P                  P                  \        4      3pRp\        P
                  ! V\        VR7       \        V 4      p Vf,   \        V 4      ^,
          V,
          pV^ 8  d   \        R4      hV^ 8  d   \        R4      hW,           pV\        V 4      ^,
          8  d   \        R4      hV RV^,            p \        V^,           V^,           V P                  R7      p\        V^,           V3V P                  R7      p\        ^V^,           4       F  pV RV RRR
1,          ) WxRV13&   K  	  \        V^,           V^,           4       F  pWV,
          V RRR
1,          ) WxR3&   K  	  \        Wg34      p	\        P                   ! W4      p
V
RV^,            p\"        R	W^,           R 3,          p\%        VRRR
1,          4      \%        VRRR
1,          4      3# )a  
Return Pade approximation to a polynomial as the ratio of two polynomials.

.. deprecated:: 1.18.0
    This function is deprecated and will be removed in SciPy 1.20.0. Use
    `mpmath.pade` instead.

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

zT`pade` is deprecated and will be removed in SciPy 1.20.0. Use `mpmath.pade` instead.)skip_file_prefixesNz.Order of q <m> must be smaller than len(an)-1.z&Order of p <n> must be greater than 0.z0Order of q+p <m+n> must be smaller than len(an).)dtype:NNNg      ?)ospathdirname__file__warningswarnDeprecationWarningr   len
ValueErrorr   r   r   ranger   r	   solver   r   )anmn_warn_skipsmsgNAkjBkjrowCpqpqs   &&&          Q/data/cameron/venvs/s3viz/lib/python3.14/site-packages/scipy/interpolate/_pade.pyr   r   
   s   X 77??8,.K$CMM#)kJ	ByGaK!Oq5MNN1uABB	A3r719}KLL	DQqSB
ac1Q3bhh
'C
1a
)CQ!}Tc(DbD))H QqS!A#a%}dd++E
 zA	a	B
4AaCA
3Q34=A!DbD'?F1TrT7O++    )N)r   r   numpyr   r   r   r   r   r   scipyr	   scipy._lib._array_apir
   __all__r    r*   r)   <module>r0      s7     	 9 9  1(d#D, $D,r*   