178 void palEl2ue (
double date,
int jform,
double epoch,
double orbinc,
179 double anode,
double perih,
double aorq,
double e,
180 double aorl,
double dm,
double u[13],
int *jstat ) {
183 const double SE=0.3977771559319137;
184 const double CE=0.9174820620691818;
188 double PHT,ARGPH,Q,W,CM,ALPHA,PHS,SW,CW,SI,CI,SO,CO,
189 X,Y,Z,PX,PY,PZ,VX,VY,VZ,DT,FC,FP,PSI,
193 if (jform < 1 || jform > 3) {
197 if (e < 0.0 || e > 10.0 || (e >= 1.0 && jform != 3)) {
205 if (jform == 1 && dm <= 0.0) {
222 PHT = epoch-(aorl-perih)/dm;
226 CM = W*W*aorq*aorq*aorq;
228 }
else if (jform == 2) {
231 PHT = epoch-aorl*sqrt(aorq*aorq*aorq)/
PAL__GCON;
250 ALPHA = CM*(e-1.0)/Q;
254 PHS = sqrt(ALPHA+2.0*CM/Q);
315 W = pow(3.0*DT+sqrt(9.0*DT*DT+8.0*Q*Q*Q), 1.0/3.0);
317 PSI = (1.0-e)*FC+e*FP;
342 palPv2ue( PV, date, CM-1.0, u, &J );
void palPv2ue(const double pv[6], double date, double pmass, double u[13], int *jstat)
void palEl2ue(double date, int jform, double epoch, double orbinc, double anode, double perih, double aorq, double e, double aorl, double dm, double u[13], int *jstat)
static const double PAL__GCON
void palUe2pv(double date, double u[13], double pv[], int *jstat)