  We perform a duality transformation that allows one to express the partition
function of the d-dimensional Ising model with random nearest neighbor coupling
in terms of new spin variables defined on the square plaquettes of the lattice.
The dual model is solved in the mean field approximation.
