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