function ising_plot ( x, y, spin )

up = find ( spin == +1 );
down = find ( spin == -1);

plot ( x(down), y(down), 'ro', x(up), y(up), 'bo' )
axis square
