Recent experiments have shown that the superconducting energy gap in some cuprates is spatially inhomogeneous. Motivated by these experiments, and using exact diagonalization of a model d-wave Hamiltonian, combined with Monte Carlo simulations of a Ginzburg-Landau free energy functional, we have calculated the single-particle local density of states LDOS ͑ , r͒ of a model high-T c superconductor as a function of temperature. Our calculations include both quenched disorder in the pairing potential and thermal fluctuations in both phase and amplitude of the superconducting gap. Most of our calculations assume two types of superconducting regions: ␣ with a small gap and large superfluid density, and  with the opposite. If the  regions are randomly embedded in an ␣ host, the LDOS on the ␣ sites still has a sharp coherence peak at T = 0, but the  component does not, in agreement with experiment. An ordered arrangement of  regions leads to oscillations in the LDOS as a function of energy. The model leads to a superconducting transition temperature T c well below the pseudogap temperature T c0 and has a spatially varying gap at very low T, both consistent with experiments in underdoped Bi2212. Our calculated LDOS ͑ , r͒ shows coherence peaks for T Ͻ T c , which disappear for T Ͼ T c , in agreement with previous work considering phase but not amplitude fluctuations in a homogeneous superconductor. Well above T c , the gap in the LDOS disappears.