This paper presents a recursive algorithm using orthogonal Laguerre polynomials to identify the parameters of a system described by a linear differential equation with coefficients of finite-order polynomials in timet. The original time-varying differential input-output model is converted to a set of linear algebraic equations which are solved recursively using Greville's algorithm to give unknown parameters and initial conditions. Numerical results are given to illustrate the algorithm.