1d data interpolation

x = 0:0.5:10;   y = sin(x);
plot(x,y,'ko-');
hold on
xi = 0:.25:10;
for i=1:size(xi,2)
    yi(i) = interp1(x,y,xi(i)); % default : linear
    % yi(i) = interp1(x,y,xi(i),'spline');
    plot(xi(i),yi(i),'rx');
    pause
end