clear

disp('Avez vous rentré dans le programme l''ordre des éléments ?')



u=['load c:\matthieu\matlab\coeurs.dat'];
eval(u);


wtox=coeurs;

[np,ne]=size(wtox);

nso=1;
nmg=2;
nal=3;
nsi=4;
nph=5;
%le soufre s'appelait vif-argent
nva=6;
npo=7;
nca=8;
nti=9;
ncr=10;
nmn=11;
nfe=12;
nni=13;

for i=1:np
if wtox(i,nsi)==0
   wtox(i,nsi)=0.01;
 disp('pas de silicium');  
disp(i)
end
%wtsi(i,:)=wtox(i,:)./wtox(i,nsi);
end


