Multidimensional interpolation using osculatory radial basis functions