现需用MATLAB构建模型来求解一任意形状的、质量分布均匀的球体绕一定轴转动的转动惯量。 模型分为以下2部分: 1.模拟一任意形状的的质量分布均匀的球体（可以自己输入半径）； 2.求这个球体绕一定轴转动的转动惯量。 最后显示模拟的图形和求得的转动惯量。请用matlab编程，给出完整代码

以下是 MATLAB 的代码实现：

% 输入球体半径和分辨率
r = input("请输入球体半径：");
n = input("请输入分辨率：");

% 生成球体的网格
[x, y, z] = sphere(n);

% 缩放球体尺寸并平移至球心
x = x * r;
y = y * r;
z = z * r;
x = x + r;
y = y + r;
z = z + r;

% 计算球体的质量分布
rho = ones(size(x));

% 计算球体的转动惯量
Ixx = sum(sum(rho .* (y.^2 + z.^2)));
Iyy = sum(sum(rho .* (x.^2 + z.^2)));
Izz = sum(sum(rho .* (x.^2 + y.^2)));
Ixy = -sum(sum(rho .* x .* y));
Ixz = -sum(sum(rho .* x .* z));
Iyz = -sum(sum(rho .* y .* z));
Iyx = Ixy;
Izx = Ixz;
Izy = Iyz;

% 输出转动惯量
I = [Ixx Ixy Ixz; Iyx Iyy Iyz; Izx Izy Izz];
disp("转动惯量矩阵：");
disp(I);

% 绘制球体图形
figure;
surf(x, y, z);
axis equal;

代码说明：

1. 用户通过 `input` 函数输入球体的半径和分辨率。
2. 使用 `sphere` 函数生成球体的网格，再通过缩放和平移操作将球体放置于坐标系原点。
3. 计算球体的质量分布，这里假设球体质量分布均匀，因此每个网格点的密度都为 1。
4. 根据转动惯量的定义，利用 `sum` 函数计算球体的转动惯量元素。
5. 输出转动惯量矩阵。
6. 使用 `surf` 函数绘制球体的图形，并通过 `axis equal` 使得三个坐标轴的尺度相同。