Non-linear system of two ODEs of 1st order

This PR contain the methods for nonlinear system of ODEs for 2 equations. This also includes all test cases for these methods and also test cases for classify_sysode.
Is there a case where you wouldn't want to compare `len(set(func))` to `len(eq)`? I'd figure you'd only want unique values for `func`.