Ndim_ := 4:
x1_ := xi:
x2_ := theta:
x3_ := phi:
x4_ := t:
sig_ := 2:
complex_ := {}:
g11_ := R^2*(exp(c-k+k*xi^2)+exp(-c+k-k*xi^2)+2)^2/(exp(c-k+k*xi^2)+exp(-c+k-k*xi^2))^2:
g22_ := R^2*xi^2*(exp(c-k+k*xi^2)+exp(-c+k-k*xi^2)+2)^2/(exp(c-k+k*xi^2)+exp(-c+k-k*xi^2))^2:
g33_ := R^2*xi^2*sin(theta)^2*(exp(c-k+k*xi^2)+exp(-c+k-k*xi^2)+2)^2/(exp(c-k+k*xi^2)+exp(-c+k-k*xi^2))^2:
g44_ := -(exp(c-k+k*xi^2)+exp(-c+k-k*xi^2)-2)/(exp(c-k+k*xi^2)+exp(-c+k-k*xi^2)+2):
constraint_ := [c=ln(sqrt(alpha*(alpha-2))+alpha-1+sqrt(2*alpha^2-4*alpha+2*sqrt(alpha*(alpha-2))*alpha-2*sqrt(alpha*(alpha-2)))),k=1/8/alpha*(2*alpha-2*sqrt(alpha^2-2*alpha))*(1+(sqrt(alpha*(alpha-2))+alpha-1+sqrt(2*alpha^2-4*alpha+2*sqrt(alpha*(alpha-2))*alpha-2*sqrt(alpha*(alpha-2))))^2)*(sqrt(alpha*(alpha-2))+alpha+sqrt(2*alpha^2-4*alpha+2*sqrt(alpha*(alpha-2))*alpha-2*sqrt(alpha*(alpha-2))))/(sqrt(alpha*(alpha-2))+alpha-1+sqrt(2*alpha^2-4*alpha+2*sqrt(alpha*(alpha-2))*alpha-2*sqrt(alpha*(alpha-2))))/(-2+sqrt(alpha*(alpha-2))+alpha+sqrt(2*alpha^2-4*alpha+2*sqrt(alpha*(alpha-2))*alpha-2*sqrt(alpha*(alpha-2))))]:
Info_ := `S.P. Goldman, ApJ., 226, 1079 (1978)`: