function c = fairec(cchar,nbt);
%genere la matrice du tenseur d'elasticite

c0 = zeros(1,10);
for i=1:10
  c0(i) = eval(cchar(i,:));
end 
c = diag(c0)*ones(10,nbt);