A statistical model of phase transition in squaric acid (C,H,O, and C,D,O,) crystals is proposed.The thermodynamical properties of the transition are investigated by a five-particle cluster approximation not taking into account the tunneling motion of protons. It is shown that the phase transition in the analysed model can be of both first and second order. A comparison with recent experimental data is carried out.