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