functional astrometry routines!

I now have functional routines that (a) find the best reverse polynomial distortion terms given forward polynomial distortion terms, and (b) invert any x,y to RA,Dec transformation by Newton's method.