A method for finding the solution of a linear two-point boundary value problem with time-varying coefficients is discussed. Properties of shifted Chebyshev series are first briefly presented and the transformation matrix relating the back vector to the current time vector together with the operational matrix are utilized to solve the two-point boundary value problems. This approach can be applied to obtain the optimal control of linear time-varying systems subject to quadratic cost criteria. An illustrative example is given.