用matlab实现：考虑一个带电系统，用位于坐标原点，半径为R的球BR来表示缅甸和密度为σ的胶粒

以下是用MATLAB实现的代码：

% 定义球半径和胶粒密度
R = 1;
sigma = 2;

% 生成三维网格
[x,y,z] = meshgrid(-R:0.1:R,-R:0.1:R,-R:0.1:R);

% 计算每个点到球心的距离
d = sqrt(x.^2+y.^2+z.^2);

% 判断点是否在球内
is_inside = (d <= R);

% 计算电势能
U = 0;
if is_inside
    U = (sigma * 4/3 * pi * R^3) / d;
end

% 显示三维图像
isosurface(x,y,z,U)
xlabel('x')
ylabel('y')
zlabel('z')

运行结果如下所示：