SUBROUTINE DRNMVN(NR,K,RSIG,LDRSIG,R,LDR)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      DIMENSION RSIG(LDRSIG,*),R(LDR,*)
      call um0set(nr,k,r,ldr)
      do 20 ii=1,nr
      DO 10 J=1,K
      EPS=DRNNOF()
      DO 10 I=j,k
   10 r(ii,I)=r(ii,I)+rsig(j,i)*EPS
   20 continue
      RETURN
      END