\section { Introduction}
Nous pr\'esentons ici, la m\'ethode des volumes finis d\'ecrite dans le 
chapitre pr\'ec\'edent appliqu\'ee aux milieux h\'et\'erog\`enes (les
coefficients  
du syst\`eme de Maxwell $\varepsilon$ et $\mu$ sont discontinus). Cela
demande une nouvelle 
formulation du flux num\'erique, qui tient compte des conditions de
saut \`a travers les interfaces de mat\'eriaux d'indices diff\'erents.
% 
Nous pr\'esentons \'egalement une 
am\'elioration significative de cette m\'ethode en abandonnant  
l'interpolation lin\'eaire dans le calcul  des gradients (comme 
pr\'esent\'ee dans \cite{DURL} et \cite{REMA1})
 au profit d'une d\'efinition de gradients 
'faibles' qui se calculent alors comme des flux.  Ce proc\'ed\'e a 
l'avantage de respecter la structure du maillage 
%qui peut \^etre alors
%hybride voire non conforme 
et de plus il est plus facile \`a 
impl\'ementer et moins co\^uteux que les proc\'ed\'es 
plus classiques.










