帮我写一个画卢瑟福散射的入射能量和散射角的关系图的matlab代码

下面是一个简单的绘制卢瑟福散射的入射能量和散射角的关系图的Matlab代码：

% 定义原子核的电荷数
Z = 79;

% 定义散射角度的范围和步长
theta = linspace(0, pi, 1000);

% 定义入射粒子的能量范围和步长
E = linspace(1, 10, 1000);

% 计算不同入射能量下的散射角度
for i = 1:length(E)
    for j = 1:length(theta)
        % 计算卢瑟福散射公式
        r = 1 / (4 * pi * 8.854e-12 * 1.602e-19^2 * Z) * E(i) / sin(theta(j) / 2)^2;
        % 计算散射角度对应的位置
        x(i, j) = -r * sin(theta(j));
        y(i, j) = r * cos(theta(j));
    end
end

% 绘制入射能量和散射角的关系图
figure;
hold on;
for i = 1:length(E)
    plot(x(i, :), y(i, :), 'LineWidth', 1.5);
end
xlabel('Scattering Angle (m)');
ylabel('Scattering Distance (m)');
title('Rutherford Scattering');

这段代码将会绘制出卢瑟福散射的入射能量和散射角的关系图，其中不同颜色的曲线表示不同入射能量下的散射情况。您可以根据自己的需要调整代码中的参数来得到更符合您需求的图形。