“…The integral is taken over M if it converges; otherwise, one integrates over an arbitrarily large, relatively compact domain Ω containing supports of variations g t with g 0 = g. The physical meaning of (0.1) has been discussed in [1] for the case of a globally hyperbolic spacetime (M 4 , g) when D = Span(N ) and hence S mix = g(N, N ) Ric N,N , and the following have been obtained: the Euler-Lagrange equations, see (3.12), called the mixed gravitational field equations; sufficient conditions for existence of solutions for empty space, the conservation law analogous to conservation law of stress-energy tensor in relativity; equations of motion of a particle on isoparametric foliations in the "mixed gravitational field"; the linearized mixed field equations about the Minkowski metric; the value of coupling constant a in the weak field and low velocity limit. In this paper, we explore (0.1) for any spacetime, the obtained mixed gravitational field equations generalize the result of [1]. In fact, we work in arbitrary number of dimensions of a pseudo-Riemannian manifold endowed with a distribution, and also generalize certain results of [2,8], where the particular case of variations of metric (called adapted variations) has been examined.…”