Some ideas on speed up: - Common subexpression elimination - numba compilation https://stackoverflow.com/questions/30738840/best-practice-for-using-common-subexpression-elimination-with-lambdify-in-sympy