I have to implement Monte Carlo integration to estimate:
theta = S h(x)f(x)dx (S is the integration symbol)
theta = S h(x)f(x)dx
where f(x) = k(x-1)^3
f(x) = k(x-1)^3