code: def calc_mse(q1, q2):
      
      # Difference in positions between two simulations
      dq = q2 - q1
      # Mean squared error in AUs
      return np.sqrt(np.mean(dq*dq))/au2m
docs: Compare the results of two simulations