import numpy as np


R_etm = 2303
R_itm = 1940
L = 3995


ge = 1-L/R_etm
gi = 1 -L/R_itm
c = 299792458
fsr = c/(2*L)
fhom = fsr * np.arccos(np.sqrt(ge*gi))/np.pi