A=[ 2 3 0 0 0 0;
    -2 -4 1 2 0 0 ;
    4 7 1 1 1  0;
    0 -1 5 4 4 2;
    0 0 -2 -7 2 1;
    0 0 0 0 2 1];
b=[2; 1; 4; 3; -4; 1];
B=LU_band(A,2,2);
x=vorw(B,2,2,b);
x=rueck(B,2,2,x);