A probabilistic production costing method for multistates and multiblocks of units is presented in this paper. The mixtures of normal distributions are used to represent the power of the system and the load. The energy, reliability and production costs are calculated directly from the analytical formulae. Numeral examples demonstrate that the algorithm presented in the paper is accurate and efficient.