// ---------------------------------------------------------
// simpsonf3.sce
// ---------------------------------------------------------
clear; 
exec integ.sci;
a = 0 ; b = 22 ; [exactinteg,err1] = intg(0,100,fex1x2,0,1e-12); 
i=0;
for N=20:10:80;
  i=i+1;
  nbN(i)=N;
  erreur1s(i) = abs(exactinteg - integsymp(fex1x2,a,b,N));
end
//
[nbN, erreur1s]