More precise determination of work function based on Fermi–Dirac distribution and Fowler formula

More precise numerical method to simulate current–voltage of metal at fixed temperature is presented in this paper. The new algorithm for the simulation has been developed via Fermi–Dirac distribution step by step. These calculated characteristics are shown to remain in excellent agreement with the experimental ones, taken for a range of different metals, which strong supports the validity of the model. It is also shown that based on the Fowler formula, higher precise work function can be determined.