% **************************************************************

%                     ROUTINE MATLAB 4.2

% **************************************************************

% PROGRAMME DE RECHERCHE DE LA MASSE UTILE MAXIMUM DANS UN CALCUL
% D'AMELIORATION DES PEERFORMANCES DU LANCEUR DIAMANT B

% **************************************************************
function y=massemax;
% masse=[];
massemaxi=0;
for x=900:5:2000;
  for y=0:5:x-788;
      if (x-787-y)>0,
         a=(x/(x-675))^1.12999*((x-787)/(x-787-y))^1.14006;
          b=1-0.1178156*exp(-0.0004811417*x)*a;
          z=(x-787-1.1365*y)*b;
          mu=x-787-1.1365*y-1.3*z;
          % masse=[masse mu];
          if mu>massemaxi
            massemaxi=mu;
            xmax=x;
            ymax=y;
          end

      end
   end
end

% RESULTAT

y=[massemaxi xmax ymax];
   
end;


% **************************************************************