The longitude development of the muonic component in the extensive air shower can help to determine mass composition of cosmic rays. By using the timing information of muons, the production positions of muons can be reconstructed .$X_{max}^{\mu}$ which is the position with muon production reaching maximum is sensitive to the mass composition of cosmic rays. In this paper, the reconstruction of muon’s geometrical production heights along the shower axis by using their timing information are studied for cosmic rays with energies from 1PeV to 10 PeV by CORSIKA data. According to our analysis, the resolution of $X_{max}^{\mu}$ is 500 $g/cm^2$ with a resolution about 32 $g/cm^2$ for showers initiated by protons with energy 10 PeV ,and zenith angle θ= 45°.