A new algorithm for Monte Carlo simulation of Ising spin systems