// donne  les points et les poids de Gauss  en fonction du nombre de points choisi pour 
// l'integration numerique
function [XG,PG] = hammer(Ng)
//
//
 XG = zeros(Ng,2); PG = zeros(Ng,1);
if Ng == 1  then
   XG(1,1) = 0.333333333333333; XG(1,2) =  0.33333333333333; PG(1) = .5;
end;
if Ng == 3  then
   XG(1,1) = 0.666666666666667; XG(1,2) = 0.166666666666667; PG(1) = 0.166666666666667;
   XG(2,1) = 0.166666666666667; XG(2,2) = 0.666666666666667; PG(2) = 0.166666666666667;         
   XG(3,1) = 0.166666666666667; XG(3,2) = 0.166666666666667; PG(3) = 0.166666666666667;
end;
if Ng == 4  then
   XG(1,1) = 0.333333333333333; XG(1,2) = 0.333333333333333; PG(1) = -0.281250;             
   XG(2,1) = 0.2              ; XG(2,2) = 0.2              ; PG(2) = 0.260416666666667;  
   XG(3,1) = 0.2              ; XG(3,2) = 0.6              ; PG(3) = 0.260416666666667;     
   XG(4,1) = 0.6              ; XG(4,2) = 0.2              ; PG(4) = 0.260416666666667;
end;

if Ng == 6  then
   XG(1,1) = 0.108103018168070; XG(1,2) = 0.445948490915965; PG(1) = 0.1116907948390055; 
   XG(2,1) = 0.445948490915965; XG(2,2) = 0.108103018168070; PG(2) = 0.1116907948390055;
   XG(3,1) = 0.445948490915965; XG(3,2) = 0.445948490915965; PG(3) = 0.1116907948390055;
   XG(4,1) = 0.816847572980459; XG(4,2) = 0.091576213509771; PG(4) = 0.0549758718276610;
   XG(5,1) = 0.091576213509771; XG(5,2) = 0.816847572980459; PG(5) = 0.0549758718276610;
   XG(6,1) = 0.091576213509771; XG(6,2) = 0.091576213509771; PG(6) = 0.0549758718276610;
end;


if Ng <> 1 & Ng <> 3 & Ng <>4 & Ng <> 6 then
  fprintf(%io(2)," Ng = %d n''est pas admis !!!!",Ng);
end;
endfunction
