=cot((arccos(1-(E_0*(E_1-E_2)/(E_1*E_2))))/2)