  In this work we investigate the order of the phase transition of the spin-1
Baxter-Wu model. We used extensive entropic simulations to describe the
behavior of quantities which reveal the order of the phase transition. We
applyied finite-sizing scaling laws for continuous and discontinuous phase
transitions. Our results show that this system exhibits an indeterminacy
regarding the order of the phase transition, i.e., the results are conclusive
for both transitions, whether continuous or discontinuous. In such a scenario
we carried out a study of the configurations in the region of the phase
transition, which confirmed that the model seems to undergo a tetracritical
transition, with the coexistence of a ferromagnetic and three ferrimagnetic
configurations, suggesting that it may be a multicritical point belonging to a
critical line of an external or a crystalline fields, where the continuous and
the discontinuous phase transitions may coexist reflecting different features
of the system
