I am developing a framework for optimization in continuous domain and I need to compute the k most probable configurations (most probable explanations) in a Gaussian network. Any suggestion of sample source code or library (R (preferred), Matlab, C++, Python, Java, etc) is welcome.