Let Ric(g) and S(g) be the Ricci tensor and Ricci scalar for the metric g.  Then the covariant form of the Einstein tensor is Ein(g) = Ric(g) - 1/2*g*S(g).  The contravariant form is obtained by raising both indices of the covariant Einstein tensor with the metric g.

The program EinsteinTensor(g, R) returns the contravariant form of the Einstein tensor.  This tensor is symmetric and its covariant divergence vanishes.

This command is part of the DifferentialGeometry:-Tensor package, and so can be used in the form EinsteinTensor(...) only after executing the command with(DifferentialGeometry) and with(Tensor) in that order.  It can always be used in the long form DifferentialGeometry:-Tensor:-EinsteinTensor.