A recursive algorithm is derived for subspace system identification which is flexible and applicable to continuous/discrete-time, time-invariant/varying stochastic systems. Efficacy of the algorithm is demonstrated by simulations, comparing with another proposed algorithm.