SUBROUTINE tilasto(arvot, keskiarvo, keskihajonta)
  IMPLICIT NONE
  REAL, DIMENSION(:), INTENT(IN) :: arvot
  REAL, INTENT(OUT) :: keskiarvo, keskihajonta
  REAL :: neliosumma

  n  = SIZE(arvot)
  keskiarvo = SUM(arvot)/n
  neliosumma = SUM(arvot**2)
  keskihajonta = SQRT((neliosumma - summa**2/n)/(n-1))
END SUBROUTINE tilasto