I'm going to leave this alone for a while, because positivity loss is still a problem with QuTiP's solvers. My own Euler-Maruyama seems to be stable, I'm going to try Milstein but that's it