The critical relaxation of the magnetization in a two-dimensional Ising model with quenched random non-magnetic impurities has been studied by numerical simulation. A squared lattice of size 4002 with spin concentrations p=1.0; 0.95, 0,9, 0.85, 0.8, 0.75 and 0.7 was considered. The dynamical critical exponent z was determined by the Monte Carlo method combined with the dynamical renormalization group method. The following values were received for z(p): z(1)=2.24+or-0.07, z(0.95)=2.24+or-0.06, z(0.9)=2.24+or-0.06, z(0.85)=2.38+or-0.05, z(0.8)=2.51+or-0.06, z(0.75)=2.66+or-0.07 and z(0.7)=2.88+or-0.06. A singular dynamic scaling behaviour with z=A' mod ln(p-pc) mod +B' was found for systems with spin concentrations p