Monte Carlo simulation of three-dimensional dilute Ising model