A finite‐difference method for solving the discrete Boltzmann equations, which are the governing equations for a model gas in which molecules have many discrete velocities, is developed. The method is applied to three fundamental problems in rarefied gas flow to study the features of discrete‐velocity gases: normal shock wave structure, heat transfer between two parallel plates, and two‐dimensional vapor deposition. Two different discrete‐velocity gas models are used.