Ndim_ :=    4   :
x1_   :=   x   :
x2_   :=   y   :
x3_   :=   z   :
x4_   :=   t   :
eta11_   :=   -1   :
eta22_   :=   -1   :
eta33_   :=   -1   :
eta44_   :=   1   :
bd11_   :=   L/(x^2+y^2+z^2)^(1/2)   :
bd22_   :=   L/(x^2+y^2+z^2)^(1/2)   :
bd33_   :=   L/(x^2+y^2+z^2)^(1/2)   :
bd44_   :=   P(x,y,z,t)   :
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))]   :