The semiconductor industry is continuously striving to improve the performance of electron devices and circuits. It implies the need for better understanding of their basic behaviour. However, the miniaturization and complexity of devices led to a breakdown of the analytical models. As a consequence numerical analysis and simulation based on fundamental differential equations has become necessary.