/h1><!--introduction--><!--/introduction--><pre class="codeinput"><span class="comment">% This file is from pmtk3.googlecode.com</span>

<span class="keyword">function</span> postDensityIntervals
</pre><pre class="codeinput">f = @(x)gaussProb(x, 0, 1) + gaussProb(x, 6, 1);
domain = -4:0.001:10;
h = plot(domain , f(domain), <span class="string">'-r'</span>, <span class="string">'LineWidth'</span>, 2.5);
axis([-4, 10, 0, 0.5]);
shade(f, 0, -1.5, 7.5, <span class="string">'b'</span>, 0.005);
set(gca, <span cla