function q = my_quad(f, a, b, h)

N = ceil((b-a)/h);
q = 0;

for i = 1:N     		% loop över samtliga delintervall
   xi = a + i*h;                % delintervallets högra ändpunkt
   q = q + feval(f, xi)*h;      % ``right-hand rectangle rule''
end
