With r being the vector of the residuals, calculated as r = Mx -y, whose norm can be used to probe the error in the approximate solution. Since W1…8, K3, K4 and K5 are functions of our unknown solutions, M and y must be calculated using an initial guess of x, and the system solved for a new vector of solutions that can used for the next iteration until a satisfying result is reached.