For the solution of this system, we used the Newton-Raphson procedure, where we update our variables.
