
%%  SBP  operator approximating d/dx to order 4 in the interior and to
%%  order 2 at the boundary with corresponding diagonal norm H.
%% 
%%  The following lines in  matlab will generate the upper left
%%  corner and the lower right corner of the matrix hD_1 and the
%%  first and last elements of the diagonal of hH. Here h is the space step.


A11=[-24/17 59/34 -4/17 -3/34 0 0;
-1/2 0 1/2 0 0 0;
4/43 -59/86 0 59/86 -4/43 0;
3/98 0 -59/98 0 32/49 -4/49];

ANN=-A(4:-1:1,6:-1:1);

H1=[17/48 59/48 43/48 49/48];
HN=H1(4:-1:1];




