The Monte Carlo simulations allowing a distinction between the 1st and the
2nd order phase transitions in the three-dimensional Ashkin-Teller spin-lat
tice model, a system with three order parameters, are described. The applie
d method allows a precise location of points on the phase diagram. Some cri
tical points on the phase boundaries have been calculated in the regions co
mplementary to the results recently described in literature.