Thanks, ZapperZ!
Thanks, ZapperZ, but this page is not exactly what I need. For now it would be enough for me to calculate without nonlocal contribution and spin-orbit interaction. I am looking for something more simple. Anyway, you are right, I need this calculation for Monte-Carlo simulation...