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
def accel_fl ( q : np . ndarray ) : dims : int = len ( q ) B : int = dims // 3 f = force ( q ) 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