IsingPL function

Pseudo-likelihood