def dblquad(func, a, b, gfun, hfun, args=(), epsabs=1.49e-08, epsrel=1.49e-08):
 def temp_ranges(*args):
  return [gfun(args[0]), hfun(args[0])]
 return nquad(func, [temp_ranges, [a, b]], args=args)

def dblquad ( func , a , b , gfun , hfun , args = ( ) , epsabs = 1.49e-08 , epsrel = 1.49e-08 ) : def temp_ranges ( * args ) : return [ gfun ( args [ 0 ] ) , hfun ( args [ 0 ] ) ] return nquad ( func , [ temp_ranges , [ a , b ] ] , args = args )
compute a double integral .