Auxiliary function approach for determining symmetry energy at suprasaturation densities