Modeling4

import random

a = 0
b = 1
N = [10, 20, 50, 100, 1000, 10000]

# 1/3 результат


def f(x_):
    return x_ * x_


def get_random():
    return a + random.uniform(0, 1) * (b - a)


for n in N:
    m = 0
    for _ in range(0, n):
        x = get_random()
        fx = f(x)
        y = get_random()
        if y < fx:
            m += 1
    print("N = {0}, результат: {1}".format(n, m / n))