"""Встроенный график аналитической зависимости потокосцепления от тока.

Вход: пять токов [А], L = 0,12 Гн; выход: график Ψ [Вб·виток] от I [А].
В FEMM GUI открывается окно графика; headless-запуск проверяет выполнение
вызова, но не внешний вид окна. Новое полевое решение не создаётся.
"""

import femm

# 1. Ось абсцисс — токи в заданном порядке, включая нулевое возбуждение.
currents_a = [0.0, 1.0, 2.0, 3.0, 4.0]

# 2. Формируем ординаты Ψ = LI одним проходом по тем же токам.
linkage_wb = [0.12 * current_a for current_a in currents_a]

# Проверяем правую конечную точку до передачи данных построителю.
assert abs(linkage_wb[-1] - 0.48) < 1e-12

# 3. Передаём два списка одинаковой длины: сначала x, затем y.
femm.plot(currents_a, linkage_wb)
print('Пять точек: по горизонтали ток, А; по вертикали потокосцепление, Вб.')
