We have obtained a deep 8 field XMM-Newton mosaic of M33 covering the galaxy out to the D_25_ isophote and beyond to a limiting 0.2-4.5keV unabsorbed flux of 5x10^-16^erg/cm2/s (L>4x10^34^erg/s at the distance of M33). These data allow complete coverage of the galaxy with high sensitivity to soft sources such as diffuse hot gas and supernova remnants (SNRs). Here, we describe the methods we used to identify and characterize 1296 point sources in the 8 fields. We compare our resulting source catalog to the literature, note variable sources, construct hardness ratios, classify soft sources, analyze the source density profile, and measure the X-ray luminosity function (XLF). As a result of the large effective area of XMM-Newton below 1keV, the survey contains many new soft X-ray sources. The radial source density profile and XLF for the sources suggest that only ~15% of the 391 bright sources with L>3.6x10^35^erg/s are likely to be associated with M33, and more than a third of these are known SNRs. The log(N)-log(S) distribution, when corrected for background contamination, is a relatively flat power law with a differential index of 1.5, which suggests that many of the other M33 sources may be high-mass X-ray binaries. Finally, we note the discovery of an interesting new transient X-ray source, which we are unable to classify.