“…Therefore, the first solution, obtained with the function eig, works as an initial guess for the inverse iteration, which tracks the next physical solution, eliminat-ing the necessity of processing the complete spectrum [31]. Following that, a version of the inverse iteration algorithm used in this work is presented.…”