55   OrbitType *OrbPtr;         /* orbit pointer */  in FindPeriastron()  local
57   OrbPtr = &Orbit;  in FindPeriastron()
59   OrbPtr->Nu = (3.0 * PI / 2.0) - (DTOR * OrbPtr->Omega);  in FindPeriastron()
62   temp = sqrt( (1.0 + OrbPtr->Excentricity) / (1.0 - OrbPtr->Excentricity));  in FindPeriastron()
63   temp = tan(OrbPtr->Nu/2.0) / temp;  in FindPeriastron()
64   OrbPtr->E  = 2.0 * atan(temp);   in FindPeriastron()
65   if (OrbPtr->E <= (-FLT_EPSILON)) OrbPtr->E = OrbPtr->E + 2.0 * PI;  in FindPeriastron()
68   OrbPtr->OmegaZero = -(OrbPtr->E - OrbPtr->Excentricity * sin(OrbPtr->E));  in FindPeriastron()
71   OrbPtr->MinDist = (1.0 - SQR(OrbPtr->Excentricity))/( 1.0 + OrbPtr->Excentricity);  in FindPeriastron()
74   OrbPtr->Nu = PI - (DTOR * OrbPtr->Omega);   in FindPeriastron()
75   temp = sqrt( (1.0 + OrbPtr->Excentricity) / (1.0 - OrbPtr->Excentricity));  in FindPeriastron()
76   temp = tan(OrbPtr->Nu/2.0) / temp;  in FindPeriastron()
77   OrbPtr->E  = 2.0 * atan(temp);  in FindPeriastron()
79   OrbPtr->MQuad = OrbPtr->E - OrbPtr->Excentricity * sin(OrbPtr->E);  in FindPeriastron()
85     printf(_(" Mean Anomaly at Secondary Eclipse %7.4f,\n"),OrbPtr->OmegaZero);  in FindPeriastron()
86     printf(_(" Minimum Distance %6.3f,\n"),OrbPtr->MinDist);  in FindPeriastron()
87     printf(_(" Mean Anomaly at Quadrature %7.4f\n"),OrbPtr->MQuad);   in FindPeriastron()
106   OrbitType *OrbPtr;                           /* orbit pointer          */  in SolveKepler()  local
108   OrbPtr = &Orbit;  in SolveKepler()
116  OrbPtr->M = ( OrbPtr->MQuad + (2.0*PI/PhaseSteps)*j );   in SolveKepler()
121  New =  OrbPtr->M + OrbPtr->Excentricity * sin(OrbPtr->M);  in SolveKepler()
125    New = OrbPtr->M + OrbPtr->Excentricity * sin(Old);  in SolveKepler()
128  OrbPtr->E = New;  in SolveKepler()
130  OrbPtr->Nu = 2.0 * atan( tan(OrbPtr->E/2.0)   in SolveKepler()
131                    * ( sqrt ( (1.0 + OrbPtr->Excentricity)   in SolveKepler()
132                      / (1.0 - OrbPtr->Excentricity))));  in SolveKepler()
135  OrbPtr->Phase = OrbPtr->Nu + (DTOR * OrbPtr->Omega) - (PI/2.0);  in SolveKepler()
136  if (OrbPtr->Phase <= (-FLT_EPSILON)) OrbPtr->Phase = OrbPtr->Phase + 2.0*PI;  in SolveKepler()
142  OrbPtr->Dist = (1.0 - SQR(OrbPtr->Excentricity) )   in SolveKepler()
143           / ( 1.0 + OrbPtr->Excentricity * cos (OrbPtr->Nu));  in SolveKepler()
146  OrbPtr->Dist = OrbPtr->Dist / OrbPtr->MinDist;  in SolveKepler()
147  if ((OrbPtr->Dist - 1.0) <= (-FLT_EPSILON)) OrbPtr->Dist = 1.0;   in SolveKepler()
176    OrbitType *OrbPtr = &Orbit;       /* orbit pointer          */  in UpdateGeometry()  local
181         Binary[Comp].Volume / (OrbPtr->Dist * OrbPtr->Dist * OrbPtr->Dist);  in UpdateGeometry()