> ReF := (Omega,Q) -> 1/((1-Omega^2)^2+(Omega/Q)^2);
> ImF := (Omega,Q) -> Omega * ((1-Omega^2) * Q-1/Q) / ((1-Omega^2)^2+(Omega/Q)^2);
> ReFO := (Omega,Q) -> 1/((1-Omega^2)^2+(Omega/Q)^2)/Q^2;
> ImFO := (Omega,Q) -> Omega * ((1-Omega^2) * Q-1/Q) / ((1-Omega^2)^2+(Omega/Q)^2)/Q^2;
> plot({[ReFO(10^x,1), ImFO(10^x,1), x=-infinity..infinity], [ReFO(10^x,3), ImFO(10^x,3), x=-infinity..infinity], [ReFO(10^x,10), ImFO(10^x,10), x=-infinity..infinity]} ,0..2.5,-1.1..0.6,numpoints=100,title=‘Q=1,3,10‘);
> plot({[ReF(10^x,1), ImF(10^x,1), x=-100..100], [ReF(10^x,3), ImF(10^x,3), x=-100..100], [ReF(10^x,10), ImF(10^x,10), x=-100..100]}, 0..100,-55..55,numpoints=10000,title=‘Q=1,3,10‘);