[SciPy-dev] Problems with quadrature.py

[Nils Wagner]

    m,abserr = integrate.quadrature(f,0.0,1.0,args=p)
  File "/usr/lib/python2.4/site-packages/scipy/integrate/quadrature.py",
line 74, in quadrature
    newval = fixed_quad(vec_func,a,b,(func,)+args,n)[0]
TypeError: unsupported operand type(s) for +: 'function' and 'int'

Is this a bug ?

Nils


p[0] = 1
p[1] = 1

def f(x,p):
   return x**2*(1.0-0.5*x)**(p[0]+p[1]+2)