It was shown that the shape of gas-liquid chromatographic peak is expressed by the Gaussian distribution multiplied by a higher order algebraic equation of time. Application of the curve fitting methods, i.e., the Newton’s divided-difference formula and the orthogonal polynomials method, could well clear up the experimentally observed gas-liquid chromatograms of trailing and leading types, respectively. The asymmetric shapes of gas-liquid chromatographic peaks are attributed to the slowness of vaporization or dissolution in the partition process.