多谢大佬点开此帖，
代码目的：在一个3维曲面（如图）的周围取3维的随机数，随机数以曲面为中心呈正态分布
所遇问题：不知如何操作
可能用到的函数：MultinormalDistribution、RandomVariate
代码：
ClearAll["`*"]
{eval, evec} = With[{ v1 = 5.922,
v2 = 1.25,
v3 = 0.1974,
m0 = -0.0047,
mp = 15,
mz = 0.15,
\[Mu]B = 5.788*10.^-5,
g
bp
d = 33.75,
gbmd = 11.25},
Function[{bz , kx, ky , kz }, #]] & /@
ToRadicals@
Eigensystem[{{m0 + (kx^2 + ky^2) mp + kz^2 mz - (bz gbpd \[Mu]B)/
2, 0, kx v1 - I ky v2, kz v3}, {0,
m0 + (kx^2 + ky^2) mp + kz^2 mz + (bz gbpd \[Mu]B)/2,
kz v3, -kx v1 - I ky v2}, {kx v1 + I ky v2,
kz v3, -m0 - (kx^2 + ky^2) mp - kz^2 mz - (bz gbmd \[Mu]B)/2,
0}, {kz v3, -kx v1 + I ky v2,
0, -m0 - (kx^2 + ky^2) mp - kz^2 mz + (bz gbmd \[Mu]B)/2}}];
ev[kx_, ky_, kz_] := Evaluate[(eval @@ {5, kx, ky, kz} )[[1]]]）（*定义曲面的函数*）
（*曲面方程为ev @@ {kx, ky, kz} ==0.15*）

什么叫“以曲面为中心的正态分布”？究竟什么变量遵循正态分布？