For details see the header of the comptsTabgb_ds.m

<<comptsTabgb_ds.m

zamiana={beta[x__]->1-Tanh[x],omega[x__]->1+Tanh[x]}

Ttt=ABGBTtt/.zamiana;

Trr=ABGBTrr/.zamiana;

Tang=ABGBTang/.zamiana;

First, let us check if the stress-energy tensor is covariantly conserved. The only nonzero component of the covariant
divergence of the tensor  is given by

CovariantDivergence=-1/2 D[f[r],r]   Ttt + 1/2 D[f[r],r]   Trr +D[Trr,r] f[r] +2/r Trr f[r]-2/r Tang f[r];

f[r_]:=1-2 M/r (1-Tanh[Q^2/(2 M r)] )- eps l^2 r^2/3

h[r_]:=1/f[r]

list={cc1,cc2,cc3,cc4,cc5,cc6,cc7,cc8,cc9,cc10}

CovariantDivergence/.zamiana//Expand;

It takes a few seconds to establish the covariant conservation.

Map[Simplify,Collect[CovariantDivergence,list]]

TimeUsed[]

To calculate the components of the stress-energy tenosor for a given spin
one has to specify coefficients cc1,...,cc10  (use:  <<scalar.m, <<spinor.m or <<vector.m ).
The result must be divided by 96 pi^2 m^2.

(* <<scalar.m *)

<<spinor.m;
(* vector.m *)

Ttt;Trr,Tang;