 The M-SHAKE algorithm solves the non-linear system of equations using Newton's method directly.In each iteration, the linear system of equations

is solved exactly using an LU decomposition.Each iteration costs formula_49 operations, yet the solution converges quadratically, requiring fewer iterations than SHAKE.