Motion of a charged particle in a static magnetic field in matlab

  1. I want to simulate this type of situation but for simplicity I'm implementing it in 2d at first.. by setting the initial velocity 0 in parallel to B.

    1. The problem statement, all variables and given/known data
    mass of particle is given assume anything
    magnitude of B
    initial velocity

    2. Relevant equations

    3. The attempt at a solution
    % script that simulates a moving particle with some initial velocity in a
    % magnetic field B

    v0 = [5 0 0]; %initial velocity
    B = [0 0 -5]; %magnitude of B
    m = 5; % mass
    q = 1; % charge on particle
    r0 = [0 0 0]; % initial position of particle
    t = 0;

    % Now we want to find the next velocity as the particle enters the magnetic
    % field and hence its new position

    r = r0;
    v = v0;
    dt = 0.00000000000000000001;

    xlim([-25 25])
    ylim([-25 25])
    hold on

    for n = 1:100
    %plot it

    % update time
    t = t+dt;

    % new position r
    dr = v*dt;
    r = r + dr;

    %find new velocity
    dv = (q/m) * cross(v,B);
    v = v + dv;


    I know the most probable reason is the quantissation error but I'm not able to figure out that how can I avoid this error and simulate the situation in matlab ..
  2. jcsd
Know someone interested in this topic? Share a link to this question via email, Google+, Twitter, or Facebook