PROGRAM pisteparit
  ! Ohjelma simuloi pisteparin sijoittamista
  ! yksikköneliöön.
  IMPLICIT NONE
  REAL, DIMENSION(2) :: a, b
  REAL :: d, s = 0.0
  INTEGER :: i, n

  WRITE (*,*) 'Anna lukuparien lkm:'
  READ (*,*) n
  WRITE (*,*) 'Pistepareja:', n, ' kpl'

  IF (n > 0) THEN
    DO i = 1, n
      CALL RANDOM_NUMBER(a)
      CALL RANDOM_NUMBER(b)
      s = s + SQRT(SUM((a - b)**2))
    END DO
    d = s/n
    WRITE (*,*) 'Etäisyyden keskiarvo:', d
  ELSE
    WRITE (*,*) 'Lukumäärä negatiivinen!'
  END IF
END PROGRAM pisteparit