function [A_el,F_el] = cal_Melm_1D_a(iel,gradP,U);
//
   global CONNEC COOR NUMER ADRESS UC Nelt Nno Npt NI;
   global XG PG Ng;
//
   x  = COOR(CONNEC(iel,:),1)
   
   A_el = zeros(Nno,Nno); F_el = zeros(Nno,1);

   U_el= U(ADRESS(iel,:));

   jac = (x(2) - x(1))/2.;

   for ig=1:Ng    //  boucle pour l'integration numerique
      xg = XG(ig);
//
      [FP,DxiFP] = P2(xg);
      DxFP = DxiFP/jac;
      
      gp =  abs(U_el'*DxFP);
      mu = cal_mu(gp);

      for i=1:Nno
         F_el(i) = F_el(i) + PG(ig)*gradP*FP(i)*jac;
         for j=1:Nno
            A_el(i,j) = A_el(i,j) + mu * PG(ig) * DxFP(j) * DxFP(i) * jac;
         end;
      end;
   end;
endfunction

