How to evaluate single integrals of multivariate functions with Python's scipy.integrate.quad?

Question

There is a function I am trying to integrate in Python using scipy.integrate.quad. This particular function takes two arguments. There is only one argument I

Answer 1

Found the answer in the scipy documentation.

You can do the following:

from scipy import integrate as integrate
def f(x,a):  #a is a parameter, x is the variable I want to integrate over
    return a*x

result = integrate.quad(f,0,1,args=(1,))

The args=(1,) argument in the quad method will make a=1 for the integral evalution.

This can also be carried to functions with more than two variables:

from scipy import integrate as integrate
def f(x,a,b,c):  #a is a parameter, x is the variable I want to integrate over
    return a*x + b + c

result = integrate.quad(f,0,1,args=(1,2,3))

This will make a=1, b=2, c=3 for the integral evaluation.

The important thing to remember for the function you want to integrate this way is to make the variable you want to integrate over the first argument to the function.

Answer 2

Use the args argument (see the scipy documentation):

result = integrate.quad(f,0,1, args=(a,))

The comma in args=(a,) is required because a tuple must be passed.