The stored energy distribution versus crystal orientation in polycrystalline copper was determined using synchrotron radiation. This distribution is an important input data for recrystallization models. The stochastic vertex model of recrystallization was used in the present work. It is a mixture of the classical vertex model and the Monte Carlo algorithm. Both grain boundary energy and stored energy are taken into account in the calculations. In each elementary step, a reasonably small, random modification of a given vertex position is generated and a corresponding total energy change of a system is calculated. A new vertex position is retained with a probability proportional to the Boltzmann factor. In such a way one avoids solving a complex system of equations. This approach is also closer to the stochastic nature of recrystallization process. The inclusion of the stored energy distribution in the above model enables a good explanation of the recrystallization process. The recrystallization textures for polycrystalline copper rolled to low and high reductions were predicted in agreement with experimental results.