python中的二重积分

Question

我想在python中做一个二重积分。在http://docs.scipy.org/doc/scipy/reference/tutorial/integrate.html的帮助下，我做了以下事情

import numpy as np
import matplotlib.pyplot as plt
from scipy import integrate


def func(x,a):
    return a*x*x

def integrl(a):
    return integrate.quad(func,0,10,args=(a))

print integrate.quad(lambda x: integrl(x),0,1)

但是我得到了以下错误

^{pr2}$

我该怎么办？在

Answer 1

将integrl的返回改为：

return integrate.quad(func,0,10,args=(a))[0]

在集成.quad默认情况下，函数返回一个包含积分值和绝对误差估计值的元组：scipy doc

Answer 2

错误在integrl函数中。我应该指定与哪个参数集成。正确的代码应该是

import numpy as np
import matplotlib.pyplot as plt
from scipy import integrate


def func(x,a):
   return a*x*x

def integrl(a):
   return integrate.quad(func,0,10,args=(a))[0]

print integrate.quad(lambda x: integrl(x),0,1)