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