function V=zapremina(a,b,tol)
h=(b-a)/2;
I1=simps(a,b,h);
h=h/2;
I2=simps(a,b,h);
while(Runge(I1,I2)>tol)
    I1=I2;
    h=h/2;
    I2=simps(a,b,h);
end
V=I2;
