n = 10; %n = input('n = ');
x = (1:n)';
y = randn(n,1);
u = (1:.01:n)';
v1 = splinetx(x,y,u);
v2 = pchiptx(x,y,u);
v3 = yzspline(x,y,u);
v4 = myspline(x,y,u);
plot(x,y,'r*',u,v1,':',u,v2,'-.',...
    u,v3,'k-',u,v4,'k-','linewidth',2);