EN
The method of computation of the latent heat based on the energy distribution histogram is applied to the standard 3D Ashkin-Teller (AT) model. Similarly as in the original method for the q-state Potts model for strong first order phase transitions, the characteristic histogram with two peaks in the critical region have been observed. Positions of two minima of negative logarithm of internal energy probability for samples of finite size show good linear scalability to the thermodynamic limit. The applicability of this method has been confirmed by proving that the latent heat values are consistent with the ones obtained by us using the analysis of the behavior of the cumulants of the type of Challa and of Lee-Kosterlitz. The presented method is far more efficient than the one based on those cumulants.