We discuss a new mechanism leading to a matrix product form for the stationary state of one-dimensional stochastic models. The corresponding algebra is quadratic and involves four different matrices. For the example of a coagulation-decoagulation model explicit four-dimensional representations are given and exact expressions for various physical quantities are recovered. We also find the general structure of $n$-point correlation functions at the phase transition.