[Scipy-tickets] [SciPy] #1763: scipy.special.hyp1f1(0.5, 1.5, -1000) fails

SciPy Trac scipy-tickets@scipy....
Sat Nov 3 15:12:45 CDT 2012


#1763: scipy.special.hyp1f1(0.5, 1.5, -1000) fails
------------------------------+---------------------------------------------
 Reporter:  unknown (google)  |       Owner:  pv         
     Type:  defect            |      Status:  new        
 Priority:  normal            |   Milestone:  Unscheduled
Component:  scipy.special     |     Version:  0.10.1     
 Keywords:  hyp1f1            |  
------------------------------+---------------------------------------------
 >>> scipy.special.hyp1f1(0.5, 1.5, -1000)
 nan

 According to wolfram alpha this is
 (1/20)*sqrt(pi/10)*erf(10/sqrt(10))
 which is essentially
 (1/20)*sqrt(pi/10)
 which is about
 0.028025

 This is correctly computed by mpmath:
 >>> mpmath.hyp1f1(0.5, 1.5, -1000)
 mpf('0.028024956081989643')

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1763>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list