  We present a new approach to multi-frequency synthesis in radio astronomy.
Using Bayesian inference techniques, the new technique estimates the sky
brightness and the spectral index simultaneously. In principle, the bandwidth
of a wide-band observation can be fully exploited for sensitivity and
resolution, currently only limited by higher order effects like spectral
curvature. Employing this new approach, we further present a multi-frequency
extension to the imaging algorithm RESOLVE. In simulations, this new algorithm
outperforms current multi-frequency imaging techniques like MS-MF-CLEAN.
