  We apply the Wang-Landau method to the study of the critical behaviour of the
three dimensional Random Field Ising Model with a bimodal probability
distribution. Our results show that for high values of the random field
intensity the transition is first order, characterized by a double-peaked
energy probability distribution at the transition temperature. On the other
hand, the transition looks continuous for low values of the field intensity. In
spite of the large sample to sample fluctuations observed, the double peak in
the probability distribution is always present for high fields
