% Schematic46.m %
clear all; close all;
fs=18; lw=1;

Nx = 10; Ny = 10;
xl = linspace(0,1,Nx+1);
yl = linspace(0,1,Ny+1);
xp = linspace((xl(1)+xl(2))/2,(xl(end-1)+xl(end))/2,Nx);
yp = linspace((yl(1)+yl(2))/2,(yl(end-1)+yl(end))/2,Ny);

for i=1:Nx
    for j=1:Ny
        w(i,j) = (2*(i+j-1))/(Nx*Ny*(Nx+Ny));
    end
end

figure(1); clf; hold on; axis image; grid on
for i = 1:Nx+1
    line([xl(i) xl(i)],[yl(1) yl(end)],'color','k','linewidth',lw)
    end
for j = 1:Ny+1
    line([xl(1) xl(end)],[yl(j) yl(j)],'color','k','linewidth',lw)
end
for ii = 1:2
plot(xp(1:end-3*(ii-1)),yp(ii),'ko-','MarkerFace','k','linewidth',lw);
line([xp(1) xp(end-3*(ii-1))],[yp(ii) yp(ii)],'Color','k','linewidth',lw)
plot(xp(end-3*(ii-1))-.025,yp(ii),'k>','MarkerFace','k');
plot(xp(end-3*(ii-1):end),yp(ii),'ko');
end
for jj = ii:Ny
    plot(xp,yp(jj),'ko');
end

axis off;

print(figure(1),'Schematic46','-dpng')