We present a numerical procedure for recovering of photon statistics of few-photon pulsed laser sources. It provides controlled level of mean error for both of photon and photocounting statistics. It is shown that procedure is correct for model data with average error level of 10−4 and fidelity 0.9998.