We present a new numerical method to study critical properties of disordered systems. For the two-dimensional diluted Ising model, the critical line and exponents are determined very accurately and the computational effort required is moderate even for impurity concentrations down to 0.60. The method can be generally applicable to other disordered models exhibiting second-order phase transition.