The authors present a completely algebraic version of the Schwarz method, including an algebraic approach to finding the optimal operator or a sparse approximation thereof. The method is applied to any banded or block banded systems of linear equations. The authors prove that the optimized Schwarz method converges in two iterations. Numerical examples are also presented and discussed.