Macierz zawierająca współrzędne węzłów jest następująca
wn = 8; x0 = 0.0; xn = 1.0;
hn = 10; y0 = 0.0; yn = 3.0;
xx = linspace(x0, xn, wn)
yy = linspace(y0, yn, hn)
it = 1
LW = wn*hn
nodes = zeros(LW,2)
for i=1:wn
for j=1:hn
nodes(it,1) = xx(i)
nodes(it,2) = yy(j)
it = it + 1
end
end
Macierz zawierająca definicje elementów jest następująca
LE = (wn-1)*(hn-1)
elements = zeros(LE,4)
it = 1
for i=1:wn-1
for j=1:hn-1
elements(it,1) = (i-1)*hn+j
elements(it,2) = i*hn+j
elements(it,3) = i*hn+j+1
elements(it,4) = (i-1)*hn+j+1
it = it + 1
end
end
Pliki zawierające współrzędne węzłów oraz macierz elementów można pobrać z następujących linków:
• pliki zawierajace numery wezłów dla górnej i dolnej krawędzi:
,
//n- liczba wierszy
//m - liczba kolumn
//name - nazwa pliku np. 'dane.txt'
data=read(name,n,m)