H = -\sum_{} J_{ij} \sigma_{i} \sigma_{j}
P(J_{ij}) = p \delta( J_{ij }+J ) + (1-p ) \delta (J_{ij}-J)
