In this paper two new recursive algorithms are presented for computing optimal control law of linear time-invariant singular systems with quadratic performance index by using the elegant properties of block-pulse functions (BPFs) and shifted Legendre polynomials (SLPs). Also a unified approach is given to solve the optimal control problem of singular systems via BPFs or SLPs. Two numerical examples are included to demonstrate the validity of the proposed algorithms and approach.