A master equation approach for interacting slow- and stationary-light polaritons