lc=400;
Point(1) = {0, 0, 0, lc};
Point(2) = {0, 10000, 0, lc};
Point(3) = {10000, 10000, 0, lc};
Point(4) = {10000, 0, 0, lc};
Point(5) = {2000,2000,0,lc};
Point(6) = {4000,3000,0,lc};
Point(7) = {5000,7000,0,lc};
Point(8) = {8500,7000,0,lc};

Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 1};
Line(5) = {5, 6};

Line Loop(1) = {1, 2, 3, 4};
Plane Surface(1) = {1};
Point{7} In Surface{1};
Point{8} In Surface{1};
Line{5} In Surface{1};

Field[1] = Attractor;
Field[1].NodesList = {7,8};
Field[1].NNodesByEdge = 20;
Field[1].EdgesList = {5};

Field[2] = MathEval;
Field[2].F = "20*Fabs(Exp(-(F1/400)^2))";


Background Field = 2;



