function y=f(x);y=1/sqrt(1+x^2);endfunction
n=input("n = ")
u=0;v=0
for k=0:n-1
    v=v+1/n*f(k/n)
    u=u+1/n*f((k+1)/n)
end
afficher ([u,v])