The TLM (Transmission Line Matrix) method, in time domain, is extended to account for the presence of dispersive and anisotropic media in electromagnetic structures or devices. The model is thoroughly constructed by using Maxwell's equations that make it a unified general TLM formulation. The theoretical derivation was revisited, starting with Maxwell's equations, without invoking circuit analogy. The procedure is general and can be applied to derive the algorithm for new extended TLM nodes. The case of anisotropic ferrite-based structures is studied.