% FIGURE 5.4.1
      function FUNA = FUNA541(X,K,I)
%                              INTEGRAND FOR (K,I)-TH ELEMENT
%                              OF JACOBIAN MATRIX
      PK  = PHI(K,0,X);
      PKx = PHI(K,1,X);
      PI  = PHI(I,0,X);
      PIx = PHI(I,1,X);
      FUNA = PKx*PIx + 0.11/X^2*PK*PI;