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