[SciPy-user] Bug in scipy.special.ivp

Takanobu Amano amano at eps.s.u-tokyo.ac.jp
Tue Jun 13 00:50:11 CDT 2006


Hello,

I found a small bug in scipy.special.ivp (scipy/Lib/special/basic.py)
this patch is for the latest svn snapshot (revision 1882).

Takanobu Amano

--- basic.py.orig       2006-06-13 14:29:05.766772544 +0900
+++ basic.py    2006-06-13 14:23:58.534420748 +0900
@@ -168,7 +168,7 @@
     if n == 0:
         return iv(v,z)
     else:
-        return (ivp(v-1,z,n-1) - ivp(v+1,z,n-1))/2.0
+        return (ivp(v-1,z,n-1) + ivp(v+1,z,n-1))/2.0

 def h1vp(v,z,n=1):
     """Return the nth derivative of H1v(z) with respect to z.



More information about the SciPy-user mailing list