function I=simps(a,b,h)
n=(b-a)/h+1;
X=linspace(a,b,n);
funkcija;
Y=f(X).^2*pi;
I=(h/3)*(Y(1)+2*sum(Y(3:2:end-1))+4*sum(Y(2:2:end-1))+Y(end));