Compute offsets between to RA,DEC positions in sexigecimal format. The first (lmfg) version is just straight gfortran interactive code:

%  RA_DEC_SEP 

                    --:--:--.-- ---:--:--.- 
Enter center RA,DEC:02:19:00.01 +57:07:41.9 

                    --:--:--.-- ---:--:--.- 
Enter center RA,DEC:02:20:25.32 +57:34:18.6 

Position 1:  02:19:00.01 +57:07:41.9
Position 2:  02:20:25.32 +57:34:18.6
a(x,y) arcmin              =     -11.4370     26.6423
Separation in Degrees      =       0.483223
Separation in arc_MINUTES  =      28.993364
Separation in arc_SECONDS  =    1739.601929