A = 10*log[10,1+epsilon^2*cosh(N*arccosh(omega_g/omega_p))^2] solve for N