Numerical techniques for Black hole solutions in f(R) gravity