CODEWITHSUNDEEP
pythonnumpyscipytypeerrornumerical-integration
Pedro Paxton
Pedro Paxton

Reputation: 23

Function integration returns “only length-1 arrays can be converted to Python scalars”

I'm trying to integrate a function into a given range that shows the flow of particles in a zero-angle (theta) direction as a function of the energy E of the particles. I've tried several ways and got different errors but there are two that persist at the end. My knowledge of Python is limited, I try to learn new ways of doing stuff as I need them but I've been around this function for days with no success.

My function at the moment looks like this:

from numpy import radians, cos, arange
from scipy.integrate import quad

def integral(self):
    theta=0
    E = arange(1, 5000, 1)
    costh = cos(radians(theta))
    a = 18 / (E * costh + 145)
    b = (E + 2.7 / costh)**-2.7
    c = (E + 5) / (E + 5 / costh)
    return a*b*c*1**4


A = quad(integral, 500, 1000)

Applying "quad" to the function like this returns:

TypeError: only length-1 arrays can be converted to Python scalars

If I don't put "self" as argument in the funtion, it returns:

TypeError: integral() takes 0 positional arguments but 1 was given

Has someone an idea on how to bypass this?

Upvotes: 2

Views: 2491

Answers (1)

Adrian W
Adrian W

Reputation: 5026

integral must have one of the signatures described in the scipy documentation. In your case, a function taking a double as argument and returning a double seems appropriate.

self is only used for member functions of a class. It does not make sense outside a class definition. Replace that with a simple name, say x that will be the input argument of the function to integrate.

The function must return a double. Since E is an array and you seem to make calculations with it, then the computed return value is likely to be an array too. It must be a scalar. Fix that and it will work. When you replace your function with something like:

def integral(x):
    return x * x

then it works. This is not what you need, of course, but this is the kind of arguments and return values that are required for quad() to work.

You may want to look here for an example on how quad() is being used.

Some recommendations:

  • use four spaces for indentation (not 5 as in the example)
  • parentheses around the result are not required. Remove for more clarity.

Upvotes: 1

Related Questions