% backward facing step problem
% \bar{Omega} = [-2,0] x [0,1] \cup [0,4] x [-1,1] 

global NMVT NMEDG NMEL
global MCOORD MVERT MEDGE MADJ

NMEL=5;
NMVT=11;
NMEDG=15;

% Initialisierung

% Matrix mit Koordinaten der Ecken
MCOORD = zeros(NMVT,2);

% Matrix mit Eck-Nummern
MVERT = zeros(NMEL,4);

% Matrix mit Kanten-Nummern
MEDGE = zeros(NMEL,4);

% Matrix mit Kantennachbarschaften
MADJ = zeros(NMEL,4);


% Matrix-Belegung

MCOORD = [  0  -1;   % Ecke 1
	    2  -1;   % Ecke 2
            4  -1;   % Ecke 3 
	    4   0;   % Ecke 4
	    2   0;   % Ecke 5
            0   0;   % Ecke 6 
           -2   0;   % Ecke 7
           -2   1;   % Ecke 8
            0   1;   % Ecke 9
            2   1;   % Ecke 10
            4   1;   % Ecke 11
         ];


MVERT = [   1   2   5   6;   % Element 1
	    2   3   4   5;   % Element 2
            7   6   9   8;   % Element 3 
	    6   5  10   9;   % Element 4
	    5   4  11  10;   % Element 5
	];


MEDGE = [ 1   4   7   5;   % Element 1 
	  2   3   8   4;   % Element 2 
	  6  11  13  12;   % Element 3
	  7  10  14  11;   % Element 4
	  8   9  15  10;   % Element 5
        ];


MADJ = [ -1   2   4  -6;   % Element 1
	 -1  -2   5   1;   % Element 2 
	 -5   4  -3  -4;   % Element 3
	  1   5  -3   3;   % Element 4
	  2  -2  -3   4;   % Element 5
       ];






