A7_fem_50.out
coordinates =
     1     3
     1     9
     6    11
    11     9
    11     3
     6     1
     6     6
elements3 =
     1     2     7
     2     3     7
     3     4     7
     4     5     7
     5     6     7
     6     1     7
elements4 =
     []
dirichlet =
     1     2
     2     3
     3     4
     4     5
     5     6
     6     1
neumann =
     []
assembly from triangles A=
   (1,1)      -1.0667
   (2,1)       0.2667
   (6,1)       0.3000
   (7,1)       0.5000
   (1,2)       0.2667
   (2,2)      -1.0667
   (3,2)       0.3000
   (7,2)       0.5000
   (2,3)       0.3000
   (3,3)      -1.3600
   (4,3)       0.3000
   (7,3)       0.7600
   (3,4)       0.3000
   (4,4)      -1.0667
   (5,4)       0.2667
   (7,4)       0.5000
   (4,5)       0.2667
   (5,5)      -1.0667
   (6,5)       0.3000
   (7,5)       0.5000
   (1,6)       0.3000
   (5,6)       0.3000
   (6,6)      -1.3600
   (7,6)       0.7600
   (1,7)       0.5000
   (2,7)       0.5000
   (3,7)       0.7600
   (4,7)       0.5000
   (5,7)       0.5000
   (6,7)       0.7600
   (7,7)      -3.5200
assembly from triangles and rectangles A=
   (1,1)      -1.0667
   (2,1)       0.2667
   (6,1)       0.3000
   (7,1)       0.5000
   (1,2)       0.2667
   (2,2)      -1.0667
   (3,2)       0.3000
   (7,2)       0.5000
   (2,3)       0.3000
   (3,3)      -1.3600
   (4,3)       0.3000
   (7,3)       0.7600
   (3,4)       0.3000
   (4,4)      -1.0667
   (5,4)       0.2667
   (7,4)       0.5000
   (4,5)       0.2667
   (5,5)      -1.0667
   (6,5)       0.3000
   (7,5)       0.5000
   (1,6)       0.3000
   (5,6)       0.3000
   (6,6)      -1.3600
   (7,6)       0.7600
   (1,7)       0.5000
   (2,7)       0.5000
   (3,7)       0.7600
   (4,7)       0.5000
   (5,7)       0.5000
   (6,7)       0.7600
   (7,7)      -3.5200
forces from triangles b=
   (1,1)       8.3333
   (2,1)       8.3333
   (3,1)       7.5758
   (4,1)       8.3333
   (5,1)       8.3333
   (6,1)       7.5758
   (7,1)      24.2424
forces from triangles and rectangles b=
   (1,1)       8.3333
   (2,1)       8.3333
   (3,1)       7.5758
   (4,1)       8.3333
   (5,1)       8.3333
   (6,1)       7.5758
   (7,1)      24.2424
using g() add neumann conditions b=
   (1,1)       8.3333
   (2,1)       8.3333
   (3,1)       7.5758
   (4,1)       8.3333
   (5,1)       8.3333
   (6,1)       7.5758
   (7,1)      24.2424
BoundNodes=
     1
     2
     3
     4
     5
     6
using u_d() add dirichlet conditions u=
   (1,1)       3.6364
   (2,1)      23.2727
   (3,1)      40.5455
   (4,1)      45.0909
   (5,1)      25.4545
   (6,1)       7.8182
apply u to b, dirichlet conditions b=
   (1,1)       3.6606
   (2,1)      20.0242
   (3,1)      42.2085
   (4,1)      37.4788
   (5,1)      21.1152
   (6,1)       9.4812
   (7,1)     -61.2412
solve for A * u = b  u=
   (1,1)       3.6364
   (2,1)      23.2727
   (3,1)      40.5455
   (4,1)      45.0909
   (5,1)      25.4545
   (6,1)       7.8182
   (7,1)      17.3981
1  x= 1  y= 3  solution= 3.6364
2  x= 1  y= 9  solution= 23.2727
3  x= 6  y= 11  solution= 40.5455
4  x= 11  y= 9  solution= 45.0909
5  x= 11  y= 3  solution= 25.4545
6  x= 6  y= 1  solution= 7.8182
7  x= 6  y= 6  solution= 17.3636