import sys
sys.version

from scipy.integrate import quad

import numpy as np


def f(x):
    return 2 * x**2 - 5 * x + 1


def g(x):
    a = np.sin(f(x - 1))
    b = f(x**2)
    return a * b

quad(g, -5, 5)

xx = np.linspace(-5, 5, 1000)
yy = g(xx)

import matplotlib.pyplot as plt
plt.plot(xx, yy)
plt.grid()

def f(a, b):
    def wrapped(r):
        return a * r + b

    return wrapped

f2 = f(2, 5)
f2(10)

def f(sigma, epsilon):
    def U(r):
        return 4*epsilon*(((sigma/r)**12)-((sigma/r)**6))
    return U

def F(r,k,T,U):
    return (r**2)*(1-(np.exp((-U(r))/(k*T))))

F(3, 1.2, 2.2, f(3.3, 4.4))

from scipy.integrate import quad
quad(lambda x : F(x,  1.2, 2.2, f(3.3, 4.4)), 0, 1000)