pro example_idl, n1, x, x1, ps=ps x=[23.5,23.2,23.8,23.7,23.4,23.5,24.7,24.4,23.1,23.9,23.8,23.3, $ 23.8,23.7,23.2,23.0,23.1,24.2,24.4,24.5,24.0,23.9,24.1,24.5,23.4] loadct,12 if keyword_set(ps) then begin set_plot,'ps' device,file='example.ps',/landscape,/color !p.charsize=1.5 endif mom=moment(x) print,mom mu=mom(0) sig=sqrt(mom(1)) x1=randomu(seed,n1) x1=3.*sig*(2.*x1-1.) x1=x1+mu !p.multi=[0,1,2] n=n_elements(x) number=indgen(n) plot,number, x, psym=1,/ynozero,title='sample 1',xtitle='n',ytitle='x(n)' plots,[0,n-1],[mu,mu],col=100 plots,[0,n-1],[mu-sig,mu-sig],lines=1,col=200 plots,[0,n-1],[mu+sig,mu+sig],lines=1,col=200 number1=indgen(n1) plot,number1,x1, psym=4,/ynozero,title='sample 2',xtitle='n1',ytitle='x1(n1)' mom1=moment(x1) mu1=mom(0) sig1=sqrt(mom(1)) plots,[0,n1-1],[mu1,mu1],col=100 plots,[0,n1-1],[mu1-sig1,mu1-sig1],lines=1,col=200 plots,[0,n1-1],[mu1+sig1,mu1+sig1],lines=1,col=200 !p.multi=0 if keyword_set(ps) then begin device,/close set_plot,'x' !p.charsize=1. endif return end