Secondly, I do not know if you realize this or not, but if you interept the answer (even though its functioning on cluster), you should be capable to obtain the NR residuals for the final three (I usually use three however , you can enter a lot more) iterations. The https://a-solution.info/