k=9e9; q=1e-9; r0=0.1; u0=k*q/r0; [X,Y,Z]=sphere(20); r0=0.1; x=r0*X(:)'; y=r0*Y(:)'; z=r0*Z(:)'; x=[x;zeros(size(x))]; y=[y;zeros(size(y))]; z=[z;zeros(size(z))]; plot3(x,y,z);hold on; axis square; shading interp; 改写一下这个代码，不改变其功能，并添加注释

% 定义常量
k = 9e9; % 真空中的电磁力常数
q = 1e-9; % 电荷量
r0 = 0.1; % 初始半径

% 计算电势能
u0 = k * q / r0; % 电势能公式

% 生成球体点坐标
[X, Y, Z] = sphere(20); % 生成球体的点坐标
r0 = 0.1; % 球体半径
x = r0 * X(:)'; % 将球体点坐标按行展开并乘以半径
y = r0 * Y(:)';
z = r0 * Z(:)';

% 添加底部平面点坐标
x = [x; zeros(size(x))]; % 在每个点坐标后添加一个零元素，生成底部平面
y = [y; zeros(size(y))];
z = [z; zeros(size(z))];

% 绘制三维图形
plot3(x, y, z); % 绘制球体
hold on; % 添加保持绘图状态，方便后续操作
axis square; % 设置坐标轴比例为1:1:1，使球体显示为正圆
shading interp; % 设置颜色插值方式为平滑插值

% 添加注释
title('Electric Potential Energy Distribution'); % 添加图标题
xlabel('X-axis'); % 添加X轴标签
ylabel('Y-axis'); % 添加Y轴标签
zlabel('Z-axis'); % 添加Z轴标签
legend('Sphere', 'Location', 'NorthEast'); % 添加图例，放置在右上角