uiopen('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/orderspline.m', true);
% Modify expression to add input arguments.
% Example:
%   a = [1 2 3; 4 5 6]; 
%   foo(a);

orderspline
orderspline
orderspline
x=linspace(0,30,11);
y=atan(x);
xx=linspace(0,30,100);
yy=quadraticspline([x(1),x],[1,y],xx);
plot(x,y,'*',xx,yy)
plot(x,y,'*',xx,yy,xx,atan(xx))
yy=quadraticspline([x,x(end)],[y,1/901],xx);
plot(x,y,'*',xx,yy,xx,atan(xx))
diary off