Create GibbsIsing4.m

Gibbs method for rectangular net