Used to generate the Wold decomposition for this model

Series[(1/((1 - p B) (1 - B))), {B, 0, 4}]