We demonstrate a new multi-wavelength technique for two-dimensional parametric modelling of galaxy surface-brightness profiles, which we have incorporated into the widely used software GALFIT. Our new method, named GALFITM, extends GALFIT3's current single-band fitting process by simultaneously using multiple images of the same galaxy to constrain a wavelength-dependent model. Each standard profile parameter may vary as a function of wavelength, with a user-definable degree of smoothness, from constant to fully free. The performance of GALFITM is evaluated by fitting elliptical Sérsic profiles to ugriz imaging data for 4026 galaxies, comprising the original SDSS imaging for 163 low redshift (v