1C Output from Public domain Ratfor, version 1.05
2      subroutine kuantiles(k,m,n,x)
3      integer i,j,k(m),m,n
4      double precision x(n)
5      j = 0
6      do23000 i = 1,m
7      call dsel05(k(i)-j,n-j,x(j+1))
8      j = k(i)
923000 continue
1023001 continue
11      return
12      end
13