DOUBLE PRECISION FUNCTION DBLINF(NRA,NCA,A,LDA,X,Y)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION A(LDA,*),X(*),Y(*),XTRAN(100),M(100)
      DO 10 J=1,NCA
         SUM=0D0
         DO 20 I=1,NRA
            SUM=SUM+X(I)*A(I,J)
 20         CONTINUE
         M(J)=SUM   
 10   CONTINUE
      SUM=0D0
      DO 30 I=1,NCA
         SUM=SUM+M(I)*Y(I)
 30      CONTINUE
      DBLINF=SUM
      RETURN
      END