When using 

with(physics)

I can enter a metric tensor. If it is the Schwarzschild metric, it is easy enough to do

Geodesics(output=solutions)

to get a general geodesic solution. But is it possible, once I have entered the metric, to numerically integrate a specific geodesic and then plot it?

My question is about the physics package: when transforming the metric tensor to new coordinates, sometimes the new coordinate can be expressed easily as a variable inside the new matrix entries,  but the new matrix is still displayed completely in terms of the old variables. How can I force it to display everything in terms of the new ones?

How do I calculate the conformal killing vectors of a metric tensor in GRTensorIII?