title: how do you generate random numbers from a bimodal Gaussian Probability Density Function in MATLAB?
body: <p>How do you generate numbers from a bimodal Gaussian PDF in MATLAB? </p>
  
  <p>For a unimodal symmetric Gaussian PDF it is really easy: </p>
  
  <pre><code>x=randn(1000,1); 
  </code></pre>
  
  <p>But now I want to draw 1000 numbers from a Gaussian PDF that has two peaks (or two means).  How do I do this? </p>