sparse_abc.java running with user file abc.java xmin=-1.000, xmax=1.000, nx=7, hx=0.333 ymin=-1.000, ymax=1.000, ny=7, hy=0.333 internal cells zeroed build_coef for point ix=2, iy=2 deriv pdxx returns bb=12, aa=-1 16 -30 16 -1 deriv pdyy returns bb=12, aa=-1 16 -30 16 -1 deriv pdx returns bb=12, aa=1 -8 0 8 -1 deriv pdy returns bb=12, aa=1 -8 0 8 -1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 - 0.750 0.000 0.000 i=1 0.000 0.000 12.000 0.000 0.000 i=2 0.000 0.000 - 22.500 0.000 0.000 i=3 0.000 0.000 12.000 0.000 0.000 i=4 0.000 0.000 - 0.750 0.000 0.000 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0 0.062 - 0.500 0.000 0.500 - 0.062 i=1 - 0.500 4.000 - 0.000 - 4.000 0.500 i=2 0.000 - 0.000 0.000 0.000 - 0.000 i=3 0.500 - 4.000 0.000 4.000 - 0.500 i=4 - 0.062 0.500 - 0.000 - 0.500 0.062 pdyy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 - 0.750 12.000 - 22.500 12.000 - 0.750 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 pdx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.250 0.000 0.000 i=1 0.000 0.000 - 2.000 0.000 0.000 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.000 0.000 2.000 0.000 0.000 i=4 0.000 0.000 - 0.250 0.000 0.000 pdy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 0.250 - 2.000 0.000 2.000 - 0.250 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 working row 6, i=2, j=2, bi=0, bj=0 a1=0.30326532985631666, b1=1.3662650602409636, c1=4.873973929071194 d1=-0.22222222222222224, e1=-0.03703703703703705 f1=-1.6666666666666667, c=-41.22199257545743 ij=0, u[ii][ij]=5.4650602409638545 ij=5, u[ii][ij]=4.083628402720245 ij=10, u[ii][ij]=-5.4650602409638545 ij=15, u[ii][ij]=0.6831325301204818 ij=1, u[ii][ij]=58.561761222928396 ij=6, u[ii][ij]=-118.15454999253565 ij=11, u[ii][ij]=58.413613074780244 ij=16, u[ii][ij]=-3.646221187544136 ij=2, u[ii][ij]=-5.4650602409638545 ij=7, u[ii][ij]=3.1947395138313555 ij=12, u[ii][ij]=5.4650602409638545 ij=17, u[ii][ij]=-0.6831325301204818 ij=3, u[ii][ij]=0.6831325301204818 ij=8, u[ii][ij]=-0.17189344183668195 ij=13, u[ii][ij]=-0.6831325301204818 ij=18, u[ii][ij]=0.08539156626506023 cc=-41.22199257545743 ii=6, ut[ii][ii]=- 118.155, ut[ii][cs]=- 58.892 working row 11, i=2, j=3, bi=0, bj=0 a1=0.423240862445307, b1=1.4, c1=4.566937378852422 d1=0.11111111111111113, e1=-0.0 f1=-1.0, c=-14.084240661089666 ij=0, u[ii][ij]=-0.7 ij=5, u[ii][ij]=5.6 ij=10, u[ii][ij]=4.856668127121462 ij=15, u[ii][ij]=-5.6 ij=20, u[ii][ij]=0.7 ij=1, u[ii][ij]=-3.425203034139316 ij=6, u[ii][ij]=54.80324854622906 ij=11, u[ii][ij]=-113.27901042919889 ij=16, u[ii][ij]=54.80324854622906 ij=21, u[ii][ij]=-3.425203034139316 ij=2, u[ii][ij]=0.7 ij=7, u[ii][ij]=-5.6 ij=12, u[ii][ij]=5.301112571565906 ij=17, u[ii][ij]=5.6 ij=22, u[ii][ij]=-0.7 ij=3, u[ii][ij]=-0.0875 ij=8, u[ii][ij]=0.7 ij=13, u[ii][ij]=-0.34520842461175805 ij=18, u[ii][ij]=-0.7 ij=23, u[ii][ij]=0.0875 cc=-14.084240661089666 ii=11, ut[ii][ii]=- 113.279, ut[ii][cs]=- 15.545 working row 16, i=2, j=4, bi=0, bj=0 a1=0.5906802064328228, b1=1.3662650602409638, c1=4.20421655312113 d1=0.4444444444444444, e1=-0.03703703703703703 f1=-0.3333333333333335, c=14.01034218585618 cc=14.01034218585618 ii=16, ut[ii][ii]=- 108.219, ut[ii][cs]= 56.619 working row 7, i=3, j=2, bi=0, bj=0 a1=0.3582656552868946, b1=1.4, c1=4.307036550218772 d1=-0.33333333333333337, e1=0.0 f1=-0.6666666666666667, c=-19.998998499103198 ij=0, u[ii][ij]=-0.7 ij=5, u[ii][ij]=-0.35203257479850425 ij=10, u[ii][ij]=0.7 ij=15, u[ii][ij]=-0.0875 ij=1, u[ii][ij]=5.6 ij=6, u[ii][ij]=4.965854530109402 ij=11, u[ii][ij]=-5.6 ij=16, u[ii][ij]=0.7 ij=2, u[ii][ij]=51.684438602625264 ij=7, u[ii][ij]=-105.63596629054416 ij=12, u[ii][ij]=51.684438602625264 ij=17, u[ii][ij]=-3.230277412664079 ij=3, u[ii][ij]=-5.6 ij=8, u[ii][ij]=3.632521196776068 ij=13, u[ii][ij]=5.6 ij=18, u[ii][ij]=-0.7 ij=4, u[ii][ij]=0.7 ij=9, u[ii][ij]=-0.1853659081318376 ij=14, u[ii][ij]=-0.7 ij=19, u[ii][ij]=0.0875 cc=-19.998998499103198 ii=7, ut[ii][ii]=- 105.636, ut[ii][cs]=- 24.979 working row 12, i=3, j=3, bi=0, bj=0 a1=0.5, b1=1.4, c1=4.0 d1=0.0, e1=0.0 f1=0.0, c=7.0 ij=0, u[ii][ij]=0.0875 ij=5, u[ii][ij]=-0.7 ij=10, u[ii][ij]=-0.375 ij=15, u[ii][ij]=0.7 ij=20, u[ii][ij]=-0.0875 ij=1, u[ii][ij]=-0.7 ij=6, u[ii][ij]=5.6 ij=11, u[ii][ij]=6.0 ij=16, u[ii][ij]=-5.6 ij=21, u[ii][ij]=0.7 ij=2, u[ii][ij]=-3.0 ij=7, u[ii][ij]=48.0 ij=12, u[ii][ij]=-101.25 ij=17, u[ii][ij]=48.0 ij=22, u[ii][ij]=-3.0 ij=3, u[ii][ij]=0.7 ij=8, u[ii][ij]=-5.6 ij=13, u[ii][ij]=6.0 ij=18, u[ii][ij]=5.6 ij=23, u[ii][ij]=-0.7 ij=4, u[ii][ij]=-0.0875 ij=9, u[ii][ij]=0.7 ij=14, u[ii][ij]=-0.375 ij=19, u[ii][ij]=-0.7 ij=24, u[ii][ij]=0.0875 cc=7.0 ii=12, ut[ii][ii]=- 101.250, ut[ii][cs]= 7.000 working row 17, i=3, j=4, bi=0, bj=0 a1=0.6978062125430446, b1=1.4, c1=3.6372791742687083 d1=0.33333333333333326, e1=0.0 f1=0.6666666666666665, c=36.320037764136224 cc=36.320037764136224 ii=17, ut[ii][ii]=- 96.873, ut[ii][cs]= 87.945 working row 8, i=4, j=2, bi=0, bj=0 a1=0.42324086244530695, b1=1.3662650602409638, c1=3.5158117000465934 d1=-0.22222222222222232, e1=0.03703703703703704 f1=0.33333333333333304, c=2.940551924356532 cc=2.940551924356532 ii=8, ut[ii][ii]=- 88.295, ut[ii][cs]= 9.761 working row 13, i=4, j=3, bi=0, bj=0 a1=0.5906802064328229, b1=1.4, c1=3.2087751498278214 d1=0.11111111111111106, e1=0.0 f1=0.9999999999999998, c=30.27883488648057 cc=30.27883488648057 ii=13, ut[ii][ii]=- 84.488, ut[ii][cs]= 42.591 working row 18, i=4, j=4, bi=0, bj=0 a1=0.8243606353500639, b1=1.3662650602409638, c1=2.8460543240965297 d1=0.4444444444444443, e1=0.037037037037037014 f1=1.6666666666666663, c=62.017461601245515 ij=6, u[ii][ij]=0.08539156626506024 ij=11, u[ii][ij]=-0.6831325301204819 ij=16, u[ii][ij]=-0.5071593654014369 ij=21, u[ii][ij]=0.6831325301204819 ij=7, u[ii][ij]=-0.6831325301204819 ij=12, u[ii][ij]=5.465060240963855 ij=17, u[ii][ij]=9.003438735311878 ij=22, u[ii][ij]=-5.465060240963855 ij=8, u[ii][ij]=-2.125281483813138 ij=13, u[ii][ij]=34.07857781508428 ij=18, u[ii][ij]=-80.91766992088168 ij=23, u[ii][ij]=34.22672596323243 ij=9, u[ii][ij]=0.6831325301204819 ij=14, u[ii][ij]=-5.465060240963855 ij=19, u[ii][ij]=10.781216513089657 ij=24, u[ii][ij]=5.465060240963855 cc=62.017461601245515 ii=18, ut[ii][ii]=- 80.918, ut[ii][cs]= 140.198 standard central cells generated build_coef for point ix=1, iy=2 deriv pdxx returns bb=12, aa=11 -20 6 4 -1 deriv pdyy returns bb=12, aa=-1 16 -30 16 -1 deriv pdx returns bb=12, aa=-3 -10 18 -6 1 deriv pdy returns bb=12, aa=1 -8 0 8 -1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 8.250 0.000 0.000 i=1 0.000 0.000 - 15.000 0.000 0.000 i=2 0.000 0.000 4.500 0.000 0.000 i=3 0.000 0.000 3.000 0.000 0.000 i=4 0.000 0.000 - 0.750 0.000 0.000 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0 - 0.188 1.500 - 0.000 - 1.500 0.188 i=1 - 0.625 5.000 - 0.000 - 5.000 0.625 i=2 1.125 - 9.000 0.000 9.000 - 1.125 i=3 - 0.375 3.000 - 0.000 - 3.000 0.375 i=4 0.062 - 0.500 0.000 0.500 - 0.062 pdyy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 - 0.750 12.000 - 22.500 12.000 - 0.750 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 pdx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 - 0.750 0.000 0.000 i=1 0.000 0.000 - 2.500 0.000 0.000 i=2 0.000 0.000 4.500 0.000 0.000 i=3 0.000 0.000 - 1.500 0.000 0.000 i=4 0.000 0.000 0.250 0.000 0.000 pdy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.250 - 2.000 0.000 2.000 - 0.250 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 working row edge i=1, j=2 working row 5, i=1, j=2, bi=-1, bj=0 a1=0.25670855951629595, b1=1.2741573033707865, c1=5.2802023121535875 d1=0.11111111111111116, e1=-0.0740740740740741 f1=-2.666666666666667, c=-57.114459389317965 ij=0, u[ii][ij]=69.88136241084513 ij=5, u[ii][ij]=-125.5996248606446 ij=10, u[ii][ij]=56.84349308084097 ij=15, u[ii][ij]=-3.1452849009899304 ij=1, u[ii][ij]=-11.467415730337079 ij=6, u[ii][ij]=1.655188517823332 ij=11, u[ii][ij]=11.467415730337079 ij=16, u[ii][ij]=-1.4334269662921348 ij=2, u[ii][ij]=3.8224719101123594 ij=7, u[ii][ij]=0.6034590118822212 ij=12, u[ii][ij]=-3.8224719101123594 ij=17, u[ii][ij]=0.47780898876404493 ij=3, u[ii][ij]=-0.6370786516853932 ij=8, u[ii][ij]=-0.1647536418594442 ij=13, u[ii][ij]=0.6370786516853932 ij=18, u[ii][ij]=-0.07963483146067415 cc=-57.114459389317965 ii=5, ut[ii][ii]=- 125.600, ut[ii][cs]=- 71.367 working row edge i=1, j=3 working row 10, i=1, j=3, bi=-1, bj=0 a1=0.3582656552868946, b1=1.4, c1=4.9731657619348155 d1=0.44444444444444453, e1=-0.0 f1=-2.0, c=-30.704761499614747 ij=0, u[ii][ij]=-4.604874321451112 ij=5, u[ii][ij]=66.67798914321779 ij=10, u[ii][ij]=-120.38132558394788 ij=15, u[ii][ij]=52.677989143217786 ij=20, u[ii][ij]=-2.8548743214511116 ij=1, u[ii][ij]=1.575 ij=6, u[ii][ij]=-12.6 ij=11, u[ii][ij]=3.612195448791026 ij=16, u[ii][ij]=12.6 ij=21, u[ii][ij]=-1.575 ij=2, u[ii][ij]=-0.5249999999999999 ij=7, u[ii][ij]=4.199999999999999 ij=12, u[ii][ij]=0.408130299194017 ij=17, u[ii][ij]=-4.199999999999999 ij=22, u[ii][ij]=0.5249999999999999 ij=3, u[ii][ij]=0.0875 ij=8, u[ii][ij]=-0.7 ij=13, u[ii][ij]=-0.1575881303540598 ij=18, u[ii][ij]=0.7 ij=23, u[ii][ij]=-0.0875 cc=-30.704761499614747 ii=10, ut[ii][ii]=- 120.381, ut[ii][cs]=- 28.077 working row edge i=1, j=4 working row 15, i=1, j=4, bi=-1, bj=0 a1=0.49999999999999983, b1=1.2741573033707865, c1=4.610444936203524 d1=0.7777777777777778, e1=-0.07407407407407406 f1=-1.3333333333333335, c=-3.318971201755444 cc=-3.318971201755444 ii=15, ut[ii][ii]=- 114.513, ut[ii][cs]= 23.941 build_coef for point ix=3, iy=2 deriv pdxx returns bb=12, aa=-1 4 6 -20 11 deriv pdyy returns bb=12, aa=-1 16 -30 16 -1 deriv pdx returns bb=12, aa=-1 6 -18 10 3 deriv pdy returns bb=12, aa=1 -8 0 8 -1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 - 0.750 0.000 0.000 i=1 0.000 0.000 3.000 0.000 0.000 i=2 0.000 0.000 4.500 0.000 0.000 i=3 0.000 0.000 - 15.000 0.000 0.000 i=4 0.000 0.000 8.250 0.000 0.000 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0 - 0.062 0.500 - 0.000 - 0.500 0.062 i=1 0.375 - 3.000 0.000 3.000 - 0.375 i=2 - 1.125 9.000 - 0.000 - 9.000 1.125 i=3 0.625 - 5.000 0.000 5.000 - 0.625 i=4 0.188 - 1.500 0.000 1.500 - 0.188 pdyy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 0.000 0.000 0.000 0.000 0.000 i=3 - 0.750 12.000 - 22.500 12.000 - 0.750 i=4 0.000 0.000 0.000 0.000 0.000 pdx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 - 0.250 0.000 0.000 i=1 0.000 0.000 1.500 0.000 0.000 i=2 0.000 0.000 - 4.500 0.000 0.000 i=3 0.000 0.000 2.500 0.000 0.000 i=4 0.000 0.000 0.750 0.000 0.000 pdy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.250 - 2.000 0.000 2.000 - 0.250 i=4 0.000 0.000 0.000 0.000 0.000 working row edge i=5, j=2 working row 9, i=5, j=2, bi=1, bj=0 a1=0.49999999999999983, b1=1.2741573033707865, c1=2.4115684681094205 d1=0.11111111111111088, e1=0.07407407407407408 f1=1.3333333333333328, c=25.104017586646457 cc=25.104017586646457 ii=9, ut[ii][ii]=- 60.149, ut[ii][cs]=- 28.906 working row edge i=5, j=3 working row 14, i=5, j=3, bi=1, bj=0 a1=0.6978062125430446, b1=1.4, c1=2.1045319178906485 d1=0.44444444444444425, e1=0.0 f1=1.9999999999999996, c=54.46724693077414 cc=54.46724693077414 ii=14, ut[ii][ii]=- 54.708, ut[ii][cs]=- 52.636 working row edge i=5, j=4 working row 19, i=5, j=4, bi=1, bj=0 a1=0.9738670205273376, b1=1.2741573033707865, c1=1.7418110921593568 d1=0.7777777777777775, e1=0.07407407407407403 f1=2.666666666666666, c=90.54140587655328 ij=6, u[ii][ij]=-0.07963483146067415 ij=11, u[ii][ij]=0.6370786516853932 ij=16, u[ii][ij]=-0.9248447098399477 ij=21, u[ii][ij]=-0.6370786516853932 ij=7, u[ii][ij]=0.47780898876404493 ij=12, u[ii][ij]=-3.8224719101123594 ij=17, u[ii][ij]=4.088267728248679 ij=22, u[ii][ij]=3.8224719101123594 ij=8, u[ii][ij]=-1.4334269662921348 ij=13, u[ii][ij]=11.467415730337079 ij=18, u[ii][ij]=0.8824015923730211 ij=23, u[ii][ij]=-11.467415730337079 ij=9, u[ii][ij]=-0.4914914859942575 ij=14, u[ii][ij]=14.3827984409102 ij=19, u[ii][ij]=-49.18764377038448 ij=24, u[ii][ij]=27.420667770914363 cc=90.54140587655328 ii=19, ut[ii][ii]=- 49.188, ut[ii][cs]=- 67.994 build_coef for point ix=2, iy=1 deriv pdxx returns bb=12, aa=-1 16 -30 16 -1 deriv pdyy returns bb=12, aa=11 -20 6 4 -1 deriv pdx returns bb=12, aa=1 -8 0 8 -1 deriv pdy returns bb=12, aa=-3 -10 18 -6 1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 - 0.750 0.000 0.000 0.000 i=1 0.000 12.000 0.000 0.000 0.000 i=2 0.000 - 22.500 0.000 0.000 0.000 i=3 0.000 12.000 0.000 0.000 0.000 i=4 0.000 - 0.750 0.000 0.000 0.000 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0 - 0.188 - 0.625 1.125 - 0.375 0.062 i=1 1.500 5.000 - 9.000 3.000 - 0.500 i=2 - 0.000 - 0.000 0.000 - 0.000 0.000 i=3 - 1.500 - 5.000 9.000 - 3.000 0.500 i=4 0.188 0.625 - 1.125 0.375 - 0.062 pdyy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 8.250 - 15.000 4.500 3.000 - 0.750 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 pdx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.250 0.000 0.000 0.000 i=1 0.000 - 2.000 0.000 0.000 0.000 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.000 2.000 0.000 0.000 0.000 i=4 0.000 - 0.250 0.000 0.000 0.000 pdy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 - 0.750 - 2.500 4.500 - 1.500 0.250 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 working row edge i=2, j=1 working row 1, i=2, j=1, bi=0, bj=-1 a1=0.2172991042535391, b1=1.2741573033707865, c1=5.133874757704843 d1=-0.5555555555555556, e1=-0.1481481481481482 f1=-2.3333333333333335, c=-66.26449653079774 ij=0, u[ii][ij]=10.089486879007513 ij=5, u[ii][ij]=-11.467415730337079 ij=10, u[ii][ij]=3.8224719101123594 ij=15, u[ii][ij]=-0.6370786516853932 ij=1, u[ii][ij]=-83.86031417424024 ij=6, u[ii][ij]=22.435769743005128 ij=11, u[ii][ij]=15.623846495336752 ij=16, u[ii][ij]=-3.8874431053156697 ij=2, u[ii][ij]=-4.874308376922574 ij=7, u[ii][ij]=11.467415730337079 ij=12, u[ii][ij]=-3.8224719101123594 ij=17, u[ii][ij]=0.6370786516853932 ij=3, u[ii][ij]=0.7722628753054761 ij=8, u[ii][ij]=-1.4334269662921348 ij=13, u[ii][ij]=0.47780898876404493 ij=18, u[ii][ij]=-0.07963483146067415 cc=-66.26449653079774 ii=1, ut[ii][ii]=- 83.860, ut[ii][cs]= 66.541 working row edge i=3, j=1 working row 2, i=3, j=1, bi=0, bj=-1 a1=0.256708559516296, b1=1.4, c1=4.566937378852422 d1=-0.6666666666666667, e1=0.0 f1=-1.3333333333333335, c=-44.64628388616851 ij=0, u[ii][ij]=-1.2341980863038888 ij=5, u[ii][ij]=1.575 ij=10, u[ii][ij]=-0.5249999999999999 ij=15, u[ii][ij]=0.0875 ij=1, u[ii][ij]=11.413836047528886 ij=6, u[ii][ij]=-12.6 ij=11, u[ii][ij]=4.199999999999999 ij=16, u[ii][ij]=-0.7 ij=2, u[ii][ij]=-75.61333660523631 ij=7, u[ii][ij]=20.551218204835898 ij=12, u[ii][ij]=13.700812136557264 ij=17, u[ii][ij]=-3.425203034139316 ij=3, u[ii][ij]=-5.252830619137781 ij=8, u[ii][ij]=12.6 ij=13, u[ii][ij]=-4.199999999999999 ij=18, u[ii][ij]=0.7 ij=4, u[ii][ij]=0.8491352470294447 ij=9, u[ii][ij]=-1.575 ij=14, u[ii][ij]=0.5249999999999999 ij=19, u[ii][ij]=-0.0875 cc=-44.64628388616851 ii=2, ut[ii][ii]=- 75.613, ut[ii][cs]=- 1.719 working row edge i=4, j=1 working row 3, i=4, j=1, bi=0, bj=-1 a1=0.30326532985631666, b1=1.2741573033707865, c1=3.775712528680243 d1=-0.5555555555555557, e1=0.14814814814814817 f1=-0.3333333333333337, c=-20.414307323516862 cc=-20.414307323516862 ii=3, ut[ii][ii]=- 64.163, ut[ii][cs]=- 23.770 build_coef for point ix=2, iy=3 deriv pdxx returns bb=12, aa=-1 16 -30 16 -1 deriv pdyy returns bb=12, aa=-1 4 6 -20 11 deriv pdx returns bb=12, aa=1 -8 0 8 -1 deriv pdy returns bb=12, aa=-1 6 -18 10 3 pdxx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 - 0.750 0.000 i=1 0.000 0.000 0.000 12.000 0.000 i=2 0.000 0.000 0.000 - 22.500 0.000 i=3 0.000 0.000 0.000 12.000 0.000 i=4 0.000 0.000 0.000 - 0.750 0.000 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0 - 0.062 0.375 - 1.125 0.625 0.188 i=1 0.500 - 3.000 9.000 - 5.000 - 1.500 i=2 - 0.000 0.000 - 0.000 0.000 0.000 i=3 - 0.500 3.000 - 9.000 5.000 1.500 i=4 0.062 - 0.375 1.125 - 0.625 - 0.188 pdyy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 - 0.750 3.000 4.500 - 15.000 8.250 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 pdx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.250 0.000 i=1 0.000 0.000 0.000 - 2.000 0.000 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.000 0.000 0.000 2.000 0.000 i=4 0.000 0.000 0.000 - 0.250 0.000 pdy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 - 0.250 1.500 - 4.500 2.500 0.750 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 working row edge i=2, j=5 working row 21, i=2, j=5, bi=0, bj=1 a1=0.8243606353500639, b1=1.2741573033707865, c1=3.775712528680243 d1=0.7777777777777777, e1=-0.1481481481481481 f1=0.33333333333333304, c=42.54374998055788 cc=42.54374998055788 ii=21, ut[ii][ii]=- 75.221, ut[ii][cs]=- 350.310 working row edge i=3, j=5 working row 22, i=3, j=5, bi=0, bj=1 a1=0.9738670205273378, b1=1.4, c1=3.2087751498278214 d1=0.6666666666666665, e1=0.0 f1=1.333333333333333, c=69.118755700485 cc=69.118755700485 ii=22, ut[ii][ii]=- 68.710, ut[ii][cs]=- 396.662 working row edge i=4, j=5 working row 23, i=4, j=5, bi=0, bj=1 a1=1.150487945446412, b1=1.2741573033707865, c1=2.417550299655643 d1=0.7777777777777776, e1=0.14814814814814806 f1=2.333333333333333, c=100.51996394802688 ij=6, u[ii][ij]=-0.07963483146067415 ij=11, u[ii][ij]=0.47780898876404493 ij=16, u[ii][ij]=-1.4334269662921348 ij=21, u[ii][ij]=0.12792679996637688 ij=7, u[ii][ij]=0.6370786516853932 ij=12, u[ii][ij]=-3.8224719101123594 ij=17, u[ii][ij]=11.467415730337079 ij=22, u[ii][ij]=5.879513272947458 ij=8, u[ii][ij]=-1.850199761778769 ij=13, u[ii][ij]=7.474873121189151 ij=18, u[ii][ij]=10.212309681783728 ij=23, u[ii][ij]=-59.445529563675215 ij=9, u[ii][ij]=-0.6370786516853932 ij=14, u[ii][ij]=3.8224719101123594 ij=19, u[ii][ij]=-11.467415730337079 ij=24, u[ii][ij]=21.732197417766432 cc=100.51996394802688 ii=23, ut[ii][ii]=- 59.446, ut[ii][cs]=- 337.037 four sides, one element in, special case, generated build_coef for point ix=1, iy=1 deriv pdxx returns bb=12, aa=11 -20 6 4 -1 deriv pdyy returns bb=12, aa=11 -20 6 4 -1 deriv pdx returns bb=12, aa=-3 -10 18 -6 1 deriv pdy returns bb=12, aa=-3 -10 18 -6 1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 8.250 0.000 0.000 0.000 i=1 0.000 - 15.000 0.000 0.000 0.000 i=2 0.000 4.500 0.000 0.000 0.000 i=3 0.000 3.000 0.000 0.000 0.000 i=4 0.000 - 0.750 0.000 0.000 0.000 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0 0.562 1.875 - 3.375 1.125 - 0.188 i=1 1.875 6.250 - 11.250 3.750 - 0.625 i=2 - 3.375 - 11.250 20.250 - 6.750 1.125 i=3 1.125 3.750 - 6.750 2.250 - 0.375 i=4 - 0.188 - 0.625 1.125 - 0.375 0.062 pdyy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 8.250 - 15.000 4.500 3.000 - 0.750 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 pdx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 - 0.750 0.000 0.000 0.000 i=1 0.000 - 2.500 0.000 0.000 0.000 i=2 0.000 4.500 0.000 0.000 0.000 i=3 0.000 - 1.500 0.000 0.000 0.000 i=4 0.000 0.250 0.000 0.000 0.000 pdy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 - 0.750 - 2.500 4.500 - 1.500 0.250 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 working row 0, i=1, j=1, bi=-1, bj=-1 a1=0.18393972058572114, b1=1.0035398230088493, c1=5.540103140787237 d1=-0.2222222222222222, e1=-0.2962962962962964 f1=-3.3333333333333335, c=-80.19852027406424 ij=0, u[ii][ij]=-81.6255560638261 ij=5, u[ii][ij]=12.30730779135968 ij=10, u[ii][ij]=20.82802820308934 ij=15, u[ii][ij]=-4.856363819045033 ij=1, u[ii][ij]=-11.462094266213809 ij=6, u[ii][ij]=20.3216814159292 ij=11, u[ii][ij]=-6.773893805309733 ij=16, u[ii][ij]=1.1289823008849553 ij=2, u[ii][ij]=4.648426831373682 ij=7, u[ii][ij]=-6.773893805309733 ij=12, u[ii][ij]=2.2579646017699107 ij=17, u[ii][ij]=-0.37632743362831844 ij=3, u[ii][ij]=-0.8207227353753772 ij=8, u[ii][ij]=1.1289823008849553 ij=13, u[ii][ij]=-0.37632743362831844 ij=18, u[ii][ij]=0.06272123893805308 cc=-80.19852027406424 ii=0, ut[ii][ii]=- 81.626, ut[ii][cs]= 213.508 build_coef for point ix=1, iy=3 deriv pdxx returns bb=12, aa=11 -20 6 4 -1 deriv pdyy returns bb=12, aa=-1 4 6 -20 11 deriv pdx returns bb=12, aa=-3 -10 18 -6 1 deriv pdy returns bb=12, aa=-1 6 -18 10 3 pdxx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 8.250 0.000 i=1 0.000 0.000 0.000 - 15.000 0.000 i=2 0.000 0.000 0.000 4.500 0.000 i=3 0.000 0.000 0.000 3.000 0.000 i=4 0.000 0.000 0.000 - 0.750 0.000 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0 0.188 - 1.125 3.375 - 1.875 - 0.562 i=1 0.625 - 3.750 11.250 - 6.250 - 1.875 i=2 - 1.125 6.750 - 20.250 11.250 3.375 i=3 0.375 - 2.250 6.750 - 3.750 - 1.125 i=4 - 0.062 0.375 - 1.125 0.625 0.188 pdyy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 - 0.750 3.000 4.500 - 15.000 8.250 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 pdx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 - 0.750 0.000 i=1 0.000 0.000 0.000 - 2.500 0.000 i=2 0.000 0.000 0.000 4.500 0.000 i=3 0.000 0.000 0.000 - 1.500 0.000 i=4 0.000 0.000 0.000 0.250 0.000 pdy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 - 0.250 1.500 - 4.500 2.500 0.750 i=2 0.000 0.000 0.000 0.000 0.000 i=3 0.000 0.000 0.000 0.000 0.000 i=4 0.000 0.000 0.000 0.000 0.000 working row 20, i=1, j=5, bi=-1, bj=1 a1=0.6978062125430445, b1=1.0035398230088495, c1=4.181940911762637 d1=1.1111111111111112, e1=-0.2962962962962962 f1=-0.666666666666667, c=22.857388882439523 cc=22.857388882439523 ii=20, ut[ii][ii]=- 83.654, ut[ii][cs]=- 277.927 build_coef for point ix=3, iy=1 deriv pdxx returns bb=12, aa=-1 4 6 -20 11 deriv pdyy returns bb=12, aa=11 -20 6 4 -1 deriv pdx returns bb=12, aa=-1 6 -18 10 3 deriv pdy returns bb=12, aa=-3 -10 18 -6 1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 - 0.750 0.000 0.000 0.000 i=1 0.000 3.000 0.000 0.000 0.000 i=2 0.000 4.500 0.000 0.000 0.000 i=3 0.000 - 15.000 0.000 0.000 0.000 i=4 0.000 8.250 0.000 0.000 0.000 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0 0.188 0.625 - 1.125 0.375 - 0.062 i=1 - 1.125 - 3.750 6.750 - 2.250 0.375 i=2 3.375 11.250 - 20.250 6.750 - 1.125 i=3 - 1.875 - 6.250 11.250 - 3.750 0.625 i=4 - 0.562 - 1.875 3.375 - 1.125 0.188 pdyy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 0.000 0.000 0.000 0.000 0.000 i=3 8.250 - 15.000 4.500 3.000 - 0.750 i=4 0.000 0.000 0.000 0.000 0.000 pdx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 - 0.250 0.000 0.000 0.000 i=1 0.000 1.500 0.000 0.000 0.000 i=2 0.000 - 4.500 0.000 0.000 0.000 i=3 0.000 2.500 0.000 0.000 0.000 i=4 0.000 0.750 0.000 0.000 0.000 pdy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 0.000 0.000 0.000 0.000 0.000 i=3 - 0.750 - 2.500 4.500 - 1.500 0.250 i=4 0.000 0.000 0.000 0.000 0.000 working row 4, i=5, j=1, bi=1, bj=-1 a1=0.3582656552868946, b1=1.0035398230088495, c1=2.67146929674307 d1=-0.2222222222222225, e1=0.29629629629629634 f1=0.6666666666666661, c=2.4651929136270514 cc=2.4651929136270514 ii=4, ut[ii][ii]=- 52.348, ut[ii][cs]=- 51.723 build_coef for point ix=3, iy=3 deriv pdxx returns bb=12, aa=-1 4 6 -20 11 deriv pdyy returns bb=12, aa=-1 4 6 -20 11 deriv pdx returns bb=12, aa=-1 6 -18 10 3 deriv pdy returns bb=12, aa=-1 6 -18 10 3 pdxx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 - 0.750 0.000 i=1 0.000 0.000 0.000 3.000 0.000 i=2 0.000 0.000 0.000 4.500 0.000 i=3 0.000 0.000 0.000 - 15.000 0.000 i=4 0.000 0.000 0.000 8.250 0.000 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0 0.062 - 0.375 1.125 - 0.625 - 0.188 i=1 - 0.375 2.250 - 6.750 3.750 1.125 i=2 1.125 - 6.750 20.250 - 11.250 - 3.375 i=3 - 0.625 3.750 - 11.250 6.250 1.875 i=4 - 0.188 1.125 - 3.375 1.875 0.562 pdyy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 0.000 0.000 0.000 0.000 0.000 i=3 - 0.750 3.000 4.500 - 15.000 8.250 i=4 0.000 0.000 0.000 0.000 0.000 pdx j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 - 0.250 0.000 i=1 0.000 0.000 0.000 1.500 0.000 i=2 0.000 0.000 0.000 - 4.500 0.000 i=3 0.000 0.000 0.000 2.500 0.000 i=4 0.000 0.000 0.000 0.750 0.000 pdy j=0 j=1 j=2 j=3 j=4 y i=0 0.000 0.000 0.000 0.000 0.000 i=1 0.000 0.000 0.000 0.000 0.000 i=2 0.000 0.000 0.000 0.000 0.000 i=3 - 0.250 1.500 - 4.500 2.500 0.750 i=4 0.000 0.000 0.000 0.000 0.000 working row 24, i=5, j=5, bi=1, bj=1 a1=1.359140914229522, b1=1.00353982300885, c1=1.31330706771847 d1=1.1111111111111107, e1=0.2962962962962961 f1=3.3333333333333326, c=137.6980001654726 ij=6, u[ii][ij]=0.06272123893805312 ij=11, u[ii][ij]=-0.3763274336283187 ij=16, u[ii][ij]=1.1289823008849562 ij=21, u[ii][ij]=-1.9243458528304505 ij=7, u[ii][ij]=-0.3763274336283187 ij=12, u[ii][ij]=2.2579646017699124 ij=17, u[ii][ij]=-6.773893805309737 ij=22, u[ii][ij]=9.507363745638418 ij=8, u[ii][ij]=1.1289823008849562 ij=13, u[ii][ij]=-6.773893805309737 ij=18, u[ii][ij]=20.32168141592921 ij=23, u[ii][ij]=-10.17368889481671 ij=9, u[ii][ij]=-1.6862667642434577 ij=14, u[ii][ij]=8.14763998388304 ij=19, u[ii][ij]=-6.713274537449779 ij=24, u[ii][ij]=-26.96274398356272 cc=137.6980001654726 ii=24, ut[ii][ii]=- 26.963, ut[ii][cs]=- 520.864 matrix initialized initial matrix i=0, j=0 .. 25 - 81.626 - 11.462 4.648 - 0.821 0.000 12.307 20.322 - 6.774 1.129 0.000 20.828 - 6.774 2.258 - 0.376 0.000 - 4.856 1.129 - 0.376 0.063 0.000 0.000 0.000 0.000 0.000 0.000 213.508 i=1, j=0 .. 25 10.089 - 83.860 - 4.874 0.772 0.000 - 11.467 22.436 11.467 - 1.433 0.000 3.822 15.624 - 3.822 0.478 0.000 - 0.637 - 3.887 0.637 - 0.080 0.000 0.000 0.000 0.000 0.000 0.000 66.541 i=2, j=0 .. 25 - 1.234 11.414 - 75.613 - 5.253 0.849 1.575 - 12.600 20.551 12.600 - 1.575 - 0.525 4.200 13.701 - 4.200 0.525 0.088 - 0.700 - 3.425 0.700 - 0.088 0.000 0.000 0.000 0.000 0.000 - 1.719 i=3, j=0 .. 25 0.000 - 1.163 11.121 - 64.163 - 3.843 0.000 1.433 - 11.467 17.657 11.467 0.000 - 0.478 3.822 11.105 - 3.822 0.000 0.080 - 0.637 - 2.795 0.637 0.000 0.000 0.000 0.000 0.000 - 23.770 i=4, j=0 .. 25 0.000 0.414 - 3.022 13.902 - 52.348 0.000 - 1.129 6.774 - 20.322 24.645 0.000 0.376 - 2.258 6.774 3.807 0.000 - 0.063 0.376 - 1.129 - 1.302 0.000 0.000 0.000 0.000 0.000 - 51.723 i=5, j=0 .. 25 69.881 - 11.467 3.822 - 0.637 0.000 - 125.600 1.655 0.603 - 0.165 0.000 56.843 11.467 - 3.822 0.637 0.000 - 3.145 - 1.433 0.478 - 0.080 0.000 0.000 0.000 0.000 0.000 0.000 - 71.367 i=6, j=0 .. 25 5.465 58.562 - 5.465 0.683 0.000 4.084 - 118.155 3.195 - 0.172 0.000 - 5.465 58.414 5.465 - 0.683 0.000 0.683 - 3.646 - 0.683 0.085 0.000 0.000 0.000 0.000 0.000 0.000 - 58.892 i=7, j=0 .. 25 - 0.700 5.600 51.684 - 5.600 0.700 - 0.352 4.966 - 105.636 3.633 - 0.185 0.700 - 5.600 51.684 5.600 - 0.700 - 0.088 0.700 - 3.230 - 0.700 0.088 0.000 0.000 0.000 0.000 0.000 - 24.979 i=8, j=0 .. 25 0.000 - 0.683 5.465 42.116 - 5.465 0.000 - 0.373 5.523 - 88.295 4.634 0.000 0.683 - 5.465 42.264 5.465 0.000 - 0.085 0.683 - 2.646 - 0.683 0.000 0.000 0.000 0.000 0.000 9.761 i=9, j=0 .. 25 0.000 0.637 - 3.822 11.467 22.420 0.000 - 0.403 1.667 1.750 - 60.149 0.000 - 0.637 3.822 - 11.467 35.458 0.000 0.080 - 0.478 1.433 - 2.624 0.000 0.000 0.000 0.000 0.000 - 28.906 i=10, j=0 .. 25 - 4.605 1.575 - 0.525 0.088 0.000 66.678 - 12.600 4.200 - 0.700 0.000 - 120.381 3.612 0.408 - 0.158 0.000 52.678 12.600 - 4.200 0.700 0.000 - 2.855 - 1.575 0.525 - 0.088 0.000 - 28.077 i=11, j=0 .. 25 - 0.700 - 3.425 0.700 - 0.088 0.000 5.600 54.803 - 5.600 0.700 0.000 4.857 - 113.279 5.301 - 0.345 0.000 - 5.600 54.803 5.600 - 0.700 0.000 0.700 - 3.425 - 0.700 0.088 0.000 - 15.545 i=12, j=0 .. 25 0.088 - 0.700 - 3.000 0.700 - 0.088 - 0.700 5.600 48.000 - 5.600 0.700 - 0.375 6.000 - 101.250 6.000 - 0.375 0.700 - 5.600 48.000 5.600 - 0.700 - 0.088 0.700 - 3.000 - 0.700 0.088 7.000 i=13, j=0 .. 25 0.000 0.088 - 0.700 - 2.407 0.700 0.000 - 0.700 5.600 38.505 - 5.600 0.000 - 0.415 6.866 - 84.488 7.310 0.000 0.700 - 5.600 38.505 5.600 0.000 - 0.088 0.700 - 2.407 - 0.700 42.591 i=14, j=0 .. 25 0.000 - 0.088 0.525 - 1.575 - 0.703 0.000 0.700 - 4.200 12.600 18.254 0.000 - 0.634 2.760 1.140 - 54.708 0.000 - 0.700 4.200 - 12.600 32.254 0.000 0.088 - 0.525 1.575 - 2.453 - 52.636 i=15, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 - 4.273 1.433 - 0.478 0.080 0.000 61.844 - 11.467 3.822 - 0.637 0.000 - 114.513 5.750 0.333 - 0.181 0.000 48.806 11.467 - 3.822 0.637 0.000 23.941 i=16, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 - 0.683 - 3.162 0.683 - 0.085 0.000 5.465 50.525 - 5.465 0.683 0.000 6.199 - 108.219 7.977 - 0.554 0.000 - 5.465 50.377 5.465 - 0.683 0.000 56.619 i=17, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.088 - 0.700 - 2.728 0.700 - 0.088 - 0.700 5.600 43.647 - 5.600 0.700 - 0.440 7.707 - 96.873 9.040 - 0.607 0.700 - 5.600 43.647 5.600 - 0.700 87.945 i=18, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.000 0.085 - 0.683 - 2.125 0.683 0.000 - 0.683 5.465 34.079 - 5.465 0.000 - 0.507 9.003 - 80.918 10.781 0.000 0.683 - 5.465 34.227 5.465 140.198 i=19, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.000 - 0.080 0.478 - 1.433 - 0.491 0.000 0.637 - 3.822 11.467 14.383 0.000 - 0.925 4.088 0.882 - 49.188 0.000 - 0.637 3.822 - 11.467 27.421 - 67.994 i=20, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 - 2.435 - 1.129 0.376 - 0.063 0.000 8.338 6.774 - 2.258 0.376 0.000 31.442 - 20.322 6.774 - 1.129 0.000 - 83.654 19.430 - 3.337 0.382 0.000 - 277.927 i=21, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.637 - 2.795 - 0.637 0.080 0.000 - 3.822 11.105 3.822 - 0.478 0.000 11.467 17.657 - 11.467 1.433 0.000 1.966 - 75.221 17.819 - 1.609 0.000 - 350.310 i=22, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 - 0.088 0.700 - 2.407 - 0.700 0.088 0.525 - 4.200 9.626 4.200 - 0.525 - 1.575 12.600 14.439 - 12.600 1.575 0.311 3.353 - 68.710 20.020 - 1.772 - 396.662 i=23, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.000 - 0.080 0.637 - 1.850 - 0.637 0.000 0.478 - 3.822 7.475 3.822 0.000 - 1.433 11.467 10.212 - 11.467 0.000 0.128 5.880 - 59.446 21.732 - 337.037 i=24, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.000 0.063 - 0.376 1.129 - 1.686 0.000 - 0.376 2.258 - 6.774 8.148 0.000 1.129 - 6.774 20.322 - 6.713 0.000 - 1.924 9.507 - 10.174 - 26.963 - 520.864 computed solution us[x,y] row i=1, at x=-0.667, first y=-0.667 - 1.407 1.667 3.704 5.148 6.444 row i=2, at x=-0.333, first y=-0.667 1.000 3.852 5.963 7.778 9.741 row i=3, at x=0.000, first y=-0.667 2.741 5.593 8.000 10.407 13.259 row i=4, at x=0.333, first y=-0.667 4.037 7.111 10.037 13.259 17.222 row i=5, at x=0.667, first y=-0.667 5.111 8.630 12.296 16.556 21.852 exact solution u(x,y) - 1.407 1.667 3.704 5.148 6.444 1.000 3.852 5.963 7.778 9.741 2.741 5.593 8.000 10.407 13.259 4.037 7.111 10.037 13.259 17.222 5.111 8.630 12.296 16.556 21.852 difference exact-computed solution u(x,y)-us[x,y] - 0.000 - 0.000 - 0.000 0.000 0.000 0.000 - 0.000 - 0.000 0.000 0.000 - 0.000 - 0.000 - 0.000 - 0.000 0.000 - 0.000 - 0.000 - 0.000 - 0.000 - 0.000 - 0.000 - 0.000 - 0.000 - 0.000 - 0.000 total error=1.3489209749195652E-13 average error=5.395683899678261E-15 max error=1.5987211554602254E-14