function [I] = mc1(f,a,b,N)

I = (b-a) * mean(f(a + (b-a).*rand(N,1)));


end