function q = romberg(fx,a,b,n)
h=(b-a)/n;
q=0
x=a
for i=1:n
    x=x+(n-1)*h;
    x0=x;
    x1=x+h/8;
    x2=x+2*h/8;
    x3=x+3*h/8;
    x4=x+4*h/8;
    x5=x+5*h/8;
    x6=x+6*h/8;
    x7=x+7*h/8;
    x8=x+8*h/8;
    zz=h/5670(441*fx(x0)+2048*fx(x1)+736*fx(x2)+2048*fx(x3)+2036*fx(x4)+2048*(x5)+736*(fx(x6)+2048*fx(x7)+441*fx(x8));
    q=q+zz;
end
q