A numerical modelling procedure in the time-domain of microstrip circuits containing non-linear devices such as varactor multipliers and halversTM is described. The two-dimensional Transmission Line Matrix (TLM) method has been extended to model elements with voltage-dependent capacitance (varactors). Waveforms and spectra computed with the TLM model compare well with measured results.