Added step-5: Iteratively solves Schroedinger and Poisson, using Fick's equation for the band edge.

Loop is implemented and there, but for some reason the eigensolver still fails after the first
update, even if only tiny updates are applied. Further investigations required.