Utilizando las relaciones :
X=D cos l cos b ; Y = D sin l cos b ; Z = D sin b
lo que a través del applet se traduce, por ejemplo para X , por:
texttt= G1 * cos(C1/180.*pi) * cos(D1/180.*pi)... encontramos la solución

tpshapley3.png
Crédito : Astrophysique sur Mesure