Faraday Effect
ManipulatePlot[
(*SpB[*)Power[E(*|*),(*|*)I \[Omega] (-1+(*SqB[*)Sqrt[1+(*FB[*)((f)(*,*)/(*,*)((-I \[Gamma]-\[Omega]) \[Omega]+(*SpB[*)Power[(\[Delta]-\[Omega]0)(*|*),(*|*)2](*]SpB*)))(*]FB*)](*]SqB*))](*]SpB*)+(*SpB[*)Power[E(*|*),(*|*)I \[Omega] (-1+(*SqB[*)Sqrt[1+(*FB[*)((f)(*,*)/(*,*)((-I \[Gamma]-\[Omega]) \[Omega]+(*SpB[*)Power[(\[Delta]+\[Omega]0)(*|*),(*|*)2](*]SpB*)))(*]FB*)](*]SqB*))](*]SpB*) /. {\[Gamma]->0.694, \[Omega]0->50.0} // Abs
, {ω, 20, 80}, {{f,5},0,100,5}, {{\[Delta],0},0,10,1}
, FrameLabel->{"wavenumber", "transmission"}
, Frame->True
]