%% skrypt calkuje dwie funkcjie exp(x^2) i tan(x^(1/2))
f1=inline('exp(x.*x)');
f2=inline('tan(sqrt(x))');
q11=zeros(15,1);
q12=zeros(15,1);
q21=zeros(15,1);
q22=zeros(15,1);

for indx=7:21
q11(indx)=quad(f1,0,5,10^(-indx));
q12(indx)=quadl(f1,0,5,10^(-indx));

q21(indx)=quad(f2,0,pi/6,10^(-indx));
q22(indx)=quadl(f2,0,pi/6,10^(-indx));
end

hold on;
subplot(2,2,1);
plot(-indx,q11,'-r');
title('exp(x^2) newton-cotes')
subplot(2,2,2);
plot(-indx,q12,'-r');
title('exp(x^2) simpson')
subplot(2,2,3);
plot(-indx,q22,'-r');
title('tan(sqrt(x)) newton-cotes')
subplot(2,2,4);
plot(-indx,q22,'-r');
title('tan(sqrt(x)) simpson')

hold off;