% Gleichungssystem auf Seite 261, Matrix A
% als Bandmatrix mit IBWL=4 und IBWR=3
% (Hauptdiagonalelemente von A in Spalte 4)

clear all

A = [0  0  0  0  0  1 ;     % Matrix A ...
     0  0  0 -1  0  1 ;     % ist als ...
     0  1 -4  6 -4  1 ;     % Bandmatrix ...
     0  1 -4  6 -4  1 ;     % gespeichert, ...
     0  1 -4  6 -4  1 ;     % b muss kein ...
     0  1 -4  6 -4  1 ;     % Spaltenvektor ...
     0  1 -4  6 -4  1 ;     % sein
    -1  2  0 -2  1  0 ;
     1 -2  1  0  0  0] ;
 
b = [0 0 1/256 1/256 1/256 1/256 1/256 0 0] ;
   
v = gabamp (A , b , 4)  