  We calculate the density of states of an inhomogeneous superconductor in a
magnetic field where the positions of vortices are distributed completely at
random. We consider both the cases of s-wave and d-wave pairing. For both
pairing symmetries either the presence of disorder or increasing the density of
vortices enhances the low energy density of states. In the s-wave case the gap
is filled and the density of states is a power law at low energies. In the
d-wave case the density of states is finite at zero energy and it rises
linearly at very low energies in the Dirac isotropic case
(\alpha_D=t/\Delta_0=1, where t is the hopping integral and \Delta_0 is the
amplitude of the order parameter). For slightly higher energies the density of
states crosses over to a quadratic behavior. As the Dirac anisotropy increases
(as \Delta_0 decreases with respect to the hopping term) the linear region
decreases in width. Neglecting this small region the density of states
interpolates between quadratic and back to linear as \alpha_D increases. The
low energy states are strongly peaked near the vortex cores.
