L=0.2;
W=0.1;
cl=0.01; // BUT OPTIONAL SO WE REMOVE IT FOR NOW
//points of the exterior boundary:
Point(1)={L,W,0,cl};
Point(2)={-L,W,0,cl};
Point(3)={-L,-W,0,cl};
Point(4)={L,-W,0,cl};
Line(11)={1,2};
Line(12)={2,3};
Line(13)={3,4};
Line(14)={4,1};
Line Loop(21)={11,12,13,14};
//points of the interior boundary:
Point(30)={0.025,0.0,0,cl};
Point(31)={0.0247947503456,0.00319692904211,0,cl};
Point(32)={0.024182371576,0.00634136459774,0,cl};
Point(33)={0.0231729189337,0.00938167512198,0,cl};
Point(34)={0.0217829676031,0.0122679388001,0,cl};
Point(35)={0.0200353405467,0.0149527632623,0,cl};
Point(36)={0.0179587337524,0.0173920637651,0,cl};
Point(37)={0.0155872450465,0.0195457870617,0,cl};
Point(38)={0.0129598142078,0.0213785690751,0,cl};
Point(39)={0.0101195835781,0.0228603155754,0,cl};
Point(40)={0.00711318966578,0.0239666963259,0,cl};
Point(41)={0.00398999737583,0.0246795445854,0,cl};
Point(42)={0.000801289439291,0.024987155405,0,cl};
Point(43)={-0.00240057564769,0.0248844778237,0,cl};
Point(44)={-0.00556302334891,0.0243731978045,0,cl};
Point(45)={-0.00863412636053,0.0234617105512,0,cl};
Point(46)={-0.011563457256,0.0221649826593,0,cl};
Point(47)={-0.0143029165031,0.0205043063649,0,cl};
Point(48)={-0.0168075222565,0.0185069499269,0,cl};
Point(49)={-0.0190361489592,0.0162057098827,0,cl};
Point(50)={-0.0209522026223,0.0136383725303,0,cl};
Point(51)={-0.0225242216976,0.0108470934779,0,cl};
Point(52)={-0.0237263936753,0.00787770545059,0,cl};
Point(53)={-0.0245389789248,0.00477896571753,0,cl};
Point(54)={-0.0249486348188,0.00160175549952,0,cl};
Point(55)={-0.0249486348188,-0.00160175549952,0,cl};
Point(56)={-0.0245389789248,-0.00477896571753,0,cl};
Point(57)={-0.0237263936753,-0.00787770545059,0,cl};
Point(58)={-0.0225242216976,-0.0108470934779,0,cl};
Point(59)={-0.0209522026223,-0.0136383725303,0,cl};
Point(60)={-0.0190361489592,-0.0162057098827,0,cl};
Point(61)={-0.0168075222565,-0.0185069499269,0,cl};
Point(62)={-0.0143029165031,-0.0205043063649,0,cl};
Point(63)={-0.011563457256,-0.0221649826593,0,cl};
Point(64)={-0.00863412636053,-0.0234617105512,0,cl};
Point(65)={-0.00556302334891,-0.0243731978045,0,cl};
Point(66)={-0.00240057564769,-0.0248844778237,0,cl};
Point(67)={0.000801289439291,-0.024987155405,0,cl};
Point(68)={0.00398999737583,-0.0246795445854,0,cl};
Point(69)={0.00711318966578,-0.0239666963259,0,cl};
Point(70)={0.0101195835781,-0.0228603155754,0,cl};
Point(71)={0.0129598142078,-0.0213785690751,0,cl};
Point(72)={0.0155872450465,-0.0195457870617,0,cl};
Point(73)={0.0179587337524,-0.0173920637651,0,cl};
Point(74)={0.0200353405467,-0.0149527632623,0,cl};
Point(75)={0.0217829676031,-0.0122679388001,0,cl};
Point(76)={0.0231729189337,-0.00938167512198,0,cl};
Point(77)={0.024182371576,-0.00634136459774,0,cl};
Point(78)={0.0247947503456,-0.00319692904211,0,cl};
Point(79)={0.025,-6.12323399574e-18,0,cl};
Line(80)={79,78};
Line(81)={78,77};
Line(82)={77,76};
Line(83)={76,75};
Line(84)={75,74};
Line(85)={74,73};
Line(86)={73,72};
Line(87)={72,71};
Line(88)={71,70};
Line(89)={70,69};
Line(90)={69,68};
Line(91)={68,67};
Line(92)={67,66};
Line(93)={66,65};
Line(94)={65,64};
Line(95)={64,63};
Line(96)={63,62};
Line(97)={62,61};
Line(98)={61,60};
Line(99)={60,59};
Line(100)={59,58};
Line(101)={58,57};
Line(102)={57,56};
Line(103)={56,55};
Line(104)={55,54};
Line(105)={54,53};
Line(106)={53,52};
Line(107)={52,51};
Line(108)={51,50};
Line(109)={50,49};
Line(110)={49,48};
Line(111)={48,47};
Line(112)={47,46};
Line(113)={46,45};
Line(114)={45,44};
Line(115)={44,43};
Line(116)={43,42};
Line(117)={42,41};
Line(118)={41,40};
Line(119)={40,39};
Line(120)={39,38};
Line(121)={38,37};
Line(122)={37,36};
Line(123)={36,35};
Line(124)={35,34};
Line(125)={34,33};
Line(126)={33,32};
Line(127)={32,31};
Line(128)={31,30};
Line(129)={30,79};
Line Loop(130)={80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129};
//Create a plane surface that can be meshed:
//Plane Surface(#)={boundary Line Loop,interior hole Line Loop}
Plane Surface(131)={21,130};
// If everything went well 25 to 29  are still available for labeling 
//Name the surfaces to allow assignment of boundary conditions
Physical Line(25)={14};// FREE exterior boundary where BC is FREE
// 26 will be exterior boundary where BC is set to Vo + interior set to 0 but both do not move by Newton method! so same physical line
Physical Line(26)={11,12,13,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129};// du=0 boundary 
//You also have to give a name to the Plane surface
Physical Surface(29)={131};
