帮我写一个画卢瑟福散射的入射粒子能量和散射角的坐标图的matlab代码

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

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

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

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

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

% 绘制入射粒子能量和散射角的坐标图
figure;
plot(x(:), y(:), '.', 'MarkerSize', 1);
xlabel('Scattering Distance (m)');
ylabel('Incident Particle Energy (MeV)');
title('Rutherford Scattering');

这段代码将会绘制出卢瑟福散射的入射粒子能量和散射角的坐标图，其中每个点表示一个散射事件，横坐标为散射距离，纵坐标为入射粒子能量。您可以根据自己的需要调整代码中的参数来得到更符合您需求的图形。