"""Постановка оптимизации: физическая модель, критерии и ограничения.

Среда: встроенный Python FEMM 5.
Проверяется только авторский аналитический эталон, полевой решатель не вызывается.
Полученные числа не являются квитанцией расчёта электромагнита или машины.
"""

# Силу задаём в ньютонах, индукцию — в теслах, глубину — в миллиметрах.
# Из магнитного давления получаем площадь в м², затем переводим её в мм² и находим ширину.
import math
mu0 = 4 * math.pi * 1e-7
force_n = 100.0
b_limit_t = 1.0
depth_mm = 20.0
area_mm2 = 2 * mu0 * force_n / b_limit_t**2 * 1e6
print("Площадь, мм²:", area_mm2)
print("Ширина, мм:", area_mm2 / depth_mm)