function y=saturmom(u);

% En entrée le moment magnétique u général, créé par la régulation ( capteurs+électronique+bobines )

% Ce moment est écrété proportionnellement en se calquant sur la composante la plus grande, pas sur le 

% vecteur complet. Donc seule, une des 3 composantes atteindra à chaque instant le seuil de saturation,
 
% fixé à 0.2 Am²

mt_mgn=u;

maxi=[abs(mt_mgn(1)),abs(mt_mgn(2)),abs(mt_mgn(3))];

test=max(maxi);

if test>=0.2

  y=0.2*mt_mgn/test;

end

if test<0.2

  y=mt_mgn;

% y=0.2*mt_mgn/test; Inutile de limiter, on peut garder la régulation naturelle

end

% Si le moment est en ligne, il sort en ligne, et s'il entre en colonne, sort en colonne

end