clear all
close all
N=100; % number of samples
r=rand(N,1); % uniform random numbers
f=r.*r;
I=mean(f)
plot(r,f,'b*')
hold on
plot(r,zeros(N,1),'b*')

