stationary_p_sparse function

Sparse version of stationary_p