We described a method of resonant mode identification in dielectric-disc loaded cylindrical cavity resonators. The characteristic equations are solved using the ContourPlot graph of Mathematica. Contour graph method uses graphical method. It is comparable with numerical method. The numerical method is very difficult mode identification. The analysis is based on the approximated electromagnetic representation which is only concentrated on the calculation of resonant frequencies, and a mode identification of resonant frequencies has not been covered. However, It is possible to calculate precise resonant frequencies and to identify the mode of resonant frequencies using the contour graph method. The contour graph method is not a method using approximated representation of electromagnetic field variation at the outer area of dielectric in the resonators. It is a method using enact representation.