% plot all the atoms

  hold off
  for iatom=i1:i2
     for jatom=top(iatom):-2:1
        plotatom(iatom,jatom,'b')
        hold on
	end
     end

  axis([i1 i2 0 (i2-i1)])