Problem 179. Monte-Carlo integration

Solution 1590931

Submitted on 24 Jul 2018 by J. S. Kowontan

This solution is locked. To view this solution, you need to provide a solution of the same size or smaller.

Test	Status	Code Input and Output
1	Pass	fun = @(x) x^2; assert(abs((dquad(1,fun) - 1/3)*3)<0.01)
2	Pass	fun = @(x) x(1)x(2); assert(abs((dquad(2,fun) - 0.25)4)<0.01)
3	Pass	fun = @(x) 1+sin(x(1)); assert(abs((dquad(50,fun) - 1.45969769)/1.45969769)<0.01)
4	Pass	d = randi([5 10],1) r = rand0.8 fun = @(x) 2^d(norm(x)>r); dball = exp(d/2log(pi)+dlog(r)-gammaln(d/2+1)); assert(abs((dquad(d,fun) - 2^d+dball)/(2^d-dball))<0.01) d = 6 r = 0.4300

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!