ymx=ceil(max(tpy(:)))+5;
ymn=floor(min(tpy(:)))-5;
NNy=ymx-ymn+1;
xmx=ceil(max(tpx(:)))+5;
xmn=floor(min(tpx(:)))-5;
NNx=xmx-xmn+1;

ms=5;

eNf=Nf-1;
eNf=46;
sNf=45;

for nf=sNf:eNf
  hold off
  plot(tpx(1:Npf(nf+1),nf+1),tpy(1:Npf(nf+1),nf+1),'ro','MarkerSize',ms);
  hold on                     
  plot(tpx(1:Npf(nf),nf),tpy(1:Npf(nf),nf),'bo','MarkerSize',ms);
	if(0)
		text(tpx(1:Npf(nf+1),nf+1),tpy(1:Npf(nf+1),nf+1),num2str((1:Npf(nf))'),'horizontalAlignment','center');
		text(tpx(1:Npf(nf),nf),tpy(1:Npf(nf),nf),num2str((1:Npf(nf))'),'horizontalAlignment','center');
	end
  %axis ij;                 
  axis equal;
  axis([xmn xmx ymn ymx]+.05*[NNx*[-1 1] NNy*[-1 1]]);
  mm=find(cc(1:Npf(nf),nf)>0);
  plot([tpx(mm,nf) tpx(cc(mm,nf),nf+1)]',...
       [tpy(mm,nf) tpy(cc(mm,nf),nf+1)]','k','linewidth',3)
  disp(nf);
  pause(.2);
end

