  A class of matrix models which arises as partition function in U(N)
Chern-Simons matter theories on three sphere is investigated. Employing the
standard technique of the 1/N expansion we solve the system beyond the planar
limit. In particular we study a case where the matrix model potential has 1/N
correction and give a general solution thereof up to the order of 1/N^2. We
confirm that the general solution correctly reproduces the past exact result of
the free energy up to the order in the case of pure Chern-Simons theory. We
also apply to the matrix model of N=2 Chern-Simons theory with arbitrary
numbers of fundamental chiral multiplets and anti-fundamental ones, which does
not admit the Fermi gas analysis in general.
