SynopsisWe describe a procedure for the solution of the time-dependent Dirac equation. The procedure is based on the relativistic generalization of the matrix iteration method. We use this procedure to study electron momentum distribution along the laser beam propagation direction for the process of the tunneling ionization of hydrogen atom. We found, in agreement with the experimental observations [1], that relativistic effects lead to appreciable deviation of the distribution from the strict left-right symmetry present in the non-relativistic case. The expectation value of the momentum along the laser beam propagation direction grows linearly with intensity and follows closely the behavior of the expectation value of the kinetic energy divided by the speed of light. These features agree with the experimental results [1].