function outv=aq(a,m,b,fa,fm,fb,tol)
%
global count;
global x0;
h=(b-a)/4;
x=[a+h,b-h];
y=myf(x);
S1=2*h*(fa+4*fm+fb)/3;
S2=h*(fa+4*y(1)+fm)/3;
S3=h*(fm+4*y(2)+fb)/3;
if (abs(S2+S3-S1)>15*tol)
    x0=[x0,m];
    S2=AQ(a,x(1),m,fa,y(1),fm,tol/2);
    S3=AQ(m,x(2),b,fm,y(2),fb,tol/2);
end
outv=S2+S3;