We present two-dimensional Monte Carlo radiative transfer models for the disc of the eclipsing binary Epsilon Aurigae by fitting its spectral energy distribution (SED) from optical to far-infrared (IR) wavelengths. We also ...