Python scipy.integrate.quand import scipy.integrate as integrate import numpy as np if __name__ == '__main__': res = integrate.quad(lambda x: 0.02 * x * np.exp(-0.02 * x), 0, float('inf')) # result and error print(res)