function IN2 = slozeneSimpsonovoPravidlo(f,a,b,N)

H = (b-a)/N;
ksik = a:H:b;
stredy = 0.5*(ksik(1:N)+ksik(2:(N+1)));

IN2 = H * 0.5*( (4/3)*sum(feval(f,stredy)) + ...
    (2/3)*sum(feval(f,ksik)) - ...
    (1/3)*(feval(f,ksik(1))+feval(f,ksik(N+1))) );