Sympy code to collect terms for various PN equations, and generate working c++ code to evolve PN systems