[SciPy-dev] Pade approximants

Nils Wagner nwagner at mecha.uni-stuttgart.de
Wed Feb 23 12:45:49 CST 2005


On Wed, 23 Feb 2005 09:43:03 -0800
  Robert Kern <rkern at ucsd.edu> wrote:
> Nils Wagner wrote:
>> Hi all,
>> 
>> Has someone written a wrapper for the FORTRAN subroutine 
>>pade
>> which is available through
>> http://lib-www.lanl.gov/numerical/bookfpdf/f5-12.pdf
>> 
>> It would be nice to have a function called pade in scipy 
>>:-) .
>> 
>> Any comments or suggestions ?
> 
> Anything wrong with the one that's already there?
> 
Oops, I have ignored that.
Anyway, it doesn't provide the user with a residual
as in the case of Numerical Recipes routine.
What is the meaning of poly1d
x = pade(cof,2)
>>> cof
array([ 1.        ,  1.        ,  0.5       ,  0.16666667, 
 0.04166667])
>>> x
(poly1d([ 0.08333333,  0.5       ,  1.        ]), poly1d([ 
0.08333333, -0.5       ,  1.        ]))


Nils
  
> -- 
> Robert Kern
> rkern at ucsd.edu
> 
> "In the fields of hell where the grass grows high
>  Are the graves of dreams allowed to die."
>   -- Richard Harter
> 
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-dev

  




More information about the Scipy-dev mailing list