Müller & Janka [208Jump To The Next Citation Point] did not specifically describe what they used, but any standard derivative discretization would work, e.g., ∂Φ ∕∂r = (Φk+1 − Φk −1)∕(rk+1 − rk−1).