clear all;
global count;
global x0;
count=0;
a=0.05;
b=2;

x0=[a,b];
m=(a+b)/2;
saq=aq(a,m,b,myf(a),myf(m),myf(b),10^-4)
count1=count
x1=[a:.01:b];
y1=myf(x1);
y0=zeros(size(x0));
plot(x0,y0,'bo',x1,y1,'k-')

count=0;
[r,A]=romberg(a,b,10^-4,15)
count2=count

count=0;
[Tm,Ta]=trapid(a,b,10^-4,20)
count4=count

count=0;
[Sm,Sa]=SimpID(a,b,10^-4,20)
count5=count

count=0;
gq=gaussquad(a,b,4,10^-4,15)
count3=count