[SciPy-user] Using scipy specfunc in integration
Fri Jun 13 13:57:40 CDT 2008
On Fri, Jun 13, 2008 at 09:21, Neal Becker <firstname.lastname@example.org> wrote:
> Any ideas on this?
> from scipy.special import erf
> from math import exp, tan
> def cot(x):
> return 1/tan(x)
> N = 8
> esnodB = 10
> Rd = 10**(.1 * esnodB)
> def F (y):
> return exp (-(y**2)) * erf (y * cot (pi/N))
> Pe = float(N-1)/float(N) - 0.5 * erf (sqrt (Rd * sin (pi/N))) - 1/(sqrt(pi))
> * quadrature(F, 0, sqrt(Rd) * sin (pi/N))
> TypeError: only length-1 arrays can be converted to Python scalars
> It seems to be complaining about erf.
> Out: <ufunc 'erf'>
I think it's complaining about exp(), actually. quadrature() is going
to pass arrays to F(), not scalars.
"I have come to believe that the whole world is an enigma, a harmless
enigma that is made terrible by our own mad attempt to interpret it as
though it had an underlying truth."
-- Umberto Eco
More information about the SciPy-user