In our calculations we see a big discrepancy in specific heat near the phase transition T$_{M}$ . Here we present a brief discussion of possible way to improve accordance to the experimental data by considering a little bit more realistic model. We assume that our order parameters are three dimensional. Due to the fact that our main interest lies in region near ferromagnetic phase transition and because T$_{f}$ >> T$_{M}$ we can consider a following simplification: ⃗ P ( ⃗r ) = P ( ⃗r ) ⃗e$_{z}$ . The main idea of proposed improvement is to take into account the angular fluctuation of ⃗ M near the phase transition. If we apply the orientation of the order parameters as presented at Fig. 12 we obtain following Landau free energy potential In the spirit of mode-mode coupling approach we can present the following two terms from B1 as: and After substitution B2 and B3 in B1 the free energy potential will have the form: 26 [1] W. Eerenstein, N.D Mathur, J.F. Scott, Nature , 442 , 756–765(2006). Δ F = M 2 z $^{{}$a 2 $^{〈}$cos $^{2}$Θ $^{〉}$+ b 4 $^{〈}$M 2 z $^{〉 〈}$cos $^{2}$Θ $^{〉}$+ γ 2 P 2 $^{〈}$cos $^{2}$Θ 〉 } (B4) + cos $^{2}$Θ$^{{}$ a 2 $^{〈}$M 2 z $^{〉}$+ b 4 $^{〈}$M 2 z 〉 2 $^{〈}$cos $^{2}$Θ $^{〉}$+ γ 2 P 2 $^{〈}$M 2 z 〉 } + P 2 $^{{}$α 2 + β 4 $^{〈}$P 2 $^{〉}$+ γ 2 $^{〈}$M 2 z $^{〉 〈}$cos $^{2}$Θ 〉 } − $^{〈}$M 2 z $^{〉 〈}$cos $^{2}$Θ 〉 $^{{}$a 2 P 2 + b 4 $^{〈}$M 2 z $^{〉 〈}$cos $^{2}$Θ 〉 } ( M$_{z}$ cos(Θ)) 4 ≈ $^{〈}$M 2 z $^{〉 〈}$cos $^{2}$Θ $^{〉 {}$M 2 z + $^{〈}$M 2 z $^{〉}$cos $^{2}$Θ − $^{〈}$cos $^{2}$Θ $^{〉 〈}$M 2 z $^{〉}}$. (B3) ( M$_{z}$ cos(Θ)) 2 ≈ $^{〈}$M 2 z $^{〉}$cos $^{2}$Θ + $^{〈}$cos $^{2}$Θ $^{〉}$M 2 z − $^{〈}$M 2 z $^{〉 〈}$cos $^{2}$Θ 〉 (B2) Δ F = a 2 M 2 z cos $^{2}$Θ + b 4 M 4 z cos $^{4}$Θ + α 2 P 2 + β 4 P 4 + γ ( PM$_{z}$ ) $^{2}$cos$^{2}$Θ . (B1) FIG. 12: (Color online). Model picture of the or- der parameters arrange- ment.