A=zeros(7,3);
calka=c(b)-c(a);
#for m=1:3
for k=1:7
  A(k,1)=asq(a,b,f,10^(-k));
  x(1,k)=log(A(2)); #liczba przedziałów
  x(2,k)=-log(abs(A(1)-calka)); #różnica dla metody asq
  #c(b)-c(a) - wartosc całki policzona analitycznie, trzeba zdefiniować 
                       #razem z a,b,f itd
  x(3,k)=abs( nsq(a,b,f,A(2))-calka);   #różnica dla metody NSQ                  
 endfor 
#endfor
  plot(A(1,:), A(2,:), ";ASQ", A(1,:), A(3,:), "'NSQ;")