performing mathematical integration

import scipy.integrate

import math

import numpy as np

from numpy import exp

f= lambda x:math.sin(x)

i = scipy.integrate.quad(f, 0, np.pi/2)

print ("integral value ",i[0])

print ("error ",i[1])

import scipy.integrate

import math

import numpy as np

f = lambda x, y : 16*x*y

g = lambda x : 0

h = lambda y : math.sqrt(4*y**2)

i = scipy.integrate.dblquad(f, 0, 1, g, h)

print ("integral value ",i[0])

print ("error ",i[1])