I doubt there is an analytical solution for your system because of the S and DS-terms in the denominators.
Use one of the numerical integrators instead, e.g. ODE45, ODE15S.
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!