SUBROUTINE MEDER(MEDIAN, STERMD, XIL, XID, NM, B)
C++
C++   NM=  NUMBER OF INCOME INTERVALS
C++    B=  SUPPLIED BY STATISTICAL METHODS DIVISION; SAMPLE 
C++    ADJUSTMENT FACTOR (APPROX. EQUAL TO INVERSE OF SAMPLING
C++    PROBABILITY.
C++   IL=  ARRAY CONTAINING UPPER LIMITS OF INCOME INTERVALS
C++   WHERE IL(1) IS THE TOTAL        
C++   ID=  ARRAY CONTAINING DISTRIBUTIONS, WHERE ID(1)=TOTAL,
C++   WHERE THE NUMBERS ARE IN THOUSANDS.
C++
      IMPLICIT INTEGER (A-Z)
      LOGICAL DEBUG
      REAL  B, P, P2, P3
      DIMENSION XIL(NM),XID(NM)
      DEBUG=.FALSE.
C++
      DO 10 I=1,NM
      XID(I)=1000.*XID(I)
   10 CONTINUE
C++
C++       CALCULATING MEDIAN
C++
      MEDIAN=KPMEDF(XIL,XID,NM,50.)
C++
C++       CALCULATING STANDARD ERROR OF MEDIAN
C++
      P=SQRT(50.*50.*(B/1./XID(1)))
      P2=50.+P
      KP2=KPMEDF(XIL,XID,NM,P2)
      P3=50.-P
      KP3=KPMEDF(XIL,XID,NM,P3)
      STERMD=(KP3-KP2)/2.0+.5
C++
C++       DEBUG PRINTING
C++
      IF(DEBUG)PRINT *,' P,P2,P3,KP2,KP3,STER=',P,P2,P3,KP2,KP3,STERMD
C++
      RETURN
      END