Mesh yapısı oluşturma

%% MESH
M = 3; N = 3;
DELTAX = W/(M-1); DELTAY = H/(N-1);
Xp = zeros(M,N); Yp = zeros(M,N);
for i = 1:M-1
for j = 1:N
if i == 1
Xp(i,j) = 0;
else
Xp(i,j) = Xp(i-1,j) + DELTAX;
end

if j == 1
Yp(i,j) = 0.0;
else
Yp(i,j) = Yp(i,j-1) + DELTAY;
end
end
end
figure
plot(Xp, Yp, 'ko', 'MarkerFaceColor', 'r', 'LineWidth', 1.5)
hold on
for i = 2:M
plot([Xp(i,1)-DELTAX/2 Xp(i,1)-DELTAX/2], [0 H], '–k', 'LineWidth', 1);
hold on
end
for j = 2:N
plot([0 W], [Yp(1,j)-DELTAY/2 Yp(1,j)-DELTAY/2], '–k', 'LineWidth', 1);
hold on
end
xlabel('X-position \rightarrow'); ylabel('Y-position \rightarrow');
Merhababalar ısı transferi sorusu çözüyorum ve öncelikle mesh yapısı oluştmam gerekiyor ama ne yazikki sıkıntı yaşıyorum 
                                       .
                             .         .         .
                             .         .         .
bu şekilde olacak biçimde yani ev gibi bir model 7 noktadan oluşan model yapmam gerekiyor ama bir türlü yapamadım kodda nereleri düzeltmem gerekir veya nasıl yeni kod oluşturabilirim?