The attached paper refers to computer programs for SI modelling. You may try to contact the author. There are plenty of programs for solution of systems of differential equations. which would give incidence rates.
There are many differential equations solver packages in R. You just need to set up the equations and the solver will do the job for you. It basically takes the form: