Abstract In this contribution we propose two numerical methods for the solution of the system of two non-linear coupled differential Ginzburg-Landau equations. These proposals are based firstly on taking a matrix view considering the quasi-linear coupled system, as a second option, considering the computational molecule with its respective restriction of the values and eigen-vectors of the matrix. We clearly and concisely obtained the eigenvalues that lead us to an optimal spatio-temporal convergence solution of said system of equations. We compare the numerical convergence times obtained by these three methods with the known time found by applying the link variable method.