n=100000;
s=0;

for i=1:n
x=rand(1,1);
y=rand(1,1) * sqrt(x);
s = s + x^2 * y;
end

I = s/n;