I have a system of stochastic differential equations that I would like to solve. I was hoping that this issue was already address. I am a bit concerned about constructing my own solver because I fear my solver would be too slow, and there could be the issues with numerical stability. 

Is there a python module for such problems?

If not, is there a standard approach for solving such systems.