Double differential neutron yields from carbon target irradiated by 5 MeV and 9 MeV deuterons were measured at the Kyushu University Tandem Accelerator Laboratory. A carbon target was thick enough to stop the incident deuterons in the target. An NE213 liquid organic scintillator was employed to detect neutrons emitted from the target and placed at several directions from 0° to 140°. To estimate the contribution of scattered neutrons from the floor and walls, neutron yields were measured with an iron shadow bar located in front of the scintillator. The measured energy spectrum was derived by an unfolding method with FORIST and MAXED codes using the response functions of the NE213 scintillator calculated with the SCINFUL-QMD code. The incident energy dependence of total thick target neutron yields for carbon was investigated over the energy range from 5 to 50 MeV and the empirical formula was obtained. The measured neutron energy spectra were also compared with the calculations by PHITS code.