Ndim_ :=    4   :
x1_   :=   x   :
x2_   :=   y   :
x3_   :=   z   :
x4_   :=   t   :
g11_   :=   L^2/(x^2+y^2+z^2)   :
g22_   :=   L^2/(x^2+y^2+z^2)   :
g33_   :=   L^2/(x^2+y^2+z^2)   :
g44_   :=   -P(x,y,z,t)^2   :
constraint_ :=   [P(x,y,z,t) = a(t)*cos(ln((x^2+y^2+z^2)^(1/2)/L))+b(t)*sin(ln((x^2+y^2+z^2)^(1/2)/L))]   :