function ShowOtherRegions()

close all
N = 30;
Region = 'SNLCDAHB';
for k=1:8
subplot(1,2,1)
fill([1 N N 1 1],[1 1 N N 1],'w')
A = numgrid(Region(k),N);
hold on
spy(A)
axis equal off
hold off
subplot(1,2,2)
B = delsq(A);
spy(B)
pause
end