Linear combination of two Gaussian functions: multimodal posterior

I think this might answer your question: Fitting a spectra of gaussians