8.2 square()

sqaure() で $[0,1]\times[0,1]$ を $10\times10$ 分割

mesh Th = square(10,10);

実は square() で長方形 $[a,b]\times[c,d]$ も分割できるようで、

mesh Th = square(m, n, [a+(b-a)*x,c+(d-c)*y])

とすれば良い。

// [1,3]×[2,5] を 20×30 分割
mesh Th = square(20,30, [1+2*x, 2+3*y]);
plot(Th);

図 4: [1,3]×[2,5] を 20×30 分割

なお、square() で作ったメッシュは、 下の辺、右の辺、上の辺、左の辺 (反時計回り) の順に 1, 2, 3, 4 というラベルがつけられている。