In well-known statistical models of global optimization only values of objective functions are taken into consideration. However, efficient algorithms of local optimization are also based on the use of gradients of objective functions. Thus, we are interested in a possibility of the use of gradients in statistical models of multimodal functions, aiming to create productive algorithms of global optimization.