find_MAP function

Find the Maximum A Posteriori Estimation