Next, I Do not know if you recognize this or not, but for those who interept the answer (even though its functioning on cluster), you have to be capable to get the NR residuals for the last three (I ordinarily use three however you can enter more) iterations. The solution https://a-solution.info/