%
%   Tema d'esame CALCOLO NUMERICO - MILANO 7 maggio 2013 - Esercizio 2
%
clear all

x(1)=-1/sqrt(3); x(2)=-x(1);
a=0; b=4;
a1=(b-a)/2; a0=(b+a)/2;
effe=inline('exp(-t.^2)','t');
z=a1*x+a0;
w(1:2)=a1;
I_erf=sqrt(pi)/2*erf(4);
I_s=sum(w.*effe(z));
ac=[a a0 b];
I_c=0;
for i=1:2
    a=ac(i); b=ac(i+1);
    a1=(b-a)/2; a0=(b+a)/2;
    z=a1*x+a0;
    w(1:2)=a1;
    I_c=I_c+sum(w.*effe(z));
end
err_s=abs(I_erf-I_s);
err_c=abs(I_erf-I_c);
fprintf('I_erf=%8.5f I(0,4)=%8.5f I(0,2)+I(2,4)=%8.5f \n',I_erf,I_s,I_c)
fprintf('err_(0,4)=%8.5f err_(0,2,4)=%8.5f \n',err_s,err_c);

%I_erf= 0.88623 I(0,4)= 0.97894 I(0,2)+I(2,4)= 0.92231 
%err_(0,4)= 0.09271 err_(0,2,4)= 0.03609 