Helper function to find differences between two equilibria

Helper function to find differences between two equilibria

compare_eq(eq1, eq2 = g[["eq"]], g, verbose = TRUE)