m=1/sqrt(2*%pi)
function y=f(x)
    y=m*exp(-x*x/2)
endfunction
x=-5+10*rand()
y=m*rand()
while y>f(x)
    x=-5+10*rand()
    y=m*rand()
end
afficher(x)