The figure below compares the discrete solution with the exact analytical solution The error is largest at the right boundary where it is equal to solution. The error is largest at the right boundary where it is equal to approximately 15% : In a real CFD solution the number of unknowns is typically thousands to millions. A lot of work goes into solving these matrices efficiency matrices efficiency. Usually an iterative procedure is d used.