% FIGURE 5.8.1
      function FUNA = FUNA581(X,K,I)
      global TKP1
%                              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(0.5*X^2*TKP1)*(PKx*PIx - X*TKP1*PK*PI);
