#!/usr/bin/python
# Uloha MAGN1
# Vytvoril Radek Podgorny <radek@podgorny.cz>

from math import *

Eps0 = 8.854187817 * 10**-12 #epsilon
mi0 = 4*pi * 10**-7


# Tady zadej svoje parametry
A = 0.022
C = 0.0428
D = 0.031
F = 0.263
G = 0.236
H = 0.0719
N1 = 96
N2 = 18
m1 = 0.000791
m2 = 0.00214
mir = 23900
Kp = 0.83
HL = 0.0241
Bm = 0.8
f = 138
In1 = 0.00022
In2 = 0.000624
Iss = 0.00491


RA = (H-m2)/(mi0*mir * A*HL) + m2/(mi0 * A*HL)
RC = (H-m1)/(mi0*mir * C*HL) + m1/(mi0 * C*HL)
RD = H/(mi0*mir * D*HL)
RF = F/(mi0*mir * A*HL)
RG = G/(mi0*mir * A*HL)

Rleft = 2*RG+RD
Rright = 2*RF+RA
R = RC + (Rleft*Rright)/(Rleft+Rright)
L1 = N1**2/R

print("Vlastni indukcnost L1 je " + str(L1))


### KURVADRAT! Celou dobu to asi fungovalo a ja delal tolik vifikundaci! Vracet uz to nebudu...
#Rright = 2*RF + RA
#R = RD + 2*RG + (RC*Rright)/(RC+Rright)
#L2 = N2**2/R
#M = (L1*L2)**0.5

### Tohle je vykradeny z nejakyho prikladu ze skript
R1 = RC
R2 = 2*RF + RA
R3 = 2*RG + RD
M = (N1*N2 * R2)/(R1*R2 + R1*R3 + R2*R3)

print("Vzajemna indukcnost M je " + str(M))


### Tohle je zatim blbe...
#Rright = 2*RF + RA
#R = RD + 2*RG + (RC*Rright)/(RC+Rright)
I = Bm*R/N1

print("Maximalni I1 je " + str(I))

#Spravny vysledky pro zbytek:
#2.42049
#0.00127938