n=input("donner le nombre de points : ");
clf;
for k=1:n
  x=rand();y=rand();
  if y>sqrt(1-x^2) then 
    plot(x,y,"+r")
  end 
end 
