Написать код программы в фортран

Электрон совершает дрейфовое движение в неоднородном магнитном поле, градиент которого направлен вдоль оси Y. Радиус циклотронного вращения R=0.5 см. В какой-то момент времени, например, после 3-х циклотронных оборотов, направление дрейфа меняется на обратное с удвоенной величиной скорости дрейфового движения. Определить поля, в которых находится частица (их направления и величины). Величина начальной скорости электрона равна 1*109 см/с. Получить график зависимости энергии частицы (в эВ) от времени.