Описывается эффективный алгоритм нахождения такого поля, которое позволяет обеспечивать как устойчивость движения, так и изохронность ускорения во всем необходимом диапазоне радиусов. Рассматривается процедура оптимального конструирования изохронного магнитного поля в ускорителях с азимутальной вариацией, обеспечивающего фазовую устойчивость. Алгоритм реализован на языке FORTRAN для двух возможных вариантов периодичности магнитной структуры: с периодом Т = 2π/N и суперпериодом Т = 2π. Для этой цели используется процедура нахождения статических равновесных орбит и периода обращения на них ускоряемых частиц. В качестве сопутствующих результатов можно получить частоты радиальных и аксиальных колебаний, характеризующие устойчивость движения, а также другие параметры, представляющие практический интерес. Разработанные пакеты программ являются совершенно автономными и могут быть поставлены на любую платформу. Отличительной особенностью алгоритма является то, что он реализован в системе координат, используемой физик