def dr_central(f, dr):
    """
    First order central difference differentiation.

    Parameters
    ----------
    f: ndarray
        The variable to be differentiated.
    dr: scalar
        Step size in r.

    Returns
    -------
    ndarray

    """

    return (f[2:, 1:-1] - f[:-2,1:-1])/(2*dr)