%http://tex.stackexchange.com/questions/99139/how-to-fill-an-area-interval-under-a-gaussian \documentclass[a4j,12pt]{jarticle} \usepackage{amsmath,amssymb} \usepackage[dvips]{graphics} \usepackage{pstricks,pst-func,pst-plot,pst-math,pstricks-add} \pagestyle{empty} \begin{document} \psset{yunit=4cm,xunit=3} \begin{pspicture}(-2,-0.2)(2,1.4) \psgrid[griddots=10,gridlabels=0pt, subgriddiv=0] \psaxes[Dy=0.25]{->}(0,0)(-2,0)(2,1.25) \uput[-90](6,0){x}\uput[0](0,1){y} \rput[lb](1,0.75){\textcolor{red}{$\sigma =0.5$}} \rput[lb](1,0.5){\textcolor{blue}{$\sigma =1$}} \rput[lb](-2,0.5){$f(x)=\dfrac{1}{\sigma\sqrt{2\pi}}\,e^{-\dfrac{(x-\mu)^2}{2\sigma{}^2}}$} \psGauss[linecolor=red, linewidth=2pt]{-2.95}{2.95}% \psGaussI[linewidth=1pt]{-2}{2}% \pscustom[fillstyle=solid,fillcolor=red!30]{% \psGauss[linewidth=1pt,nue=5]{0.5}{1}% \psline(1,0)(0.5,0)} \end{pspicture} \end{document}