Monte Carlo simulation of a classical quadrupole solid