2
n = 10; // the density of discretisation
4
Point(40) = { 0, 0, 0, h};
5
Point(41) = { R, 0, 0, h};
6
Point(42) = { 0,-R, 0, h};
7
Point(43) = {-R, 0, 0, h};
8
Point(44) = { 0, R, 0, h};
9
Circle(21) = {41, 40, 42};
10
Circle(22) = {42, 40, 43};
11
Circle(23) = {43, 40, 44};
12
Circle(24) = {44, 40, 41};
13
Line Loop (51) = {21,22,23,24};
15
Plane Surface(61) = {51} ;
17
Mesh.RecombinationAlgorithm = 0;
18
Recombine Surface{61} = angle;
19
Physical Line("boundary") = {21,22,23,24};
20
Physical Surface("interior") = {61};