fem_check21_tri.c running Given 3 ux + 2 uy + u = 2 x + 3 y + 13 Analytic solution u(x,y) = 1 + 2 x + 3 y about to read triangles from 22_tri.tri triangles read from 22_tri.tri tri 0 has vertices 0 5 4 tri 1 has vertices 0 1 5 tri 2 has vertices 1 6 5 tri 3 has vertices 1 2 6 tri 4 has vertices 2 7 6 tri 5 has vertices 2 3 7 tri 6 has vertices 4 9 8 tri 7 has vertices 4 5 9 tri 8 has vertices 5 10 9 tri 9 has vertices 5 6 10 tri 10 has vertices 6 11 10 tri 11 has vertices 6 7 11 tri 12 has vertices 8 13 12 tri 13 has vertices 8 9 13 tri 14 has vertices 9 14 13 tri 15 has vertices 9 10 14 tri 16 has vertices 10 15 14 tri 17 has vertices 10 11 15 subtracting minvert=0 from all vertices, using base zero finding unique of nold=54 found unique of n=16 about to read boundary from 22_tri.bound Dirichlet boundaries from 22_tri.bound boundary segment 0 has vertices 1 2 boundary segment 1 has vertices 2 3 boundary segment 2 has vertices 3 4 boundary segment 3 has vertices 0 4 boundary segment 4 has vertices 4 8 boundary segment 5 has vertices 8 12 boundary segment 6 has vertices 12 13 boundary segment 7 has vertices 13 14 boundary segment 8 has vertices 14 15 boundary segment 9 has vertices 3 7 boundary segment 10 has vertices 7 11 boundary segment 11 has vertices 11 15 subtracting minvert=0 from all vertices, using base zero finding unique of nold=24 found unique of n=12 freevert na=16, nb=12 freevert nc free = 4 number of free vertices is 4 free vertex 5 free vertex 6 free vertex 9 free vertex 10 about to read coordinates from 22_tri.coord coordinates read from 22_tri.coord coordinate 0 at 0.00 , 0.00 coordinate 1 at 0.00 , 0.33 coordinate 2 at 0.00 , 0.66 coordinate 3 at 0.00 , 1.00 coordinate 4 at 0.33 , 0.00 coordinate 5 at 0.33 , 0.33 coordinate 6 at 0.33 , 0.66 coordinate 7 at 0.33 , 1.00 coordinate 8 at 0.66 , 0.00 coordinate 9 at 0.66 , 0.33 coordinate 10 at 0.66 , 0.66 coordinate 11 at 0.66 , 1.00 coordinate 12 at 1.00 , 0.00 coordinate 13 at 1.00 , 0.33 coordinate 14 at 1.00 , 0.66 coordinate 15 at 1.00 , 1.00 xmin=0, xmax=1, ymin=0, ymax=1 nvert=16, nbound=12, nfree=4, ntri=18 compute global stiffness matrix computing local stiffness matrix for triangle 0 nodes i1=0, i2=5, i3=4, area=0.054450 coord (0.000000,0.000000) (0.330000,0.330000) (0.330000,0.000000) cm matrix of am cm = ym solve for cm 1.00 -3.03 -0.00 cm matrix of am cm = ym solve for cm -0.00 0.00 3.03 cm matrix of am cm = ym solve for cm -0.00 3.03 -3.03 tri_int1 running with np=30, nv=496 integral=0.119007, at i=5, j=5 intg fg =0.253829, at i=5 integral=-0.160429, at i=0, j=5 integral=0.059571, at i=4, j=5 finished k=0 computing local stiffness matrix for triangle 1 nodes i1=0, i2=1, i3=5, area=0.054450 coord (0.000000,0.000000) (0.000000,0.330000) (0.330000,0.330000) cm matrix of am cm = ym solve for cm 1.00 0.00 -3.03 cm matrix of am cm = ym solve for cm -0.00 -3.03 3.03 cm matrix of am cm = ym solve for cm -0.00 3.03 0.00 tri_int1 running with np=30, nv=496 integral=0.174007, at i=5, j=5 intg fg =0.255337, at i=5 integral=-0.105429, at i=0, j=5 integral=-0.050429, at i=1, j=5 finished k=1 computing local stiffness matrix for triangle 2 nodes i1=1, i2=6, i3=5, area=0.054450 coord (0.000000,0.330000) (0.330000,0.660000) (0.330000,0.330000) cm matrix of am cm = ym solve for cm 1.00 -3.03 -0.00 cm matrix of am cm = ym solve for cm -1.00 0.00 3.03 cm matrix of am cm = ym solve for cm 1.00 3.03 -3.03 tri_int1 running with np=30, nv=496 integral=0.119007, at i=6, j=6 intg fg =0.271797, at i=6 integral=-0.160429, at i=1, j=6 integral=0.059571, at i=5, j=6 integral=0.064007, at i=5, j=5 intg fg =0.267406, at i=5 integral=-0.160429, at i=1, j=5 integral=0.114571, at i=6, j=5 finished k=2 computing local stiffness matrix for triangle 3 nodes i1=1, i2=2, i3=6, area=0.054450 coord (0.000000,0.330000) (0.000000,0.660000) (0.330000,0.660000) cm matrix of am cm = ym solve for cm 2.00 0.00 -3.03 cm matrix of am cm = ym solve for cm -1.00 -3.03 3.03 cm matrix of am cm = ym solve for cm -0.00 3.03 0.00 tri_int1 running with np=30, nv=496 integral=0.174007, at i=6, j=6 intg fg =0.273306, at i=6 integral=-0.105429, at i=1, j=6 integral=-0.050429, at i=2, j=6 finished k=3 computing local stiffness matrix for triangle 4 nodes i1=2, i2=7, i3=6, area=0.056100 coord (0.000000,0.660000) (0.330000,1.000000) (0.330000,0.660000) cm matrix of am cm = ym solve for cm 1.00 -3.03 -0.00 cm matrix of am cm = ym solve for cm -1.94 0.00 2.94 cm matrix of am cm = ym solve for cm 1.94 3.03 -2.94 tri_int1 running with np=30, nv=496 integral=0.069280, at i=6, j=6 intg fg =0.294164, at i=6 integral=-0.165290, at i=2, j=6 integral=0.114710, at i=7, j=6 finished k=4 computing local stiffness matrix for triangle 5 nodes i1=2, i2=3, i3=7, area=0.056100 coord (0.000000,0.660000) (0.000000,1.000000) (0.330000,1.000000) cm matrix of am cm = ym solve for cm 2.94 -0.00 -2.94 cm matrix of am cm = ym solve for cm -1.94 -3.03 2.94 cm matrix of am cm = ym solve for cm -0.00 3.03 0.00 tri_int1 running with np=30, nv=496 finished k=5 computing local stiffness matrix for triangle 6 nodes i1=4, i2=9, i3=8, area=0.054450 coord (0.330000,0.000000) (0.660000,0.330000) (0.660000,0.000000) cm matrix of am cm = ym solve for cm 2.00 -3.03 -0.00 cm matrix of am cm = ym solve for cm 0.00 0.00 3.03 cm matrix of am cm = ym solve for cm -1.00 3.03 -3.03 tri_int1 running with np=30, nv=496 integral=0.119007, at i=9, j=9 intg fg =0.265808, at i=9 integral=-0.160429, at i=4, j=9 integral=0.059571, at i=8, j=9 finished k=6 computing local stiffness matrix for triangle 7 nodes i1=4, i2=5, i3=9, area=0.054450 coord (0.330000,0.000000) (0.330000,0.330000) (0.660000,0.330000) cm matrix of am cm = ym solve for cm 1.00 0.00 -3.03 cm matrix of am cm = ym solve for cm 1.00 -3.03 3.03 cm matrix of am cm = ym solve for cm -1.00 3.03 0.00 tri_int1 running with np=30, nv=496 integral=-0.045993, at i=5, j=5 intg fg =0.264389, at i=5 integral=-0.105429, at i=4, j=5 integral=0.169571, at i=9, j=5 integral=0.174007, at i=9, j=9 intg fg =0.267316, at i=9 integral=-0.105429, at i=4, j=9 integral=-0.050429, at i=5, j=9 finished k=7 computing local stiffness matrix for triangle 8 nodes i1=5, i2=10, i3=9, area=0.054450 coord (0.330000,0.330000) (0.660000,0.660000) (0.660000,0.330000) cm matrix of am cm = ym solve for cm 2.00 -3.03 -0.00 cm matrix of am cm = ym solve for cm -1.00 0.00 3.03 cm matrix of am cm = ym solve for cm 0.00 3.03 -3.03 tri_int1 running with np=30, nv=496 integral=-0.155993, at i=5, j=5 intg fg =0.276458, at i=5 integral=0.114571, at i=10, j=5 integral=0.059571, at i=9, j=5 integral=0.119007, at i=10, j=10 intg fg =0.283776, at i=10 integral=-0.160429, at i=5, j=10 integral=0.059571, at i=9, j=10 integral=0.064007, at i=9, j=9 intg fg =0.279385, at i=9 integral=-0.160429, at i=5, j=9 integral=0.114571, at i=10, j=9 finished k=8 computing local stiffness matrix for triangle 9 nodes i1=5, i2=6, i3=10, area=0.054450 coord (0.330000,0.330000) (0.330000,0.660000) (0.660000,0.660000) cm matrix of am cm = ym solve for cm 2.00 0.00 -3.03 cm matrix of am cm = ym solve for cm -0.00 -3.03 3.03 cm matrix of am cm = ym solve for cm -1.00 3.03 0.00 tri_int1 running with np=30, nv=496 integral=-0.100993, at i=5, j=5 intg fg =0.277966, at i=5 integral=-0.050429, at i=6, j=5 integral=0.169571, at i=10, j=5 integral=-0.045993, at i=6, j=6 intg fg =0.282357, at i=6 integral=-0.105429, at i=5, j=6 integral=0.169571, at i=10, j=6 integral=0.174007, at i=10, j=10 intg fg =0.285285, at i=10 integral=-0.105429, at i=5, j=10 integral=-0.050429, at i=6, j=10 finished k=9 computing local stiffness matrix for triangle 10 nodes i1=6, i2=11, i3=10, area=0.056100 coord (0.330000,0.660000) (0.660000,1.000000) (0.660000,0.660000) cm matrix of am cm = ym solve for cm 2.00 -3.03 -0.00 cm matrix of am cm = ym solve for cm -1.94 0.00 2.94 cm matrix of am cm = ym solve for cm 0.94 3.03 -2.94 tri_int1 running with np=30, nv=496 integral=-0.160720, at i=6, j=6 intg fg =0.303489, at i=6 integral=0.114710, at i=11, j=6 integral=0.064710, at i=10, j=6 integral=0.069280, at i=10, j=10 intg fg =0.306506, at i=10 integral=-0.165290, at i=6, j=10 integral=0.114710, at i=11, j=10 finished k=10 computing local stiffness matrix for triangle 11 nodes i1=6, i2=7, i3=11, area=0.056100 coord (0.330000,0.660000) (0.330000,1.000000) (0.660000,1.000000) cm matrix of am cm = ym solve for cm 2.94 -0.00 -2.94 cm matrix of am cm = ym solve for cm -0.94 -3.03 2.94 cm matrix of am cm = ym solve for cm -1.00 3.03 0.00 tri_int1 running with np=30, nv=496 integral=-0.100720, at i=6, j=6 intg fg =0.305185, at i=6 integral=-0.055290, at i=7, j=6 integral=0.174710, at i=11, j=6 finished k=11 computing local stiffness matrix for triangle 12 nodes i1=8, i2=13, i3=12, area=0.056100 coord (0.660000,0.000000) (1.000000,0.330000) (1.000000,0.000000) cm matrix of am cm = ym solve for cm 2.94 -2.94 0.00 cm matrix of am cm = ym solve for cm -0.00 -0.00 3.03 cm matrix of am cm = ym solve for cm -1.94 2.94 -3.03 tri_int1 running with np=30, nv=496 finished k=12 computing local stiffness matrix for triangle 13 nodes i1=8, i2=9, i3=13, area=0.056100 coord (0.660000,0.000000) (0.660000,0.330000) (1.000000,0.330000) cm matrix of am cm = ym solve for cm 1.00 0.00 -3.03 cm matrix of am cm = ym solve for cm 1.94 -2.94 3.03 cm matrix of am cm = ym solve for cm -1.94 2.94 -0.00 tri_int1 running with np=30, nv=496 integral=-0.042387, at i=9, j=9 intg fg =0.284837, at i=9 integral=-0.108623, at i=8, j=9 integral=0.169710, at i=13, j=9 finished k=13 computing local stiffness matrix for triangle 14 nodes i1=9, i2=14, i3=13, area=0.056100 coord (0.660000,0.330000) (1.000000,0.660000) (1.000000,0.330000) cm matrix of am cm = ym solve for cm 2.94 -2.94 0.00 cm matrix of am cm = ym solve for cm -1.00 -0.00 3.03 cm matrix of am cm = ym solve for cm -0.94 2.94 -3.03 tri_int1 running with np=30, nv=496 integral=-0.155720, at i=9, j=9 intg fg =0.297366, at i=9 integral=0.118043, at i=14, j=9 integral=0.056377, at i=13, j=9 finished k=14 computing local stiffness matrix for triangle 15 nodes i1=9, i2=10, i3=14, area=0.056100 coord (0.660000,0.330000) (0.660000,0.660000) (1.000000,0.660000) cm matrix of am cm = ym solve for cm 2.00 0.00 -3.03 cm matrix of am cm = ym solve for cm 0.94 -2.94 3.03 cm matrix of am cm = ym solve for cm -1.94 2.94 -0.00 tri_int1 running with np=30, nv=496 integral=-0.104053, at i=9, j=9 intg fg =0.298826, at i=9 integral=-0.046957, at i=10, j=9 integral=0.169710, at i=14, j=9 integral=-0.042387, at i=10, j=10 intg fg =0.303350, at i=10 integral=-0.108623, at i=9, j=10 integral=0.169710, at i=14, j=10 finished k=15 computing local stiffness matrix for triangle 16 nodes i1=10, i2=15, i3=14, area=0.057800 coord (0.660000,0.660000) (1.000000,1.000000) (1.000000,0.660000) cm matrix of am cm = ym solve for cm 2.94 -2.94 0.00 cm matrix of am cm = ym solve for cm -1.94 0.00 2.94 cm matrix of am cm = ym solve for cm 0.00 2.94 -2.94 tri_int1 running with np=30, nv=496 integral=-0.160439, at i=10, j=10 intg fg =0.325596, at i=10 integral=0.118186, at i=15, j=10 integral=0.061519, at i=14, j=10 finished k=16 computing local stiffness matrix for triangle 17 nodes i1=10, i2=11, i3=15, area=0.057800 coord (0.660000,0.660000) (0.660000,1.000000) (1.000000,1.000000) cm matrix of am cm = ym solve for cm 2.94 0.00 -2.94 cm matrix of am cm = ym solve for cm 0.00 -2.94 2.94 cm matrix of am cm = ym solve for cm -1.94 2.94 -0.00 tri_int1 running with np=30, nv=496 integral=-0.103772, at i=10, j=10 intg fg =0.327246, at i=10 integral=-0.051814, at i=11, j=10 integral=0.174853, at i=15, j=10 finished k=17 tri 0 has vertices 0 5 4 tri 1 has vertices 0 1 5 tri 2 has vertices 1 6 5 tri 3 has vertices 1 2 6 tri 4 has vertices 2 7 6 tri 5 has vertices 2 3 7 tri 6 has vertices 4 9 8 tri 7 has vertices 4 5 9 tri 8 has vertices 5 10 9 tri 9 has vertices 5 6 10 tri 10 has vertices 6 11 10 tri 11 has vertices 6 7 11 tri 12 has vertices 8 13 12 tri 13 has vertices 8 9 13 tri 14 has vertices 9 14 13 tri 15 has vertices 9 10 14 tri 16 has vertices 10 15 14 tri 17 has vertices 10 11 15 vertices, coordinates and analytic values coordinate 0 at 0.00 , 0.00, uana=1.000000 coordinate 1 at 0.00 , 0.33, uana=1.990000 coordinate 2 at 0.00 , 0.66, uana=2.980000 coordinate 3 at 0.00 , 1.00, uana=4.000000 coordinate 4 at 0.33 , 0.00, uana=1.660000 coordinate 5 at 0.33 , 0.33, uana=2.650000 coordinate 6 at 0.33 , 0.66, uana=3.640000 coordinate 7 at 0.33 , 1.00, uana=4.660000 coordinate 8 at 0.66 , 0.00, uana=2.320000 coordinate 9 at 0.66 , 0.33, uana=3.310000 coordinate 10 at 0.66 , 0.66, uana=4.300000 coordinate 11 at 0.66 , 1.00, uana=5.320000 coordinate 12 at 1.00 , 0.00, uana=3.000000 coordinate 13 at 1.00 , 0.33, uana=3.990000 coordinate 14 at 1.00 , 0.66, uana=4.980000 coordinate 15 at 1.00 , 1.00, uana=6.000000 k computed stiffness matrix, if debug K(0,0)=1.000000 K(0,1)=0.000000 K(0,2)=0.000000 K(0,3)=0.000000 K(0,4)=0.000000 K(0,5)=0.000000 K(0,6)=0.000000 K(0,7)=0.000000 K(0,8)=0.000000 K(0,9)=0.000000 K(0,10)=0.000000 K(0,11)=0.000000 K(0,12)=0.000000 K(0,13)=0.000000 K(0,14)=0.000000 K(0,15)=0.000000 K(1,0)=0.000000 K(1,1)=1.000000 K(1,2)=0.000000 K(1,3)=0.000000 K(1,4)=0.000000 K(1,5)=0.000000 K(1,6)=0.000000 K(1,7)=0.000000 K(1,8)=0.000000 K(1,9)=0.000000 K(1,10)=0.000000 K(1,11)=0.000000 K(1,12)=0.000000 K(1,13)=0.000000 K(1,14)=0.000000 K(1,15)=0.000000 K(2,0)=0.000000 K(2,1)=0.000000 K(2,2)=1.000000 K(2,3)=0.000000 K(2,4)=0.000000 K(2,5)=0.000000 K(2,6)=0.000000 K(2,7)=0.000000 K(2,8)=0.000000 K(2,9)=0.000000 K(2,10)=0.000000 K(2,11)=0.000000 K(2,12)=0.000000 K(2,13)=0.000000 K(2,14)=0.000000 K(2,15)=0.000000 K(3,0)=0.000000 K(3,1)=0.000000 K(3,2)=0.000000 K(3,3)=1.000000 K(3,4)=0.000000 K(3,5)=0.000000 K(3,6)=0.000000 K(3,7)=0.000000 K(3,8)=0.000000 K(3,9)=0.000000 K(3,10)=0.000000 K(3,11)=0.000000 K(3,12)=0.000000 K(3,13)=0.000000 K(3,14)=0.000000 K(3,15)=0.000000 K(4,0)=0.000000 K(4,1)=0.000000 K(4,2)=0.000000 K(4,3)=0.000000 K(4,4)=1.000000 K(4,5)=0.000000 K(4,6)=0.000000 K(4,7)=0.000000 K(4,8)=0.000000 K(4,9)=0.000000 K(4,10)=0.000000 K(4,11)=0.000000 K(4,12)=0.000000 K(4,13)=0.000000 K(4,14)=0.000000 K(4,15)=0.000000 K(5,0)=-0.265857 K(5,1)=-0.210857 K(5,2)=0.000000 K(5,3)=0.000000 K(5,4)=-0.045857 K(5,5)=0.054042 K(5,6)=0.064143 K(5,7)=0.000000 K(5,8)=0.000000 K(5,9)=0.229143 K(5,10)=0.284143 K(5,11)=0.000000 K(5,12)=0.000000 K(5,13)=0.000000 K(5,14)=0.000000 K(5,15)=0.000000 K(6,0)=0.000000 K(6,1)=-0.265857 K(6,2)=-0.215719 K(6,3)=0.000000 K(6,4)=0.000000 K(6,5)=-0.045857 K(6,6)=0.054861 K(6,7)=0.059420 K(6,8)=0.000000 K(6,9)=0.000000 K(6,10)=0.234281 K(6,11)=0.289420 K(6,12)=0.000000 K(6,13)=0.000000 K(6,14)=0.000000 K(6,15)=0.000000 K(7,0)=0.000000 K(7,1)=0.000000 K(7,2)=0.000000 K(7,3)=0.000000 K(7,4)=0.000000 K(7,5)=0.000000 K(7,6)=0.000000 K(7,7)=1.000000 K(7,8)=0.000000 K(7,9)=0.000000 K(7,10)=0.000000 K(7,11)=0.000000 K(7,12)=0.000000 K(7,13)=0.000000 K(7,14)=0.000000 K(7,15)=0.000000 K(8,0)=0.000000 K(8,1)=0.000000 K(8,2)=0.000000 K(8,3)=0.000000 K(8,4)=0.000000 K(8,5)=0.000000 K(8,6)=0.000000 K(8,7)=0.000000 K(8,8)=1.000000 K(8,9)=0.000000 K(8,10)=0.000000 K(8,11)=0.000000 K(8,12)=0.000000 K(8,13)=0.000000 K(8,14)=0.000000 K(8,15)=0.000000 K(9,0)=0.000000 K(9,1)=0.000000 K(9,2)=0.000000 K(9,3)=0.000000 K(9,4)=-0.265857 K(9,5)=-0.210857 K(9,6)=0.000000 K(9,7)=0.000000 K(9,8)=-0.049052 K(9,9)=0.054861 K(9,10)=0.067615 K(9,11)=0.000000 K(9,12)=0.000000 K(9,13)=0.226087 K(9,14)=0.287753 K(9,15)=0.000000 K(10,0)=0.000000 K(10,1)=0.000000 K(10,2)=0.000000 K(10,3)=0.000000 K(10,4)=0.000000 K(10,5)=-0.265857 K(10,6)=-0.215719 K(10,7)=0.000000 K(10,8)=0.000000 K(10,9)=-0.049052 K(10,10)=0.055696 K(10,11)=0.062896 K(10,12)=0.000000 K(10,13)=0.000000 K(10,14)=0.231229 K(10,15)=0.293039 K(11,0)=0.000000 K(11,1)=0.000000 K(11,2)=0.000000 K(11,3)=0.000000 K(11,4)=0.000000 K(11,5)=0.000000 K(11,6)=0.000000 K(11,7)=0.000000 K(11,8)=0.000000 K(11,9)=0.000000 K(11,10)=0.000000 K(11,11)=1.000000 K(11,12)=0.000000 K(11,13)=0.000000 K(11,14)=0.000000 K(11,15)=0.000000 K(12,0)=0.000000 K(12,1)=0.000000 K(12,2)=0.000000 K(12,3)=0.000000 K(12,4)=0.000000 K(12,5)=0.000000 K(12,6)=0.000000 K(12,7)=0.000000 K(12,8)=0.000000 K(12,9)=0.000000 K(12,10)=0.000000 K(12,11)=0.000000 K(12,12)=1.000000 K(12,13)=0.000000 K(12,14)=0.000000 K(12,15)=0.000000 K(13,0)=0.000000 K(13,1)=0.000000 K(13,2)=0.000000 K(13,3)=0.000000 K(13,4)=0.000000 K(13,5)=0.000000 K(13,6)=0.000000 K(13,7)=0.000000 K(13,8)=0.000000 K(13,9)=0.000000 K(13,10)=0.000000 K(13,11)=0.000000 K(13,12)=0.000000 K(13,13)=1.000000 K(13,14)=0.000000 K(13,15)=0.000000 K(14,0)=0.000000 K(14,1)=0.000000 K(14,2)=0.000000 K(14,3)=0.000000 K(14,4)=0.000000 K(14,5)=0.000000 K(14,6)=0.000000 K(14,7)=0.000000 K(14,8)=0.000000 K(14,9)=0.000000 K(14,10)=0.000000 K(14,11)=0.000000 K(14,12)=0.000000 K(14,13)=0.000000 K(14,14)=1.000000 K(14,15)=0.000000 K(15,0)=0.000000 K(15,1)=0.000000 K(15,2)=0.000000 K(15,3)=0.000000 K(15,4)=0.000000 K(15,5)=0.000000 K(15,6)=0.000000 K(15,7)=0.000000 K(15,8)=0.000000 K(15,9)=0.000000 K(15,10)=0.000000 K(15,11)=0.000000 K(15,12)=0.000000 K(15,13)=0.000000 K(15,14)=0.000000 K(15,15)=1.000000 f computed forcing function and boundary, if debug F(0)=1.000000, ug[0]=0.000000 F(1)=1.990000, ug[1]=0.000000 F(2)=2.980000, ug[2]=0.000000 F(3)=4.000000, ug[3]=0.000000 F(4)=1.660000, ug[4]=0.000000 F(5)=1.595385, ug[5]=0.000000 F(6)=1.730299, ug[6]=0.000000 F(7)=4.660000, ug[7]=0.000000 F(8)=2.320000, ug[8]=0.000000 F(9)=1.693537, ug[9]=0.000000 F(10)=1.831759, ug[10]=0.000000 F(11)=5.320000, ug[11]=0.000000 F(12)=3.000000, ug[12]=0.000000 F(13)=3.990000, ug[13]=0.000000 F(14)=4.980000, ug[14]=0.000000 F(15)=6.000000, ug[15]=0.000000 ug computed Galerkin, Ua analytic, error ug[0]= 1.000, Ua= 1.000, err=0 ug[1]= 1.990, Ua= 1.990, err=0 ug[2]= 2.980, Ua= 2.980, err=0 ug[3]= 4.000, Ua= 4.000, err=0 ug[4]= 1.660, Ua= 1.660, err=0 ug[5]= 2.650, Ua= 2.650, err=-2.22045e-15 ug[6]= 3.640, Ua= 3.640, err=-2.22045e-15 ug[7]= 4.660, Ua= 4.660, err=0 ug[8]= 2.320, Ua= 2.320, err=0 ug[9]= 3.310, Ua= 3.310, err=2.66454e-15 ug[10]= 4.300, Ua= 4.300, err=7.10543e-15 ug[11]= 5.320, Ua= 5.320, err=0 ug[12]= 3.000, Ua= 3.000, err=0 ug[13]= 3.990, Ua= 3.990, err=0 ug[14]= 4.980, Ua= 4.980, err=0 ug[15]= 6.000, Ua= 6.000, err=0 maxerr=7.10543e-15, avgerr=8.88178e-16