
a=-1; b=1;

f=inline('exp(x)');
I_esatto=exp(1)-exp(-1);

N_vett=2:8;
ajac=0; bjac=0;


for index=1:length(N_vett)
    N=N_vett(index);
   [I_jac(index,1),x_jac,w_jac]=integrazione_gauss_jacobi(N,ajac,bjac,a,b,f);
   err_assoluto(index,1)=abs(I_jac(index,1)-I_esatto);
   fprintf('\n \t [N]: %3.0f [ERR.ASS.]: %2.2e',N,err_assoluto(index,1));
   
end


