The goal of the paper is to represent iterated Stratonovich stochastic integrals using orthonormal expansions of functions by the Walsh series, to construct an algorithm for modeling iterated Stratonovich stochastic integrals, and to test this algorithm by numerical experiments. Obtained results may be used in the constructing high-order numerical methods based on the Taylor –Stratonovich expansion and the Taylor– It^o expansion. Therefore, they can be applied in a lot of problems of physics, engineering, economics and finance, when it is necessary to solve numerically stochastic differential equations.