We show that the ABJM theory, which is an
superconformal U(N) × U(N) Chern-Simons matter theory, can be studied for arbitrary N at arbitrary coupling constant by applying a simple Monte Carlo method to the matrix model derived by using the localization method. Here we calculate the free energy, and show that some results obtained by the Fermi gas approach can be clearly understood from the constant map contribution obtained by the genus expansion.