Master equation approach for interacting slow- and stationary-light polaritons