load data.mat
ydata(end) = [];
slope0 = struct('p', 1, 'x', [min(xdata) max(xdata)], 'v', [0 0]);
pp = BSFK(xdata,ydata, 2, 3, [], struct('KnotRemoval', 'none', 'pntcon', slope0));
figure
xi = linspace(min(xdata),max(xdata));
yi = ppval(pp,xi)
plot(xdata,ydata,'b',xi,yi,'r');
for xb=pp.breaks
xline(xb);
end
grid on
Best Answer