In this paper the problem of optimization of multivariate multimodal functions observed with random error is considered. Using the random function for a statistical model of the objective function the minimization procedure is suggested. This algorithm is convergent on a discrete set. To avoid computational difficulties, the modified algorithm is defined by substituting the parameters of minimization procedure by their estimates.