A method of probabilistic approximation of the operator e−itH, where H=−12d2dx2+V(x), V∈L∞(R), in the strong operator topology is proposed. The approximating operators have the form of expectations of functionals of sums of independent identically distributed random variables.