x=r.cos D cos A
y=r.cos D sen A
z=r.sen D
cos d=senD1 senD2+cosD1 cosD2 cos (A1-A2)
Un applet de Java permite calcular la distancia angular
Análoga expresión es útil si las coordenadas son eclípticas, los correspondientes papeles los desarrollan la longitud y latitud.
También se puede hallar el ángulo a partir de los cosenos directores (a,b,c) de cada astro
cos d= (a1,b1,c1).(a2,b2,c2) siendo . el producto escalar.
Los cosenos directores de cada astro son:
a=x/r=cos D cos A
b=y/r=cos D sen A
c=z/r=sen D


El programa calcula el ángulo en ° ' " que forman dos astros cuyas coordenadas se dan en ecuatoriales. Establece la diferencia de calcular el ángulo entre dos astros móviles o entre una estrella fija y un cuerpo que se mueve . En este último caso se prescinde de ingresar por segunda vez las coordenadas del astro fijo. Puede ser útil para calcular las conjunciones planetarias o de un planeta y una estrella o incluso las ocultaciones de una estrella por la Luna si se dispone de las coordenadas topocéntricas lunares. El programa se ha probado con los datos del Practical Astronomy (Minaert) : &31 The angle between two celestial objets , &9 Angular separations.
10 REM CALCULO ANGULO ENTRE DOS ASTROS 20 CLS:SCREEN 0:PI=3.1415928# 30 PRINT "ANGULO ENTRE DOS ASTROS" 40 DEF FN AC(X)=-ATN (X/SQR(-X*X+1))+PI/2 50 PRINT "ELIJA ENTRE OPCIONES : 1.DOS ASTROS MOVILES" 60 PRINT " 2.PRIMER ASTRO FIJO Y EL OTRO MOVIL" 70 INPUT X 80 PRINT :PRINT "INGRESE COORDENADAS ECUATORIALES DEL PRIMER ASTRO" 90 GOSUB 230 100 F=A:G=B:H=C 110 PRINT :PRINT "INGRESE COORDENADAS ECUATORIALES DEL SEGUNDO ASTRO" 120 GOSUB 230 130 A1=FNAC(F*A+G*B+H*C)*180/PI 140 GR=INT (A1) 150 M=INT((A1-INT(A1))*60) 160 S=((A1-INT(A1))*60-M)*60 170 PRINT :PRINT "EL ANGULO ENTRE LOS DOS ASTROS ES: "; 180 PRINT GR;" G ";M;" M ";S;" S " 190 INPUT "OTRO ANGULO (S/N) ";A$ 200 IF A$="N" OR A$="n" THEN 225 210 ON X GOTO 80,110 225 END 230 INPUT "ASCENSION RECTA (H,M,S)=";HH,MM,SS 240 E=(HH+MM/60+SS/3600)*PI/12 250 INPUT "DECLINACION (G,M,S)=";GG,MM,SS 260 IF GG<0 THEN D=-(-GG+MM/60+SS/3600)*PI/180 270 IF GG>=0 THEN D=(GG+MM/60+SS/3600)*PI/180 280 A=COS(D)*COS(E) 290 B=COS(D)*SIN(E) 300 C=SIN(D) 310 RETURN