Secondly, I Don't know if you recognize this or not, but if you interept the answer (even though its jogging on cluster), you ought to be able to get the NR residuals for the final three (I usually use 3 but you can enter all the more) iterations. The solution https://a-solution.info/