demo_medianval Program

     program demo_medianval
     ! return median value
     use M_orderpack, only : medianval
     implicit none
     character(len=*),parameter :: g='(*(g0,1x))'
        write(*,g)'real   ',&
        medianval( [80.0,70.0,20.0,10.0,1000.0] )
        write(*,g)'integer',&
        medianval( [11, 22, 33, 44, 55, 66, 77, 88] )
        write(*,g)'double ',&
        medianval( [11.0d0, 22.0d0, 33.0d0, 66.0d0, 77.0d0, 88.0d0] )
     end program demo_medianval