A method of interpolation is given in which the function and its first derivative are continuous at all points, but which is simpler than spline interpolation. It is applicable to functions of any number of variables, when the prescribed values are evenly or irregularly distributed.