def accel_fl(q: np.ndarray):
    """Accelaration in the earth-sun system using Fluxion potential energy"""
    # Infer number of dimensions from q
    dims: int = len(q)
    # Number of celestial bodies
    B: int = dims // 3

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

    # The accelerations from this force
    a = np.zeros(dims)
    for i in range(B):
        a[slices[i]] = f[slices[i]] / mass[i]
    
    return a
Accelaration in the earth-sun system using Fluxion potential energy