We designed a single polyethylene bonner sphere with several thermo-luminescence dosimeters (TLD), for measurement of neutron energy spectrum. For the separation of the neutron dosage in the neutron-gamma mixed field, we used 21 ea TLD-600s and TLD-700s, respectively. Because, TLD-600 is sensitive to neutron and gamma rays, and, TLD-700 is sensitive only to gamma-rays, we could determine the each dose by neutron and gamma rays. The neutron response function of the bonner sphere with TLDs was calculated by MCNPX (ver. 2.5.0) Monte Carlo simulation in the energy range from $10^{-1}$ to 20 MeV. For the Cf-252 standard neutron source in KRISS, we could estimate the neutron energy spectrum by unfolding method using the response function.