﻿from math import *
from random import *

def M_C(N):
    L=0
    for k in range(1,N+1):
        x=2*random()
        y=random()
        if y < exp(-x**2):
            L=L+1
    S=2*L/N
    return S



