使用蒙特卡洛积分法计算多重积分，采用MATLAB编程

蒙特卡洛求解积分时，与积分重数无关,下面是蒙特卡洛期望法：

例子：所需求解积分

clc;clear all;close all;
fh3 = @(x)sqrt(x(:,1).*x(:,2)) .* log(x(:,3)) + sin(x(:,4)./x(:,2));
n = 1e8;%采样个数
m = 1/n;

%生成采样点
x(:,1) = unifrnd(1,2,1,n);
x(:,2) = unifrnd(1,6,1,n);
x(:,3) = unifrnd(1,24,1,n);
x(:,4) = unifrnd(2,98,1,n);

ind3 = x(:,2) >= x(:,1) & x(:,2) <= 3*x(:,1) &...
    x(:,3) >= x(:,1).*x(:,2) & x(:,3) <= 2*x(:,1).*x(:,2)&...
    x(:,4) >= (x(:,1) + x(:,1).*x(:,3)) & x(:,4) <= (x(:,1) + 2*x(:,1).*x(:,3));%索引
fh3sum = fh3(x(ind3,:));
i3 = (2-1)*(6-1)*(24-1)*(98-2)*m*sum(fh3sum)