function y=by_nulle(u);

% Cette fonction est destinée à éviter une division par 0 lorsque By s'annule.
% Ce qui est imprévisible en tumbling

% U est By en entrée, composante du champ sur l'axe de tangage

if u>1e-15

  y=1/u;

else 

  y=1e15;

end   % *********** fin du test *****

end   % *********** fin programme *****