We discuss and develop a compact MOSFET modeling scheme in order to obtain accurate descriptions of the drain current and its derivatives up to the 5th order. We have analyzed the physical effects which govern the behaviour of the 3rd derivative in long and deep-subrnicron channel MOSFETs. Our modeling agrees well with experimental data and describes continuous transitions between operating regimes, thanks to the use of continuous functions, which do not introduce any artificial peaks.