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