I give here those two analytic solutions:First primitive function is following:e=2.117198112*EllipticF(sin(P),0.2601014 318), r=(0.1627132736e34*sin(P)^2-0.2405124138e35)/(0.2726824096e18*sin(P)^2-0.3378949220e18), definition areas are 0=P=Pi/2 and 0.7117964733e17=r=0.3438620569e18. Second primitive function is following: ee=2.117198112*EllipticF(sin(P),0.9655812990), rr=(0.2051871297*sin(P)^2)/(0.5967134948e16*sin(P)^2+0.3378949220e18), definition areas are 0=P=Pi/2 and 0=rr=0.5967134948e16. Maple 9 used and cgs unit system with c=1 Weinberg's used units. Best Regards Hannu Poropudas, Finland