The numerical model of the elongation kinetics of the single-stranded DNA has been proposed and studied. The proposed model is based on the principles of the phenomenological modelling and implemented as diffusion cellular automata. The analytical and numerical estimates of the elongation reaction rate for different concentrations of nucleotides has been obtained. The possibility of nucleotide binding delays with reducing the concentration of one of the nucleotides in the solution is shown.