В работе рассматривается эффективный метод оптимизации формы полюсов в мультипольных магнитах с ферромагнитным ярмом, целью которого является формирование заданного распределения индукции магнитного поля в рабочей области магнита. Предполагается установка данного квадрупольного магнита в фокусирующих элементах накопительного кольца CR, входящего в состав нового ускорительного комплекса FAIR, создаваемого в настоящее время в Дармштадте, Германия. Основной идеей метода является представление характеристик магнитного поля в виде разложения в ряд Фурье и последующая компенсация нежелательных гармонических составляющих напряженности магнитного поля за счет деформации поверхности полюсного наконечника. Высокая скорость сходимости и устойчивость оптимизационной процедуры объясняется значительной корреляцией между амплитудами гармонических составляющих разложения магнитного поля в рабочей области и специфическими характеристиками формы поверхности полюса. Разработанная технология поиска оптимальной формы встроена в программу численного моделирования магнитных полей, использующую метод конечных элементов второго порядка и предназначенную для расчета характеристик магнитного поля в электромагнитных устройствах различного назначения.