%% script: mcldemo


load mclExamp.mat

rs = [1.3, 1.35, 1.4, 1.5, 1.53, 1.55, 1.57, 1.6];

for r = rs,
    
m = mcl(a,r,1);

viewmcl(a,xy,m)
spymcl(m)
figure(1);
title(['r = ', num2str(r)]);

pause

end
