NOTE: This function only works with Cantera 1.7 or later
This second method is based on the fact that the equilibrium sound speed is realized
at the sonic point for a CJ detonation. The algorithm is different than the originally
developed method which uses a minimization proceedure.