In the paper we discussed the evolution equations for diffractive production in the framework of CGC/saturation approach, and found the analytical solutions for several kinematic regions. The most impressive features of these solutions are, that diffractive production does not manifest geometric scaling behaviour i.e. being a function of one variable. Based on these solutions, we suggest an impact parameter dependent saturation model, which is suitable for describing diffraction production both deep in the saturation region, and in the vicinity of the saturation scale. Using the model we attempted to fit the combined data on diffraction production from H1 and ZEUS collaborations. We found that we are able describe both xIP and β dependence, as well as Q behavior of the measured cross sections. In spite of the sufficiently large χ2/ d. o. f. we believe that our description provides an initial impetus to find a fit of the experimental data, based on the solution of the CGC/saturation equation, rather than on describing the diffraction system in simplistic manner, assuming that only quark-antiquark pair and one extra gluon, are produced.