function ShowCSinterp()
% Generate some periodic data with some noise...
n = 50;
m = n/2;
h = 2*pi/n;
t = linspace(0,2*pi-h,n)';
fVals = 5*cos(9.1*t) - 6*sin(4.3*t) + .3*randn(n,1);
plot(t,fVals,'o')
h0 = title('The Data','fontsize',14);
shg
pause(3)

% Fit the data with a sum of cosine and sine functions...
% F = CSinterpSlow(fVals);
F = CSinterpFast(fVals);

% Display the interpolant...
z = linspace(0,2*pi,1000)';
y = CSeval(F,z);
hold on
delete(h0)
title('The Interpolant','fontsize',14)
plot(z,y,'r')
hold off
shg