v=0.4*rand();

while v>sqrt(4-u^2)/(2*%pi),

u=4*rand()-2;

v=0.4*rand();

end

disp(u)