function yi=interp(x,y,xi,method)
% 一維非單調數據插值
% 用法參考interp1，唯一區別是interp不要求節點是單調的

% L.J.HU 4-15-1999

x=x(:);y=y(:);
z=[x,y];
z=sortrows(z);
yi=interp1(z(:,1),z(:,2),xi,method);