Для бесконечного прямолинейного проводника с током I магнитное поле равно B = (μμ
0/4π)(2I/ρ), где ρ - расстояние до проводника в цилиндрической системе координат. Сила Лоренца равна dp/dt = q[v,B]. Таким образом, получается система трёх ДУ (по компонентам момента), которые можно решать численно, например, по методу Эйлера.
Можно вместо поля прямолинейного проводника взять поле диполя, B = (μμ
0/4π)(3(m,r)r/r
5 - m/r
3), и промоделировать, таким образом, движение заряженных частиц в
магнитном поле Земли (также
магнитосфера). Заряженные частицы от Солнца (
солнечный ветер) могут захватываться магнитным полем (
радиационный пояс), т.е. двигаться в нём практичеки по окружностям вдоль силовых линий.
Северное сияние, кому интересно, может порождаться, насколько я понял, как непосредственно излучением вращающихся вдоль магнитных силовых линий частиц, так и излучением атомов атмосферы, возбуждённых при столкновении с этими частицами.