We calculate the cross correlations of the order parameters, which take the form: We see that the correlation length ξ$_{p}$ near and above the magnetic phase transition is enhanced, whereas the cross correlations exhibit an oscillatory behavior and evolve contin- uously through T$_{RM}$ . In this paper we examined a simple Landau approach for a system with two single- component order parameters representing ferroelectricity and ferromagnetism, respectively. The obtained results are consistent with experiment. Magnetoelectric coupling introduces a strong renormalization of the ferromagnetic transition temperature (enhancing it by a factor of 50) causing the magnetic phase transition to be observed at T$_{RM}$ = 100 . 5 K . A number of coupling effects such as the enhancement of electric polarization, the anomaly in the dielectric susceptibility, and a fairly large negative magnetocapacitance, occur concomitantly in the vicinity of T$_{RM}$ . We introduced a simple extension of the Arrot plot which we called in text Arrot planes . We consider this concept to be potential useful in evaluation of the phase transition temperature while using only field dependence of unbounded with it order parameter. We have also estimated the contribution of the Gaussian fluctuations to the specific heat and have noted that short-range-order effects are not accounted properly. It would be important to understand the microscopic reasons of such a strong enhancement of the magnetic exchange interactions caused by a monoclinic lattice distortion which leads to the appearance of the ferroelectric dipole moments. 22 V. CONCLUSIONS g$_{pm}$ ( ⃗ r$_{1}$, ⃗ r$_{2}$ ) = 〈 δm ( ⃗ r$_{1}$ ) δp ( ⃗ r$_{2}$ ) 〉 (62) = ∑ k 〈 δm$_{1k}$δp$_{1k}$ + δm$_{2k}$δp$_{2k}$ 〉 e i ⃗ k · ⃗ r = 1 2 π 2 · k$_{B}$T C$_{m}$γ$_{me}$r ∫ k$_{max}$ 0 k sin( kr ) dk = k$_{B}$T 2 π $^{2}$C$_{m}$ γ$_{me}$r $_{3}$(sin( k$_{max}$r ) − rk$_{max}$ cos( k$_{max}$r )) .