  A new set of the next-to-leading order (NLO) and the next-to-next-to-leading
order (NNLO) low-energy constants $L_i^r$ and $C_i^r$ in chiral perturbation
theory is obtained. These values are computed using the new experimental data
with a new calculation method. This method combines the traditional global fit
and Monte Carlo method together. The higher order contributions are estimated
with this method. The theoretical values of the observables provide good
convergence at each chiral dimension, except for the NNLO values of the $\pi K$
scattering lengths $a_0^{3/2}$ and $a_0^{1/2}$. The fitted values for $L_i^r$
at NLO are close to their results with the new method at NNLO; i.e., these
$L_i^r$ are nearly order-independent in this method. The estimated ranges for
$C_i^r$ are consistent with those in the literature, and their possible upper
or/and lower boundaries are given. The values of some linear combinations of
$C_i^r$ are also given, and they are more reliable. If one knows a more exact
value $C_i^r$, another $C_i^r$ can be obtained by these values.
