"""Проверить энергию двухкатушечной матрицы.

Среда: встроенный MicroPython FEMM 5.
Самостоятельный аналитический контроль; полевой решатель не вызывается.
"""
# Собственная симметричная линейная матрица; индуктивности в Гн, токи в А.
l11, l22, mutual = .020, .080, .032
i1, i2 = 2.0, -.5
psi1 = l11*i1 + mutual*i2
psi2 = mutual*i1 + l22*i2
# Проверяем энергию общего режима и положительность определителя.
energy = .5*(i1*psi1 + i2*psi2)
print("Потокосцепления, Вб:", psi1, psi2)
print("Энергия, Дж:", energy)
print("Определитель, Гн²:", l11*l22 - mutual*mutual)