  The partition function of the 2d Ising model with random nearest neighbor
coupling is expressed in the dual lattice made of square plaquettes. The dual
model is solved in the the mean field and in different types of Bethe-Peierls
approximations, using the replica method.
