monte_carlo_integration_py_script.py

import random 

a = 0
b = 1
N = 1000

def f( x ):
    return 2 * x**5

def uniform():
    return random.uniform( a , b )

def compute():
    M = f( uniform() )
    for k in range( N ):
        M += (( f( uniform() ) - M ) / (k + 1) )
    return M * ( b - a )

integral = compute()
print( integral )