pro p1b,x

; sheet 2 - problem 1b
  
x=[  79,  89,  77,  83,  76,  91,  87,  89,  84,  68,$
     77,  73,  71,  77,  76,  77,  64,  82,  65,  46,$
     21,  77,  66,  72,  74,  65,  70,  70,  73,  76,$
     37,  63,  63,  50,  68,  64,  73,  44,  46,  33,$
     69,  61,  65,  54,  52,  77,  60,  65,  61,  26,$
     64,  49,  62,  53,  55,  62,  43,  50,  51,  54,$
     24,  37,  30,  41,  34,  31,  49,  43,  39,  48,$
     89,  33,  23,  35,  35,  14,  21,  30,  33,  26]   


nx=n_elements(x)

mom=moment(x)

print,'raw data'
print,'mean   = ',mom(0)
print,'median = ',median(x,/even)
print
sigma=sqrt(mom(1))
sigma_N=sigma*sqrt((nx-1)/float(nx)) ; mom(1) = Var = sig^2 defined w.r.t. N-1
print,'sigma_N = ',sigma_N
print
print,'gam1_N = ',mom(2)*(sigma/sigma_N)^3  

print
b=sort(x)
print,x(b)

return
end