We consider nonlocal modification of the Einstein theory of gravity in framework of the pseudo- Riemannian geometry. The nonlocal term has the form H(R)F(□)G(R), where H and G are differentiable functions of the scalar curvature R, and F(□) = ∞Σn=0 fn□n is an analytic function of the d’Alambert operator □. Using calculus of variations of the action functional, we derived the corresponding equations of motion. The variation of action is induced by variation of the gravitational field, which is the metric tensor gμv. Cosmological solutions are found for the case when the Ricci scalar R is constant.