Комментарий

В правой части уравнения записан вектор действующей на частицу силы. Этот вектор равен сумме (умноженной на заряд) двух векторов: вектора электрического поля и вектора, являющегося векторным произведением скорости частицы на магнитную индукцию. Это векторное произведение записано с помощью описанной в самом начале процедуры vprod().

Далее записанное уравнение следует решить. Для этого формируем последовательность из трех элементов-уравнений. Каждое такое уравнение описывает динамику частицы вдоль соответствующей координатной оси. Чтобы составить эти уравнения, выбираем из правой и левой частей исходного векторного уравнения элементы списков с одинаковыми индексами и приравниваем их.

Доступ к левой части векторного уравнения осуществляется с помощью команды lhs(VecEq), а к правой - rhs(VecEq). На операнды эти части разбиваются посредством процедуры ор(). Первым операндом как в правой, так и в левой части уравнения являются скалярные множители: для левой части - это масса m, для правой - заряд е. На эти скаляры будут также множиться и новые уравнения. Вторые операнды правой и левой частей векторного уравнения - списки. Доступ к элементам списков реализуется путем указания индекса этих элементов, например op(lhs(VecEq))[2][3] - третий элемент второго операнда левой части уравнения VecEq.

После последовательности из уравнений, формируемой процедурой seq(), указана переменная IniCon, в которой записаны все начальные условия. Последовательность уравнений и начальных условий заключена в фигурные скобки, после которых указано множество функций параметра t, относительно которых следует решать систему.

> dsolve({seq(op(lhs(VecEq))[1]*op(lhs(VecEq))[2][i]=op(rhs(VecEq))[1]*op(rhs(VecEq))[2][i],i=1..3),IniCon},{x(t),y(t),z(t)});