# -*- coding: utf-8 -*-

from scipy import integrate

gaussian = lambda x: ((6.73*x)+4.3025*10**(-7))/((2.316*10**(-11))*x)
result = integrate.romberg(gaussian, 0.61*(10**-6), 1.22*(10**-6), show=True)


#import numpy as np
#from scipy import integrate
#
#gaussian = lambda x: np.sin(x)
#result = integrate.romberg(gaussian, 0, np.pi, show=True)