I have the following differential equation array, $$ \frac{\mathrm d W}{\mathrm d t} = [1-aq(t-R_0) +b]{c\over q(t)+d } - [aq(t-R_0) - b]{W(t)\over q(t)+d}, ~~~~~~~~~~~~~~~~~~~~~~~(1)$$
$$ ~~~~~~~~~~~~~~~~\frac{\mathrm d q}{\mathrm d t} = {W(t)\over q(t)+d } -e, ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~(2)$$
where $a, b, c, d, e, R_0$ are constants. $q(t-R_0)$ denotes $q(t)$ with time delay $R_0$.
How can I solve it? I don't expect any closed-form expression for $W(t)$ and $q(t)$, numerical solutions with tools like matlab are welcome.