Some simple models (such as integrate and fire with current-based synapses) can be solved explicitly and don't need a general ODE solver. Sympy's dsolve function could be used to check whether a model can be solved for and then give an explicit solution to the model thereby improving the simulation time (probably should be part of lib9ML actually)