A simpler way of doing this is to use the nquad
function instead of dblquad
. Example code:
from scipy.integrate import nquad
options={'limit':100}
integral=nquad(func,[[xmin,xmax],[ymin,ymax]],
args=(other_arg,),opts=[options,options])
Note that several of the arguments are lists. The elements of these lists apply to each of the coordinates in order. See the documentation for nquad here.