Singular systems have been the subject of interest over the last two decades due to their many practical applications. On the other hand, there has been considerable interest in the application of various orthogonal functions such as Walsh, block pulse, Hermite and recently wavelets family in signal processing due to their great properties to approximate signals. This paper proposed a numerical method for analyzing time-invariant singular systems via Haar wavelet
