code: def accel_fl(q: np.ndarray):
      
      # 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
docs: Accelaration in the earth-sun system using Fluxion potential energy