The Ising model, in the presence of a random field, is investigated within the mean-field approximation based on Landau expansion. The random field is drawn from the bimodal probability distribution P(h)=pδ(h−h0)+(1−p)δ(h+h0), where the probability p assumes any value within the interval...