INSTALL @lib$+"SORTLIB"
      Sort% = FN_sortinit(0,0)

      DIM a(6), b(5)
      a() = 4.1, 5.6, 7.2, 1.7, 9.3, 4.4, 3.2
      b() = 4.1, 7.2, 1.7, 9.3, 4.4, 3.2

      PRINT "Median of a() is " ; FNmedian(a())
      PRINT "Median of b() is " ; FNmedian(b())
      END

      DEF FNmedian(a())
      LOCAL C%
      C% = DIM(a(),1) + 1
      CALL Sort%, a(0)
      = (a(C% DIV 2) + a((C%-1) DIV 2)) / 2