
% run after running SLIT_FN

%z = polyfit([(m - centre + 1):(m - centre + max_len)],s_fn(m,1:max_len),30);
z = polyfit([(m - 15):(m + 20)],s_fn(m,centre-15:centre+20),30);
%z1 = polyval(z,[(m - centre + 1):(m - centre + max_len)]);
z1 = polyval(z,[(m - 15):(m + 20)]);
figure(2)
clf
plot(z1,'k');
hold on
plot(s_fn(m,1:max_len),'g');
hold off;
