//   calcule U dans le cas lineaire 
//   a partir de Fext  et la foncion cal_Melm
//
function U = cal_sol_lin(cal_Melm,Fext,U)
//
   global CONNEC COOR NUMER ADRESS UC Nelt Nno Npt NI;
   global XG PG Ng;
    
   [A,F]= assemblage(cal_Melm,Fext,U);

   M11 = A(1:NI,1:NI);
   FI =  F(1:NI);
   M12 = A(1:NI,NI+1:$);
   FI = FI - M12*UC;

   [h,rk] = lufact(M11);
   UI = lusolve(h,FI);
   ludel(h);
   U = [UI ; UC];

endfunction
