The reliable prediction of the circuit response prior to manufacturing demands realistic element models and accurate analysis methods. This analysis has added difficulties in the case of nonlinear circuits, such as those based on transistors and diodes, in which the superposition principle does not hold. This gives rise to a dependence of the circuit response on the input amplitude and the generation of harmonic and intermodulation frequencies, among other phenomena. The special characteristics of nonlinear circuits are often undesired, but they also make them the only choice to obtain essential functions such as frequency mixing, frequency multiplication, oscillation, and frequency division. This article presents a detailed summary of numerical simulation methods in both the time and frequency domains, highlighting their advantages, limitations, and providing application examples.