In most physical systems the dynamical behaviour is best described using differential equations and physical considerations often suggest that a discrete observation model is appropriate. This paper gives a numerical method for the recursive estimation of the parameters of a linear time-invariant differential equation model from discrete observations. The numerical method in an alternate formulation gives a simple way of solving the differential equation.