r = sqrt((X(i,j,k)^2+Y(i,j,k)^2+Z(i,j,k)^2)); % 场点到长方体中心距离位置 3 处的索引超出数组边界(不能超出 6)。以上程序怎么修改

你可以在计算 r 的时候加上边界判断条件，避免出现数组越界的情况。具体的代码修改如下：

if i >= 1 && i <= 6 && j >= 1 && j <= 6 && k >= 1 && k <= 6
    r = sqrt((X(i,j,k)^2+Y(i,j,k)^2+Z(i,j,k)^2)); % 场点到长方体中心距离位置 3
else
    r = NaN; % 超出数组边界，将距离设为 NaN
end

这样，如果索引超出了数组边界，就会将距离 r 的值设为 NaN，避免程序出现错误。