把上面代码改成可以生成图形

如果你想将上述计算结果可视化为环形电流产生的磁场图，你可以使用MATLAB的`polar`函数结合`surf`或`quiver`等绘图函数。这里提供一个简化的例子，展示如何生成一个简单的环形磁场图像：

function plot_magnetic_field(i, r, mu0)
    % 定义环形电流区域的半径范围
    [r_min, r_max] = linspace(0, r, 100);

    % 创建角度向量
    theta = linspace(0, 2*pi, 100);

    % 计算每个角度对应的磁场强度
    [r_values, theta_values] = meshgrid(r_min, theta);
    B_values = magnetic_field(i, r_values, mu0);

    % 绘制磁场分布
    figure;
    colormap('parula'); % 更改颜色映射
    surfc(theta_values, r_values, B_values, 'EdgeColor', 'none');
    xlabel('\theta (弧度)');
    ylabel('r (米)');
    zlabel('B (特斯拉)');
    title(['环形电流磁场 (i = ', num2str(i), ' A, r = ', num2str(r), ' m)']);
end

在这个示例中，我们创建了一个二维网格，代表了磁场强度随半径r和角度θ的变化。然后，调用`surfc`函数绘制三维表面图，并显示各个轴的标签和图表标题。

运行`plot_magnetic_field(i, r, mu0)`即可看到磁场分布的图形。