Iteration updated to find rho using diracsolver

U,rho, U1 and rho1 saved for each iteration and errors saved. No plot
commands and specific imports from diracsolver should prevent crashing
on Andrey's box.