% FIGURE 5.11.1
      function FUNA = FUNA5111(X,K,I)
%                              INTEGRAND FOR (K,I)-TH ELEMENT
%                              OF "A" MATRIX
      PK  = PHI(K,0,X);
      PKx = PHI(K,1,X);
      PI  = PHI(I,0,X);
      PIx = PHI(I,1,X);
      FUNA = exp(-2*X)*(PKx*PIx - PK*PI);
