% Schematic37.m %
clear; clf;
nt= 160; nx = nt; x = linspace(1,nx+1,nx+1);
h = x(2) - x(1); sh = 16; hold on 
for i = 1:sh:nx+1
    line([x(1) x(nx+1)],[x(i) x(i)],'color','k','linewidth',0.5)
    line([x(i) x(i)],[x(1) x(nx+1)],'color','k','linewidth',0.5)
end
for i = 9:sh:nx+1
    %left BC
    plot(x(1)-(0.0)*x(sh), x(i)+8,...
        'ko','MarkerEdgeColor','k','markersize',6)
    %right BC
    plot(x(nx+1)+0.0*x(sh), x(i)+8,'ko',...
        'MarkerEdgeColor','k','markerfacecolor','k','markersize',6);
    %bottom BC
    plot(x(i)+8,x(1)-0.0*x(sh),...
        'ko','MarkerEdgeColor','k','markersize',6);
    %top BC
    plot(x(i)+8,x(nx+1)+0.0*x(sh),'ko',...
        'MarkerEdgeColor','k','markerfacecolor','k','markersize',6);
end
plot(x(1)-0.0*x(sh), x(1)-0.0*x(sh),...
    'ko','MarkerEdgeColor','k','markerfacecolor','k','markersize',6);
plot(x(nx+1)+0.0*x(sh), x(nx+1)+0.0*x(sh),...
    'ko','MarkerEdgeColor','k','markerfacecolor','k','markersize',6);
plot(x(1)-0.0*x(sh), x(nx+1)+0.0*x(sh),...
    'ko','MarkerEdgeColor','k','markerfacecolor','k','markersize',6);
plot(x(nx+1)+0.0*x(sh), x(1)-0.0*x(sh),...
    'ko','MarkerEdgeColor','k','markerfacecolor','k','markersize',6);
axis image
line([x(1) x(nx+1) x(nx+1) x(1) x(1)],[x(1) x(1) x(nx+1) ...
x(nx+1) x(1)],'linew',1.8,'color','k');
axis([x(1)-x(sh) x(nx+1)+x(sh) x(1)-x(sh) x(nx+1)+x(sh)]); axis off

text('Interpreter','latex','String','$x_0$','Position',[-3 ...
-8.02],'FontSize',15)
text('Interpreter','latex','String','$x_1$','Position',[14 ...
-8.02],'FontSize',15)
text('Interpreter','latex','String','$x_2$','Position',[29 ...
-8.02],'FontSize',15)
text('Interpreter','latex','String','$\cdots$','Position',[80 ...
-8.02],'FontSize',15)
text('Interpreter','latex','String','$x_{N_x}$','Position',[155 ...
-8.02],'FontSize',15)
text('Interpreter','latex','String','$y_0$','Position',[-12 ...
1],'FontSize',15)
text('Interpreter','latex','String','$y_1$','Position',[-12 ...
16],'FontSize',15)
text('Interpreter','latex','String','$y_2$','Position',[-12 ...
32],'FontSize',15)
text('Interpreter','latex','String','$\vdots$','Position',[-9 ...
85],'FontSize',15)
text('Interpreter','latex','String','$y_{N_y}$','Position',[-16 ...
160],'FontSize',15)
set(gcf,'paperpositionmode','auto')
