;
h = plot(domain, f(domain), '-r', 'LineWidth', 2.5);
axis([-4, 10, 0, 0.5]);
shade(f, 0, -1.5, 1.5, 'b', 0.005);
shade(f, 0,  4.5, 7.5, 'b', 0.005);
line([-4;10], [f(-1.5), f(-1.5)], 'Color', 'b', 'LineWidth', 2);
set(gca,'XTick', [], 'YTick', f(-1.5), 'YTickLabel', 'pMIN', 'fontsize', 20);
uistack(h, 'top');
printPmtkFigure HPD;
end

##### SOURCE END #####
--></body></html>