//
// methode de quadrature de Simpson
// application à l'approcimation de pi
//
deff('y=f(x)','y=1/(1+x^2)');
eps=evstr(x_dialog('approcimation souhaitée','1E-5'));
n=int((1/2880*24/eps)^0.25)+1;
a=0;
b=1;
I=0;
for i=1:n
    xi=a+(i-1)*(b-a)/n;
    xi1=a+(i)*(b-a)/n
I=I+1/6*(f(xi)+4*f((xi+xi1)/2)+f(xi1));
end
I=I*(b-a)/n;
format(20)

disp('erreur souhaitee');disp(eps);
disp('valeur de n utilisée:');disp(n);
disp('valeur approchée de pi approchée');disp(4*I);
disp('erreur obtenue');disp(%pi-4*I);
