A Bayesian analysis of simultaneous equation models by combining recursive analytical and numerical approaches