We all know that GRtensor works only with Maple V or older versions of it. Now, Maple have introduced the tensor with somewhat supplements the GRtensor package. Since I just started, I know that it can do the basic tensor calculations. It might also do a little bit more like defining new tensors and taking covariant derivatives. so finding extrinsic curvature etc, should not be a problem.
I have included a Maple command file below which demonstrates the ability of the tensor package to calculate the Ricci scalar from a given metric. Just edit the 2nd and 3rd to input the desired metric in 4D. The rest works fine.
=========================================
> with(tensor);
> gg := array(symmetric, sparse, 1 .. 4, 1 .. 4);
> gg[1, 1] := -f(r); gg[2, 2] := -1/gg[1, 1]; gg[3, 3] := r^2; gg[4, 4] := r^2*sin(theta)^2;
> eval(gg);
> gh := create([-1, -1], eval(gg));
> ginv := invert(gh, ‘detg’);
> coord := [t, r, theta, phi];
> d1g := d1metric(gh, coord);
> cfl1 := Christoffel1(d1g);
> cfl2 := Christoffel2(ginv, cfl1);
> d2g := d2metric(d1g, coord);
> RIEMN := Riemann(ginv, d2g, cfl1);
> displayGR(Riemann, RIEMN);
The Riemann Tensor
non-zero components :
1 d / d \
R1212 = – — |— f(r)|
2 dr \ dr /
> riciten := Ricci(ginv, RIEMN);
> rici := Ricciscalar(ginv, riciten);
[compts = (r^2*(diff(f(r), r, r))+4*(diff(f(r), r))*r-2+2*f(r))/r^2, index_char = []]