function y=storcham(u);

% En entrée les 3 composantes du champ magnétique réel u(1:3)
% ainsi que le test ( 1 ou 0 )de situation de t entre to et to+Ta u(4)

% En sortie la variable le champ magnétique dans cet intervalle

champ0(1:3)=u(1:3);

test=u(4);

if test==1
   
  champ=champ0;

  save champB.mat champ;

  y=champ0;

else

  load champB.mat;

  y=champ(1:3);

end

end