%% prdemo

load nov7
a = dista(a,xy);

d = diag(1./sum(a));
m = d*a;

al = 0.1;

s = zeros(1,100); s(1) = 1;
v = s;
si = s;
for i = 1:40,
si = si * m;
v = v + (1-al)^i * si;
end

[x,y] = gplot(a,xy);
plot(x,y,'-','Color',[.9 .9 .9])

[jnk,ord] = sort(-v);

for i = 1:100,
    

  plot(xy(ord(i),1),xy(ord(i),2),'ro');

  pause
end
