Ndim_ := 4: x1_ := xi: x2_ := theta: x3_ := phi: x4_ := t: sig_ := 2: complex_ := {}: g11_ := 2*R^2*(3*alpha-7+xi^2)*((3*alpha-7+4*xi^2)*(3*alpha-7))^(1/2)*alpha/(6*((3*alpha-7+4*xi^2)*(3*alpha-7))^(1/2)*alpha^2-14*((3*alpha-7+4*xi^2)*(3*alpha-7))^(1/2)*alpha-9*xi^2*((alpha-1)*(3*alpha-7))^(1/2)*3^(1/2)*alpha-((3*alpha-7+4*xi^2)*(3*alpha-7))^(1/2)*xi^2*alpha+24*xi^2*((alpha-1)*(3*alpha-7))^(1/2)*3^(1/2)): g22_ := R^2*xi^2: g33_ := R^2*xi^2*sin(theta)^2: g44_ := -(3*alpha-7+xi^2)^3/(3*alpha-7)^3: constraint_:=[]: Info_:=`H. Heintzmann, Z. Phys., 228, 489 (1969)`: