% Differenzenverfahren, Beispiel auf den Seiten 260 bis 262,
% Koeffizientenmatrix wird als "Sparse Matrix" gespeichert

clear all

ipos = [1  2 2 3  3 3  3 3 4  4 4  4 4 5  5 5  5 5 6  6 6  6 6 7  7 7  7 7  8 8  8 8 9  9 9] ;
jpos = [3  2 4 1  2 3  4 5 2  3 4  5 6 3  4 5  6 7 4  5 6  7 8 5  6 7  8 9  5 6  8 9 6  7 8] ; 
AA   = [1 -1 1 1 -4 6 -4 1 1 -4 6 -4 1 1 -4 6 -4 1 1 -4 6 -4 1 1 -4 6 -4 1 -1 2 -2 1 1 -2 1] ;

A = sparse (ipos,jpos,AA,9,9) ;

b = [0 ; 0 ; 1/256 ; 1/256 ; 1/256 ; 1/256 ; 1/256 ; 0 ; 0] ;
   
v = A\b  