% Guess_PCF_Demo.m

load NonCompact.dat;
qrange = [ 0.02 0.17]; % Range in which data is good.
D = 250; % Maximum correlation length
N = 20; % Initially make p(r) up of 20 segments

% Start a rough fit
[r, p, res_err] = Guess_PCF(NonCompact, qrange, D, N);

for k=1:3
% Now round out the rough edges from the first go.
p = smooth(p); p(1)=0; p(length(p))=0; p = p';
% Fit again
size(r), size(p)
[r, p, res_err] = Guess_PCF(NonCompact, qrange, D, N,r,p);

end

p = smooth(p); p(1)=0; p(length(p))=0; p=p';
figure(2); plot(r,p,'g');