Accelaration in the earth-sun system using Fluxion potential energy

def accel_ES_fl(q: np.ndarray):
    """Accelaration in the earth-sun system using Fluxion potential energy"""
    # Number of celestial bodies
    num_bodies: int = 2
    # Number of dimensions in arrays; 3 spatial dimensions times the number of bodies
    dims = 3 * num_bodies

    # The force given the positions q of the bodies
    f = force_ES(q)

    # The accelerations from this force
    a = np.zeros(dims)
    a[slices[0]] = f[slices[0]] / mass[0]
    a[slices[1]] = f[slices[1]] / mass[1]
    
    return a