On Monday, June 1, 2020 at 4:56:29 PM UTC+3, wrote:
I notice that ta has typo error. I have written this long formula by hand with my mobile phone. I try to post ta again when I get proper internet connection which I did not have when I posted these. Our libraries are all closed their internet connection which I have used before due this coronavirus pandemia. Sorry that typo error. Hannu

CORRECT ONE:

# I give this complicated primitive function here (Maple 9 used): Analytic solution (Weinberg S. Gravitation and Cosmology, 1972, General Relativity definitions used from this book).

# Defintion areas: u=1.843608337*10^(-15), 0r=5.424145573*10^(-14) cm, u=1/r,0=P=Pi.

ta := P- 0.14273426e18*arctan(19556.51583*(0.4000000000e15-1506485618*sin(P)^2)^(1/2)*sin(P)/(-0.2000000000e15+753242809*sin(P)^2))-0.1021331259e48*(0.4000000000e15-1506485618*sin(P)^2)^(1/2)*sin(P)/((-0.2000000000e15+753242809*sin(P)^2)*(0.3774722128e 26+0.1443670077e35*(0.4000000000e15-1506485618*sin(P)^2)*sin(P)^2/(-0.2000000000e15+753242809*sin(P)^2)^2))+0.19991518 56e-30*((-0.2000000000e15+753242809*sin(P)^2)*(-1+sin(P)^2))^(1/2)*(0.3908757305e66*(1-sin(P)^2)^(1/2)*(0.4000000000e15-1506485618*sin(P)^2)^(1/2)*EllipticF(sin(P), 0.1940673606e-2)*sin(P)^2+0.6697771796e73*(1-sin(P)^2)^(1/2)*(0.4000000000e15-1506485618*sin(P)^2)^(1/2)*EllipticF(sin(P), 0.1940673606e-2)+0.2468340123e75*sin(P)+0.9296297238e69*sin(P)^5-0.2468349419e75*sin(P)^3+0.1234170061e68*(1-sin(P)^2)^(1/2)*(0.4000000000e15-1506485618*sin(P)^2)^(1/2)*EllipticE(sin(P), 0.1940673606e-2)*sin(P)^2+0.2114787075e75*(1-sin(P)^2)^(1/2)*(0.4000000000e15-1506485618*sin(P)^2)^(1/2)*EllipticE(sin(P), 0.1940673606e-2)-0.1273248266e68*(1-sin(P)^2)^(1/2)*(0.4000000000e15-1506485618*sin(P)^2)^(1/2)*EllipticPi(sin(P), -0.5835906962e-7, 0.1940673606e-2)*sin(P)^2-0.2181748740e75*(1-sin(P)^2)^(1/2)*(0.4000000000e15-1506485618*sin(P)^2)^(1/2)*EllipticPi(sin(P), -0.5835906962e-7, 0.1940673606e-2))/((0.2000000000e15+753242809*sin(P)^4-0.2000007532e15*sin(P)^2)^(1/2)*(0.2202892715e19*sin(P)^2+0.3774722128e26)*cos( P)*(0.4000000000e15-1506485618*sin(P)^2)^(1/2));

u:=P-(-8.905287066*10^(-19)*cos(P)+3.686326146*10^(-15)) / (-cos(P) +1);