In the Starobinsky inflationary model inflation is driven by quantum corrections to the vacuum Einstein equation. We reduce the Wheeler-DeWitt equation corresponding to the Starobinsky model to a Schroedinger form containing time. The Schroedinger equation is solved with a Gaussian ansatz. Using the prescription for the normalization constant of the wavefunction given in our previous work, we show that the Gaussian ansatz demands Hawking type initial conditions for the wavefunction of the universe. The wormholes induce randomness in initial states suggesting a basis for time-contained description of the Wheeler-DeWitt equation.