s=0
n=input("nombre de points = ")
for i =1:n
    x=rand()
    y=rand()
    if y<=exp(-x^2/2) then s=s+1
    end
end
afficher(s/n)