# Open source solver for continuous-time stochastic non-linear DAEs (SDAEs)

I am trying to solve a system of non-linear index-1 DAEs in which the derivatives of the state variables, $x(t)$ are corrupted by additive noise, $w(t)$ (whose co-variance matrix is known).

$\dot x(t) = f(x(t),t,u(t)) + w(t)$

Is there a known (open source) solver/library/package/toolbox in MATLAB/python/Julia (with or without adaptive time-stepping capability) that can solve the system?

I will caution that this is pretty new so I would test it on a non-DAEs, i.e. using the mass-matrix functionality and then adding $M^{-1}$ to your SDE equations and checking to make sure that they calculate the same thing. If you notice any discrepancies, please file an issue. The architecture for the implicit handling here is based on the ODE solvers which are very well-tested, but since I see we don't explicitly has a mass matrix SDE test in the continuous integration suite yet I would just do this test before trusting it.