% Data used to run the subroutine SubSpline 
clear;

%x=[1.2 1.3 1.4 1.5];
%a=[2.57 3.6 5.79 14.10];
%N=4;

% This data correspond to the ruddy duck of Burden-Faires Section 3.4

x=[0.9 1.3 1.9 2.1 2.6 3 3.9 4.4 4.7 5 6 7 8 9.2 10.5 11.3 11.6 12 ...
    12.6 13 13.3];
a=[1.3 1.5 1.85 2.1 2.6 2.7 2.4 2.15 2.05 2.1 2.25 2.3 2.25 1.95 ...
      1.4 0.9 0.7 0.6 0.5 0.4 0.25];
N=21;
Subspline(N,x,a);