Numerical interpolation MATLAB code