I have a set of coupled differential equations of the form $$ \omega y_1 \frac{\partial q_1}{\partial y_1}+\omega y_2 \frac{\partial q_1}{\partial y_2}+\omega y_3 \frac{\partial q_1}{\partial y_3}+g\frac{\partial q_1}{\partial y_1}+t q_{3}e^{-g/\omega (y_{3}-y_1)}-t q_{2}e^{-g/\omega (y_{2}-y_3)}...