[SciPy-user] Prob with romberg integration?
Neal Becker
ndbecker2 at gmail.com
Thu Oct 12 14:49:17 EDT 2006
This is scipy-0.5.1:
---------------------------
from scipy.integrate import romberg
from math import sqrt, log, exp
def F (x):
return x
def G (v):
return exp (v) * F (exp (v))
gain = 1/sqrt (2*romberg (G, log (1e-7), log (0.5)))
------------------------
Traceback (most recent call last):
File "<stdin>", line 1, in ?
File "/usr/tmp/python-MSHhoK.py", line 11, in ?
gain = 1/sqrt (2*romberg (G, log (1e-7), log (0.5)))
File "/usr/lib64/python2.4/site-packages/scipy/integrate/quadrature.py", line 377, in romberg
ordsum = ordsum + _difftrap(vfunc, interval, n)
File "/usr/lib64/python2.4/site-packages/scipy/integrate/quadrature.py", line 326, in _difftrap
s = sum(function(points),0)
File "/usr/lib64/python2.4/site-packages/scipy/integrate/quadrature.py", line 54, in vfunc
output = empty((n,), dtype=y0.dtype)
AttributeError: 'float' object has no attribute 'dtype'
More information about the SciPy-User
mailing list