Accurate ab initio calculations of the ratio of double-to-single ionization of Ne atoms in strong laser fields are difficult due to the many-electron nature of the target. Here, with accurate total cross sections carefully evaluated using the state-of-the-art many-electron R-matrix theory for both electron impact ionization and electron impact excitation of Ne + , we simulate the total double ionization yields of Ne 2+ in strong laser fields at 780 nm and 800 nm for pulse durations in the range from 7.5 fs to 200 fs based on the improved quantitative rescattering model. The corresponding single ionization yields of Ne + are calculated within the nonadiabatic tunneling model of Perelomov, Popov, and Terent'ev. The ratio of double-to-single ionization of Ne is then obtained from the calculated double and single ionization yields. By normalizing the ratio to the one calculated from solving the time-dependent Schrödinger equation for a few-cycle short pulse, we make quantitative comparisons of our results with experimental data to show that our model predicts the experimental findings very well. Finally, we analyze the pulse-duration dependence of the double-to-single ionization ratio.