In a previous article in this Journal (J. Chem. Educ. 2000, 77, 1495) we introduced a numerical method, namely, the finite-difference boundary-value method, for the solution of the one-dimensional Schrödinger equation and illustrated its application to the evaluation of energy levels and wave functions for hindered internal rotations. Here the method is used to determine, in combination with vibrational spectroscopy and statistical thermodynamics, the torsional potential in ethane. In particular two distinct approaches have been exploited: the first approach is based on the experimental frequency of torsional mode, and the second, less direct but historically more relevant, approach is based on the experimental heat capacity of ethane at various temperatures and on the frequencies of the other normal modes of vibration. The two approaches provide energy barriers in good agreement with each other, 12.35 and 11.74 kJ mol–1, respectively, and with the literature values. It is shown that the finite-difference boundary-value method, providing a great number of accurate energy levels, is ideally suited for the calculation of both energy transitions and the partition function for internal rotation. The latter is used to calculate the contribution of torsional mode to thermodynamic functions, such as heat capacity, entropy, and enthalpy. The results are in excellent agreement with those obtained from the tables of Pitzer (Lewis, G. N.; Randall, M.; Pitzer, K. S.; Brewer, L. Thermodynamics, 2nd ed.; McGraw-Hill: New York, 1961; Chapter 27).