Newton’s theory does not include inductive phenomena, but a relativistic theory of gravitation should include them. Indeed, under the relativistic mass-energy equivalence, not only the mass is a source of gravitational field but any kind of energy also is. Therefore, a body creates gravitational field not only by mass but also by their kinetic energy, i.e. by their movement. And this, ultimately, is what it means induction: the production of forces by moving bodies  .
In the general relativity theory, Einstein predicted the existence of gravitational induction phenomena, such phenomena are appointed by Einstein as gravitomagnetism. It can be showed that Jefimenko equations are also derived from linearized Einstein equations (see, for example, pp 47 and 48 in  ).
Since in 2004, NASA has orbited the “Gravity Probe B”, whose purpose was to prove the existence of gravitomagnetism, (see   ).
In order to describe the time-dependent gravitational systems, the Jefimenko’s generalized theory of gravitation is based on postulating of retarded expressions for the accustomed gravitational field g and the Heaviside’s or cogravitational field K (Heaviside  was the first who supposed the existence of this field making an analogy between gravitational and electromagnetic fields), where the field g acts to and arises from motionless as well as moving masses, and the field K acts to and arises from exclusively moving masses. Let us from this point call this theory “gravitodynamics”, and the complex of the fields g and K call the “gravitodynamical field”. Jefimenko taking into account mentioned retarded expressions for g and K, obtains the system of equations for which these expressions are solutions. These equations are analogous to Maxwell’s equations. In principle, this method was proposed to eliminate the possibility of instantaneous solutions from the discussion. There against in this work we are going to postulate the system of differential equations rather than solutions for the gravitational dynamics, and we will obtain both retarded and instantaneous solutions for the fields g and K. The gravitational field g behaves analogously to the electric field in the Maxwell’s electromagnetic theory, and cogravitational field K analogously to the magnetic field.
First of all, we write the equations describing time-dependent gravitational systems  (see p. 120) and 
where c is the velocity of propagation of the fields, which is supposed equal to the velocity of light for the retarded component, G is the constant of gravitation, is the density of mass and is the density of mass current.
There are some differences between Maxwell’s equations of electrodynamics and the Jefimenko’s equations of gravitation, i.e. the analogy is not perfect. For example, we have two kinds of electric charges, positives and negatives, which repel each other if the charges are equal and attract each other if they are different, whereas while we have only one type of mass, and if we have a system of two masses in repose, they always attract each other. While the electric field is directed from positive charges generating this field and is directed to the negative charges, the gravitational field is always directed to the masses by which is created. Another difference is that the magnetic field is always right- handed relative to the electric current by which is created, while the cogravitational field is always left-handed relative to the mass current by which is created.
In the analogy between electrodynamics and the so-called gravitodynamics, following the Jefimenko’s book  we can resume the correspondence between electromagnetic and gravitodynamic symbols and constants in the following Table 1.
2. The Gravitodynamical Potentials
Here, we introduce as is made in electrodynamics, the gravitodynamical potentials. If
Table 1. Corresponding electromagnetic and gravitodynamic symbols and constants.
the cogravitational field K satisfies Equation (2), we can always write it as the curl of some other vector quantity,
where is the gravitodynamical vector potential. Substituting Equation (5) in (3), we obtain
The quantity within the parentheses can be written as the gradient of a gravitody- namical scalar potential:
Substituting the expressions (5) and (8) for the fields and, in the inhomogeneous Equations (1) and (4), we obtain
Equations (9) and (10) can be decoupled choosing the appropriate form of the potentials and. Moreover, if we simultaneously make the transformations
in (5) and (8), we get the same original fields and. Here, is an arbitrary scalar function. We can choose this function in order to impose an additional condition over and, in a similar way like the Lorentz or Coulomb gauge in the electromagnetic field, namely,
which allows us to separate Equations (9) and (10) for the potentials and. These potentials depend on the gauge condition we chose.
3. Jefimenko’s Equations for the Solenoidal and Irrotational Components
Following the ideas of the work of Chubykalo et al.  -  and using the analogy between the Maxwell’s equations and the Jefimenko’s ones   , we will apply the Helmholtz’s theorem to define potentials that are independent of gauge transformations.
The Helmholtz’s theorem claims that under certain conditions all vector fields can be represented as the sum of an irrotational and a solenoidal components. We will use this theorem to separate the fields and.
Therefore, here we state the Helmholtz’s theorem as  :
If the divergence and a curl of a vector function are specified, and if they both go to zero faster than as, and if itself tends to zero as, then is uniquely given by
1For systems localized in a finite region of space, it is evident that the fields g y K depend on r as 1/r2.
We are going to suppose that all conditions of this theorem are satisfied by the fields and defined by Equations (1) to (4)1, and then, we apply Helmholtz’s theorem to these quantities, including. Thus, we obtain
where the indices “i” and “s” mean irrotational and solenoidal components of the vectors, respectively.
We are going to substitute and given by the Equations (17)-(19) into the Jefimenko’s Equations (1)-(4) and then, we obtain for the irrotational part:
and the next equations for the solenoidal part:
4. The Gravitodynamical Potentials from Helmholtz’s Theorem
By definition, for the irrotational component of the gravitational field we can define the scalar potential as
and if we substitute this relation into Equation (22), we obtain the Poisson’s equation
Apparently, we need to take into account that is not completely defined only by the Poisson’s Equation (29), because we have another differential equation for, which can be obtained by substituting (28) into Equation (23)
We show now that Equation (30) is equivalent to the law of conservation of mass. Indeed, let us take the divergence of the Equation (23), then we obtain as the result
But from Equation (22) and because, Equation (31) becomes the conservation mass law or the continuity equation
Now, we will demonstrate that the solution of Equation (30), indeed, is the same solution of the Poisson’s Equation (29). To do this, we note that the irrotational component of can be written as
where the potential is defined as
and where, if we relate Equations (32), (23), (28) y (33) and the fact that, we have
And from (36) and (34), we obtain
which is the solution of the Poisson’s Equation (29). So we have found that the Poisson’s equation given by Equation (29), completely defines the potential, together with its boundary conditions.
Since by definition, then Equations (24) and (25) have the trivial solution.
Let us now apply the Helmholtz’s theorem to the vector potential. The Helmholtz theorem it is also known as the fundamental theorem of vector calculus (see Section III), and allows us to decompose every vectorial field in two components, an irrotational and a solenoidal one. Intuitively, it says that every vector function can be written as the sum of a divergence-free function (like) and a curl-free function (like), so that there exist scalar and vector potentials. So that, because this is the form in which we can easily solve the system formed by Equations (26) and (27). Supposing that
and taking into account (26) we obtain
One can substitute Equations (39) and (40) into (27), and we obtain
where we used the vector identity for any arbitrary vector V.
We have found that system of Equations (1)-(4) reduces to Equations (29) and (41), applying the Helmholtz’s theorem. Therefore, we obtain separated equations for vector and scalar potentials, namely,
5. Invariance of the Potentials Gs and Τ under Gauge Transformations
Now, we will show that the potentials and (which are gravitodynamical countertypes of the electromagnetic potentials and from  ) are invariant under gauge transformations (11) and (12), in common with the fields g and K. This will be the most prominent property of these potentials.
If we apply the Helmholtz theorem to the gravitational and cogravitational fields in terms of the ordinary potentials given by (5) and (8) without taking into account any gauge condition, we have
then, by definition, is given by
where is an scalar function. If we substitute (48) into (44) then
and from Equations (49) and (22) we have the relation between and.
Now, we apply the Helmholtz theorem and the gauge transformations (11) and (12) and from
comparing the solenoidal parts we obtain
If we seek the transformation law for, then we can obtain the other transformation law for.
From Equation (51), we have the irrotational part of
and including Equation (48) in (53) we get
At last, we can use Equation (50) for
Hence, we have checked that and are invariants under gauge transformation and we can see that any gauge transformation is irrelevant if we use the Helmholtz’s theorem.
It is convenient to remark that the fields g and K are generated only by and given by (28), (39) and (40), so we can consider and as the potentials generating the gravitodynamical field g and K.
And so, we have shown that it is possible to define vector as well as scalar gravitody- namical potentials, which are invariant under gauge transformation. These potentials are defined uniquely from their differential Equations (42) and (43). For this reason, we have arguments for supposing the physical reality of these potentials, similarly to the fields g and K and unlike the gravitational potentials introduced by Jefimenko in  , which are only used as a mathematical tool for calculating gravitational and cogravitational fields.
Our scalar potential T is a generator of the so-called instantaneous action at a distance in gravitation, and the vector potential can propagate with the velocity of light and it is responsible for the retarded action of the gravitodynamical field. So, one can conclude that retarded interaction in gravitodynamics takes place not instead but together with instantaneous action at a distance.