title: Fitting a bimodal distribution to a set of values
body: <p>Given a 1D array of values, what is the simplest way to figure out what the best fit bimodal distribution to it is, where each 'mode' is a normal distribution? Or in other words, how can you find the combination of two normal distributions that bests reproduces the 1D array of values?</p>
  
  <p>Specifically, I'm interested in implementing this in python, but answers don't have to be language specific.</p>
  
  <p>Thanks!</p>