Application of the piecewise linear polynomial functions expansion is extended to linear time-varying systems. With the treatment of the product of two time functions, two types of operational matrix are developed. By applying these operational matrices, the dynamic equations are transformed into a set of algebraic equations. A recursive algorithm is derived and the system equations can be solved with very low dimensional matrix inversions. This represents a considerable saving of computer memory capacity and computing time.