In plasma spraying, the size and initial velocity of indivdual feedstock particles injected into the flame are essentially random variables within the fairly wide ranges of the injection conditions and powder characteristics. As a result, the motion and heating of particles cannot be specified in a deterministic manner and so a simulation based on a single particle cannot provide a realistic prediction of the deposition process. In the present study, these random elements are simulated using a statistical method. The distributions of the motion and heating of a whole population of particles are simulated rather than those of a single particle.