It is not very difficult to integrate this DE numerically. For parameter estimation it is a good idea for stability to use a semi-implicit formulation. The idea is described here. http://otter-rsch.com/admodel/cc4.html