fixing a bug in interp1.m