demo_median.f90 Source File

     program demo_median
     ! calculate median value
     use M_orderpack, only : median
     implicit none
     character(len=*),parameter :: g='(*(g0,1x))'

        write(*,g) 'real   ',median(&
        [80.0,70.0,20.0,10.0,1000.0] )

        write(*,g) 'integer',median(&
        [11, 22, 33, 44, 55, 66, 77, 88] )

        write(*,g) 'double ',median(&
        [11.0d0,22.0d0,33.0d0,66.0d0,77.0d0,88.0d0])

     end program demo_median