Ndim_ := 4:
x1_ := xi:
x2_ := theta:
x3_ := phi:
x4_ := t:
sig_ := 2:
complex_ := {}:
g11_ := -112*R^2*(1+C*xi^2)^3*(1+6*C*xi^2)^(1/3)/(-112*(1+6*C*xi^2)^(1/3)+309*C*xi^2*(1+6*C*xi^2)^(1/3)+54*C^2*xi^4*(1+6*C*xi^2)^(1/3)+8*C^3*xi^6*(1+6*C*xi^2)^(1/3)-112*K*C*xi^2):
g22_ := xi^2*R^2:
g33_ := xi^2*R^2*sin(theta)^2:
g44_ := -(1+C*xi^2)^5:
constraint_:=[C=1/((-11+5*alpha)),K=-125/112*5^(1/3)*((-1+alpha)*(-11+5*alpha)^2)^(1/3)*(-1792+2097*alpha-792*alpha^2+95*alpha^3)/(-11+5*alpha)^3/alpha]:
Info_:=`M. Durgapal, J. Phys. A, 15,  2637 (1982), M. Korkina, Sov. Phys. J. 24,  468 (1981)`: