% Generate a spline of raw data and plot it for a
% larger number of x coordinates.  

clear
close all

load intdata2.dat
x = intdata2(:,1);
y = intdata2(:,2);
% Compute a not-a-knot spline for the data
S = spline(x,y);

xo = linspace(x(1),x(end),101)';
yo = ppval(xo,S);

figure(1)
    hold on
    plot(x,y,'b*')
    plot(xo,yo,'ro')
    xlabel('x')
    ylabel('y')
    title('Not-a-knot Spline Interpolation')