% Nx=1700;
% Ny=1500;
eNf=Nf-1;
sNf=1;

for nf=sNf:eNf
  hold off
  plot(px(1:Npf(nf+1),nf+1),py(1:Npf(nf+1),nf+1),'r.','MarkerSize',15);
  hold on                     
  plot(px(1:Npf(nf),nf),py(1:Npf(nf),nf),'b.','MarkerSize',15);
%  text(px(1:Npf(nf+1),nf+1),py(1:Npf(nf+1),nf+1),num2str((1:Npf(nf))'),'horizontalAlignment','center');
%  text(px(1:Npf(nf),nf),py(1:Npf(nf),nf),num2str((1:Npf(nf))'),'horizontalAlignment','center');
%   axis ij;                 
  axis equal;
  axis([-.1*Nx Nx*1.1 -.1*Ny Ny*1.1]);
  mm=find(cc(1,1:Npf(nf),nf)>0);
  plot([px(mm,nf) px(cc(1,mm,nf),nf+1)]',...
       [py(mm,nf) py(cc(1,mm,nf),nf+1)]','k','linewidth',3)
  disp(nf);
	if(length(mm)<256)
		pause;
	end
	drawnow;
end

