The Yang-Baxter-like matrix equation AXA = XAX is reconsidered, and an infinite number of solutions that commute with any given complex square matrix A are found. Our results here are based on the fact that the matrix A can be replaced with its Jordan canonical form. We also discuss the explicit structure of the solutions obtained.

We obtain the coefficient matrices of the finite element (FE), finite volume (FV) and finite difference (FD) methods based on P1-conforming elements on a quasi-uniform mesh, in order to approximately solve a boundary value problem involving the elliptic Poisson equation. The three methods are shown to possess the same H1-stability and convergence. Some numerical tests are made, to compare the numerical results from the three methods and to review our theoretical results.