clear clc close all y = [1 2 2.2 3 2.8 4.5 5 8 9]; x = 1:length(y); plot(x, y,'r+') x1 = 1:.1:length(y); y1 = interp1(x, y, x1); hold on plot(x1, y1,'.') sply = spline(x, y, x1) plot(x1, sply, 'k')