[求助] 隐函数曲面求积分

ClearAll["Global`*"];n=125000;SeedRandom[0];
f[x_]:=(Cos/@x//Total)+51/100*Dot[Cos/@RotateRight@x,Cos/@x]+1
Count[Flatten@(f/@RandomReal[{0,2Pi},{n,3}]),u_/;u<=0]/n*(2.0 Pi)^3

2楼: Originally posted by Mr__Right at 2017-05-23 21:56:56
Monte Carlo
计算近似值

ClearAll;n=125000;SeedRandom;
f:=(Cos/@x//Total)+51/100*Dot+1
Count/n*(2.0 Pi)^3

67.3446...