We analytically calculate to second order the correction to the asymptotic form of quasinormal frequencies of four dimensional Schwarzschild black holes based on the monodromy analysis proposed by Motl and Neitzke. Our results are in good agreement with those obtained from numerical calculation.