function y = f(x);y=1/(1+x^2)
endfunction
n=input("n = ")
T=0;S=0
for k=0:n-1
    a=k/n;b=(k+1)/n
    T=T+(b-a)*(f(a)+f(b))/2
    S=S+(b-a)/6*(f(a)+4*f((a+b)/2)+f(b))
    ecartT=abs(T-%pi/4)
    ecartS=abs(S-%pi/4)
end
afficher ("méthode des trapèzes : " + string(T)+" ;écart à pi/4 :" + string(ecartT))
afficher ("méthode de Simpson : " + string(S)+" ;écart à pi/4 :" + string(ecartS))