A16_fem_50d.out enter debug level:0=none, 1=input and matrices, 2=assembly 2 debug = 2 A16_coord.dat = 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 A16_elem3.dat = 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 A16_elem4.dat = A16_dirich.dat = 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 5 A16_neum.dat = MM = stima3(VV) at j=1 nodes(1,2,4) VV= 2 3 3 3 2 2 MM= -1.0000 0.5000 0.5000 0.5000 -0.5000 0 0.5000 0 -0.5000 MM = stima3(VV) at j=2 nodes(2,3,4) VV= 3 3 3 2 2 2 MM= -0.5000 0.5000 0 0.5000 -1.0000 0.5000 0 0.5000 -0.5000 MM = stima3(VV) at j=3 nodes(5,6,16) VV= 1 4 2 4 1 3 MM= -1.0000 0.5000 0.5000 0.5000 -0.5000 0 0.5000 0 -0.5000 MM = stima3(VV) at j=4 nodes(6,1,16) VV= 2 4 2 3 1 3 MM= -0.5000 0.5000 0 0.5000 -1.0000 0.5000 0 0.5000 -0.5000 MM = stima3(VV) at j=5 nodes(6,7,1) VV= 2 4 3 4 2 3 MM= -1.0000 0.5000 0.5000 0.5000 -0.5000 0 0.5000 0 -0.5000 MM = stima3(VV) at j=6 nodes(7,2,1) VV= 3 4 3 3 2 3 MM= -0.5000 0.5000 0 0.5000 -1.0000 0.5000 0 0.5000 -0.5000 MM = stima3(VV) at j=7 nodes(7,8,2) VV= 3 4 4 4 3 3 MM= -1.0000 0.5000 0.5000 0.5000 -0.5000 0 0.5000 0 -0.5000 MM = stima3(VV) at j=8 nodes(8,9,2) VV= 4 4 4 3 3 3 MM= -0.5000 0.5000 0 0.5000 -1.0000 0.5000 0 0.5000 -0.5000 MM = stima3(VV) at j=9 nodes(2,9,3) VV= 3 3 4 3 3 2 MM= -1.0000 0.5000 0.5000 0.5000 -0.5000 0.0000 0.5000 0.0000 -0.5000 MM = stima3(VV) at j=10 nodes(9,10,3) VV= 4 3 4 2 3 2 MM= -0.5000 0.5000 0 0.5000 -1.0000 0.5000 0 0.5000 -0.5000 MM = stima3(VV) at j=11 nodes(3,10,12) VV= 3 2 4 2 3 1 MM= -1.0000 0.5000 0.5000 0.5000 -0.5000 0 0.5000 0 -0.5000 MM = stima3(VV) at j=12 nodes(10,11,12) VV= 4 2 4 1 3 1 MM= -0.5000 0.5000 0 0.5000 -1.0000 0.5000 0 0.5000 -0.5000 MM = stima3(VV) at j=13 nodes(4,3,13) VV= 2 2 3 2 2 1 MM= -1.0000 0.5000 0.5000 0.5000 -0.5000 0 0.5000 0 -0.5000 MM = stima3(VV) at j=14 nodes(3,12,13) VV= 3 2 3 1 2 1 MM= -0.5000 0.5000 0.0000 0.5000 -1.0000 0.5000 0.0000 0.5000 -0.5000 MM = stima3(VV) at j=15 nodes(15,4,14) VV= 1 2 2 2 1 1 MM= -1.0000 0.5000 0.5000 0.5000 -0.5000 0 0.5000 0 -0.5000 MM = stima3(VV) at j=16 nodes(4,13,14) VV= 2 2 2 1 1 1 MM= -0.5000 0.5000 0 0.5000 -1.0000 0.5000 0 0.5000 -0.5000 MM = stima3(VV) at j=17 nodes(16,1,15) VV= 1 3 2 3 1 2 MM= -1.0000 0.5000 0.5000 0.5000 -0.5000 0 0.5000 0 -0.5000 MM = stima3(VV) at j=18 nodes(1,4,15) VV= 2 3 2 2 1 2 MM= -0.5000 0.5000 0 0.5000 -1.0000 0.5000 0 0.5000 -0.5000 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 analytic solution= 4.1818 2 x= 3 y= 3 analytic solution= 5.0909 3 x= 3 y= 2 analytic solution= 3.7273 4 x= 2 y= 2 analytic solution= 2.8182 5 x= 1 y= 4 analytic solution= 5.5455 6 x= 2 y= 4 analytic solution= 6.0909 7 x= 3 y= 4 analytic solution= 7 8 x= 4 y= 4 analytic solution= 8.2727 9 x= 4 y= 3 analytic solution= 6.3636 10 x= 4 y= 2 analytic solution= 5 11 x= 4 y= 1 analytic solution= 4.1818 12 x= 3 y= 1 analytic solution= 2.9091 13 x= 2 y= 1 analytic solution= 2 14 x= 1 y= 1 analytic solution= 1.4545 15 x= 1 y= 2 analytic solution= 2.2727 16 x= 1 y= 3 analytic solution= 3.6364