“…Note that the resolution of the complex linear system can be performed by other procedures, such as Cholesky decomposition, conjugate gradient, tridiagonal matrix algorithm (TDMA), modified strongly implicit procedure (MSI), among others. We chose the Relaxation Gauss-Seidel method because of its mathematical simplicity and easy computational implementation ( Smith, 2004;Sperandio, Mendes, & Monken, 2006;Romeiro, Castro, Cirilo, & Natti, 2011;Ferreira et al, 2012;Pardo, Natti, Romeiro, & Cirilo, 2011;Ladeia, Romeiro, Natti, & Cirilo, 2013;Saita et al, 2017;Romeiro, Mangili, Costanzi, Cirilo, & Natti, 2017;Cirilo, Barba, Romeiro, & Natti, 2018;Cirilo, Petrovskii, Romeiro, & Natti, 2019;Agrawal, 2019). We show that our procedure converges.…”