Calculation of work function