We have performed Monte Carlo simulations of the Ising model coupled to three-dimensional quantum gravity based on a summation over dynamical triangulations. These were done both in the microcanonical ensemble, with the number of points in the triangulation and the number of Ising spins fixed, and in the grand canoncal ensemble. We have investigated the two possible cases of the spins living on the vertices of the triangulation (``diect'' case) and the spins living in the middle of the tetrahedra (``dual'' case). We observed phase transitions which are probably second order, and found that the dual implementation more effectively couples the spins to the quantum gravity.

11 pages