In this paper Haar functions are developed to approximate the solutions of continuous time linear system. Properties of Haar functions are first presented, and an explicit expression for the inverse of the Haar operational matrix is presented. Using the inverse of the Haar operational matrix the full order Stein equation should be solved in terms of the solutions of pure algebraic matrix equations, which reduces the computation time remarkably. Finally a numerical example is illustrated to demonstrate the validity of the proposed algorithm.