forms the 7-diagonal matrix instead of the matrix consisting of 4 3-diagonal blocks and solves the system of equations with the help of the internal octave solver for banded matrices. Works in the case alpha(x)=-x, and delta(x)=0