add a new fermi level calculation method in band.py and its example

[AsymmetryChou]

Considering the fermi level estimation method in band.py is not very accurate and might fail in the case of metal (like Au primitive cell), I add an another fermi level calculation method, which gets fermi level iteratively by sampling kpoints in Brillouinzone with or without time-reversal-symmetry.