Algorithmic product-form approximations of interacting stochastic models