wolffd@0: function h = plotgauss1d(mu, sigma2)
wolffd@0: % function h = plotgauss1d(mu, sigma^2)
wolffd@0: % Example
wolffd@0: % plotgauss1d(0,5); hold on; h=plotgauss1d(0,2);set(h,'color','r')
wolffd@0: 
wolffd@0: sigma = sqrt(sigma2);
wolffd@0: x = linspace(mu-3*sigma, mu+3*sigma, 100);
wolffd@0: p = gaussian_prob(x, mu, sigma2);
wolffd@0: h = plot(x, p, '-');