Theoretical treatment of ferroelectric composites using Monte Carlo calculations