| Server IP : 35.80.110.71 / Your IP : 216.73.216.221 Web Server : Apache/2.4.58 (Ubuntu) System : Linux ip-172-31-21-44 6.17.0-1019-aws #19~24.04.1-Ubuntu SMP Tue Jun 23 18:53:06 UTC 2026 x86_64 User : ubuntu ( 1000) PHP Version : 8.3.31 Disable Function : NONE MySQL : OFF | cURL : ON | WGET : ON | Perl : ON | Python : OFF | Sudo : ON | Pkexec : OFF Directory : /lib/python3/dist-packages/scipy/interpolate/__pycache__/ |
Upload File : |
�
x�f# � �> � d dl mZmZmZmZmZmZ d dlmZ dgZ dd�Z
y)� )�zeros�asarray�eye�poly1d�hstack�r_)�linalg�padeNc � � t | � } |�!t | � dz
|z
}|dk r t d� �|dk r t d� �||z }|t | � dz
kD r t d� �| d|dz } t |dz |dz | j �� }t |dz |f| j �� }t
d|dz � D ] }| d| ddd� ||d|�f< � t
|dz |dz � D ] }| ||z
| ddd� ||dd�f< � t ||f� }t j || � }|d|dz } t d ||dz d f }
t | ddd� � t |
ddd� � fS )
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
N� r z.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���g �?) r �len�
ValueErrorr r
r �ranger r �solver r ) �an�m�n�N�Akj�Bkj�row�C�pq�p�qs �9/usr/lib/python3/dist-packages/scipy/interpolate/_pade.pyr
r
s� � �N
���B��y� ��G�a�K�!�O�� �q�5� � M�N� N��1�u��A�B�B� �A��A��3�r�7�1�9�}��K�L�L� �D�Q�q�S��B�
�a��c�1�Q�3�b�h�h�
'�C�
��1��a�����
)�C��Q��!��}� *�� �T�c�(�D�b�D�)�)��C����H�
�*��Q�q�S�!�A�#�� ,���#�a�%��}�d��d�+�+��C��E�
�,���S�z��A� � � �a�� �B�
�4�A�a�C��A�
�3��1�Q�3�4��=��A� �!�D�b�D�'�?�F�1�T�r�T�7�O� +�+� )N) �numpyr r r r r r �scipyr �__all__r
� r r �<module>r$ s �� 9� 9� � �(��<,r