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