编写程序，用Matlab数值模拟的方法，模拟磁镜现象，即沿带电粒子束进入方向，磁场中间弱两端强。自行设定参数，观察并记录带电粒子束的运动轨迹图。

以下是一个简单的 Matlab 磁镜模拟代码：

% 定义常数
m = 9.10938356e-31; % 电子质量
q = -1.60217662e-19; % 电子电荷
B0 = 0.5; % 磁场强度
B1 = 1.0; % 磁场强度

% 定义初始条件
x0 = 0;
y0 = 0;
z0 = 0;
vx0 = 1e5;
vy0 = 1e5;
vz0 = 0;

% 定义时间步长和总时间
dt = 1e-10;
t = 1e-5;

% 初始化位置和速度数组
x = zeros(1, t/dt);
y = zeros(1, t/dt);
z = zeros(1, t/dt);
vx = zeros(1, t/dt);
vy = zeros(1, t/dt);
vz = zeros(1, t/dt);

% 初始条件
x(1) = x0;
y(1) = y0;
z(1) = z0;
vx(1) = vx0;
vy(1) = vy0;
vz(1) = vz0;

% 计算电子运动轨迹
for i = 1:t/dt-1
    if x(i) < -1e-4
        B = B1;
    elseif x(i) > 1e-4
        B = B1;
    else
        B = B0;
    end
    ax = q*B/m*vy(i);
    ay = -q*B/m*vx(i);
    az = 0;
    vx(i+1) = vx(i) + ax*dt;
    vy(i+1) = vy(i) + ay*dt;
    vz(i+1) = vz(i) + az*dt;
    x(i+1) = x(i) + vx(i+1)*dt;
    y(i+1) = y(i) + vy(i+1)*dt;
    z(i+1) = z(i) + vz(i+1)*dt;
end

% 绘制电子运动轨迹
figure;
plot(x, y);
axis equal;
title('Electron Trajectory in Magnetic Mirror');
xlabel('x');
ylabel('y');

在这个代码中，我们将磁场设置为两端强中间弱的磁镜结构，即中间为 $0.5T$ 的磁场，两端为 $1.0T$ 的磁场。根据电子在磁场中的运动方程，我们可以在计算电子的运动轨迹时根据位置调整磁场强度。运行代码后，程序将绘制出电子在磁镜中的运动轨迹。你可以根据需要调整初始条件和磁场参数，观察不同情况下电子的运动轨迹。