.mu');
figure;
contour(grid1, grid2, p);
hold <span class="string">on</span>
ndx = 1:1:Nsamples;
plot(samples.mu(1,ndx), samples.mu(2,ndx), <span class="string">'.'</span>);
line([35 65], [35 65], <span class="string">'color'</span>, <span class="string">'r'</span>, <span class="string">'linewidth'</span>, 3);
<span class="comment">%set(gca, 'xlim', mu1range);</span>
<span class="comment">%set(gca, 'ylim', mu2range);</span>
xlabel(sprintf(<span class="string">'%s'</span>, <span class="string">'\mu_1'</span>