function  [ AN, rhoL, gammaG, hLiquid0,  ...
            time, yRocket,  ...
            hLiquid, FTrust, FDrag ] ...
          = best_rocket( ParametersDataFileName )
