Various types of magnetic phases in two dimensional systems with competing interactions are studied by means of a simulation method. We consider only first (J1) and second (J2) antiferromagnetic neighbor interactions. The simulation is based on the Crystal Growth Probability Method using the cluster probabilities obtained within the Kramer–Wannier approximation. We focus our attention on the evolution of the magnetic ordering as a function of temperature and of the ratioV(≡J2/J1). We notice that the pattern of the excitations changes drastically close toV=0.5, and is very different from systems withVvalues distant from that multicritical point.