Most numerical MOS-simulators are based on Gummel's nonlinear relaxation method, for which severe convergence problems arise above threshold. A new hierarchical approach for the numerical simulation of MOS-transistors (MOST's) is presented, which overcomes these shortcomings resulting in an improvement of convergence speed of more than one order of magnitude.