A Monte Carlo method has been used to study a simple cubic Ising ferromagnet in a random quenched magnetic field. The Hamiltonian for this model is ℋ==JΣ(ij)σiσj−ΣiHiσi, where σi,σj=±1, J is the nearest-neighbor interaction constant, and the field Hi=tH is fixed at each site with ti=±1 at random and Σti=0. L×L×L lattices with periodic boundary conditions have been studied for a range of H and T. As expected we find a ferromagnetic ordered state which for small H undergoes a second order phase change to the paramagnetic state with increasing temperature. A finite size scaling analysis of the preliminary data suggests that the critical exponent β is substantially smaller (β∼0.2) than the usual 3-dim Ising value of 0.31. Results obtained for small lattices indicate that below kT/J-2 the transition becomes first order suggesting that a tricritical point appears on the critcal field curve.