A16_fem_50.out coordinates = 2 3 3 3 3 2 2 2 1 4 2 4 3 4 4 4 4 3 4 2 4 1 3 1 2 1 1 1 1 2 1 3 elements3 = 1 2 4 2 3 4 5 6 16 6 1 16 6 7 1 7 2 1 7 8 2 8 9 2 2 9 3 9 10 3 3 10 12 10 11 12 4 3 13 3 12 13 15 4 14 4 13 14 16 1 15 1 4 15 elements4 = [] dirichlet = 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 5 neumann = [] assembly from triangles A= (1,1) -4.0000 (2,1) 1.0000 (4,1) 1.0000 (6,1) 1.0000 (16,1) 1.0000 (1,2) 1.0000 (2,2) -4.0000 (3,2) 1.0000 (7,2) 1.0000 (9,2) 1.0000 (2,3) 1.0000 (3,3) -4.0000 (4,3) 1.0000 (9,3) 0.0000 (10,3) 1.0000 (12,3) 1.0000 (13,3) 0.0000 (1,4) 1.0000 (3,4) 1.0000 (4,4) -4.0000 (13,4) 1.0000 (15,4) 1.0000 (5,5) -1.0000 (6,5) 0.5000 (16,5) 0.5000 (1,6) 1.0000 (5,6) 0.5000 (6,6) -2.0000 (7,6) 0.5000 (2,7) 1.0000 (6,7) 0.5000 (7,7) -2.0000 (8,7) 0.5000 (7,8) 0.5000 (8,8) -1.0000 (9,8) 0.5000 (2,9) 1.0000 (3,9) 0.0000 (8,9) 0.5000 (9,9) -2.0000 (10,9) 0.5000 (3,10) 1.0000 (9,10) 0.5000 (10,10) -2.0000 (11,10) 0.5000 (10,11) 0.5000 (11,11) -1.0000 (12,11) 0.5000 (3,12) 1.0000 (11,12) 0.5000 (12,12) -2.0000 (13,12) 0.5000 (3,13) 0.0000 (4,13) 1.0000 (12,13) 0.5000 (13,13) -2.0000 (14,13) 0.5000 (13,14) 0.5000 (14,14) -1.0000 (15,14) 0.5000 (4,15) 1.0000 (14,15) 0.5000 (15,15) -2.0000 (16,15) 0.5000 (1,16) 1.0000 (5,16) 0.5000 (15,16) 0.5000 (16,16) -2.0000 assembly from triangles and rectangles A= (1,1) -4.0000 (2,1) 1.0000 (4,1) 1.0000 (6,1) 1.0000 (16,1) 1.0000 (1,2) 1.0000 (2,2) -4.0000 (3,2) 1.0000 (7,2) 1.0000 (9,2) 1.0000 (2,3) 1.0000 (3,3) -4.0000 (4,3) 1.0000 (9,3) 0.0000 (10,3) 1.0000 (12,3) 1.0000 (13,3) 0.0000 (1,4) 1.0000 (3,4) 1.0000 (4,4) -4.0000 (13,4) 1.0000 (15,4) 1.0000 (5,5) -1.0000 (6,5) 0.5000 (16,5) 0.5000 (1,6) 1.0000 (5,6) 0.5000 (6,6) -2.0000 (7,6) 0.5000 (2,7) 1.0000 (6,7) 0.5000 (7,7) -2.0000 (8,7) 0.5000 (7,8) 0.5000 (8,8) -1.0000 (9,8) 0.5000 (2,9) 1.0000 (3,9) 0.0000 (8,9) 0.5000 (9,9) -2.0000 (10,9) 0.5000 (3,10) 1.0000 (9,10) 0.5000 (10,10) -2.0000 (11,10) 0.5000 (10,11) 0.5000 (11,11) -1.0000 (12,11) 0.5000 (3,12) 1.0000 (11,12) 0.5000 (12,12) -2.0000 (13,12) 0.5000 (3,13) 0.0000 (4,13) 1.0000 (12,13) 0.5000 (13,13) -2.0000 (14,13) 0.5000 (13,14) 0.5000 (14,14) -1.0000 (15,14) 0.5000 (4,15) 1.0000 (14,15) 0.5000 (15,15) -2.0000 (16,15) 0.5000 (1,16) 1.0000 (5,16) 0.5000 (15,16) 0.5000 (16,16) -2.0000 forces from triangles b= (1,1) 0.9091 (2,1) 0.9091 (3,1) 0.9091 (4,1) 0.9091 (5,1) 0.1515 (6,1) 0.4545 (7,1) 0.4545 (8,1) 0.3030 (9,1) 0.4545 (10,1) 0.4545 (11,1) 0.1515 (12,1) 0.4545 (13,1) 0.4545 (14,1) 0.3030 (15,1) 0.4545 (16,1) 0.4545 forces from triangles and rectangles b= (1,1) 0.9091 (2,1) 0.9091 (3,1) 0.9091 (4,1) 0.9091 (5,1) 0.1515 (6,1) 0.4545 (7,1) 0.4545 (8,1) 0.3030 (9,1) 0.4545 (10,1) 0.4545 (11,1) 0.1515 (12,1) 0.4545 (13,1) 0.4545 (14,1) 0.3030 (15,1) 0.4545 (16,1) 0.4545 using g() add neumann conditions b= (1,1) 0.9091 (2,1) 0.9091 (3,1) 0.9091 (4,1) 0.9091 (5,1) 0.1515 (6,1) 0.4545 (7,1) 0.4545 (8,1) 0.3030 (9,1) 0.4545 (10,1) 0.4545 (11,1) 0.1515 (12,1) 0.4545 (13,1) 0.4545 (14,1) 0.3030 (15,1) 0.4545 (16,1) 0.4545 BoundNodes= 5 6 7 8 9 10 11 12 13 14 15 16 using u_d() add dirichlet conditions u= (5,1) 5.5455 (6,1) 6.0909 (7,1) 7.0000 (8,1) 8.2727 (9,1) 6.3636 (10,1) 5.0000 (11,1) 4.1818 (12,1) 2.9091 (13,1) 2.0000 (14,1) 1.4545 (15,1) 2.2727 (16,1) 3.6364 apply u to b, dirichlet conditions b= (1,1) -8.8182 (2,1) -12.4545 (3,1) -7.0000 (4,1) -3.3636 (5,1) 0.8333 (6,1) 6.3636 (7,1) 7.2727 (8,1) 1.8939 (9,1) 6.5455 (10,1) 5.1818 (11,1) 0.3788 (12,1) 3.1818 (13,1) 2.2727 (14,1) -0.3788 (15,1) 2.4545 (16,1) 3.8182 solve for A * u = b u= (1,1) 4.1818 (2,1) 5.0909 (3,1) 3.7273 (4,1) 2.8182 (5,1) 5.5455 (6,1) 6.0909 (7,1) 7.0000 (8,1) 8.2727 (9,1) 6.3636 (10,1) 5.0000 (11,1) 4.1818 (12,1) 2.9091 (13,1) 2.0000 (14,1) 1.4545 (15,1) 2.2727 (16,1) 3.6364 1 x= 2 y= 3 solution= 4.1818 2 x= 3 y= 3 solution= 5.0909 3 x= 3 y= 2 solution= 3.7273 4 x= 2 y= 2 solution= 2.8182 5 x= 1 y= 4 solution= 5.5455 6 x= 2 y= 4 solution= 6.0909 7 x= 3 y= 4 solution= 7 8 x= 4 y= 4 solution= 8.2727 9 x= 4 y= 3 solution= 6.3636 10 x= 4 y= 2 solution= 5 11 x= 4 y= 1 solution= 4.1818 12 x= 3 y= 1 solution= 2.9091 13 x= 2 y= 1 solution= 2 14 x= 1 y= 1 solution= 1.4545 15 x= 1 y= 2 solution= 2.2727 16 x= 1 y= 3 solution= 3.6364