Calculation of the susceptibility of an eight-site order-disorder model for ferroelectrics