Routines for numerical interpolation, with 1st and 2nd order differentiation, having non-uniformly spaced points, out to 3 dimensions - CPC, Mendeley (Mendeley, 1987)