pde_abc_eq.c running with user files abc.h abc.c xmin=-1, xmax=1, nx=7, hx=0.333333 ymin=-1, ymax=1, ny=7, hy=0.333333 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.00 0.00 -0.75 0.00 0.00 i=1, 0.00 0.00 12.00 0.00 0.00 i=2, 0.00 0.00 -22.50 0.00 0.00 i=3, 0.00 0.00 12.00 0.00 0.00 i=4, 0.00 0.00 -0.75 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.06 -0.50 0.00 0.50 -0.06 i=1, -0.50 4.00 -0.00 -4.00 0.50 i=2, 0.00 -0.00 0.00 0.00 -0.00 i=3, 0.50 -4.00 0.00 4.00 -0.50 i=4, -0.06 0.50 -0.00 -0.50 0.06 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, -0.75 12.00 -22.50 12.00 -0.75 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.25 0.00 0.00 i=1, 0.00 0.00 -2.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 2.00 0.00 0.00 i=4, 0.00 0.00 -0.25 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.25 -2.00 0.00 2.00 -0.25 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row 6, i=2, j=2, bi=0, bj=0 ii=6, ut[ii][ii]=-118.155, ut[ii][cs]=-58.8925 working row 11, i=2, j=3, bi=0, bj=0 ii=11, ut[ii][ii]=-113.279, ut[ii][cs]=-15.5446 working row 16, i=2, j=4, bi=0, bj=0 ii=16, ut[ii][ii]=-108.219, ut[ii][cs]=56.6188 working row 7, i=3, j=2, bi=0, bj=0 ii=7, ut[ii][ii]=-105.636, ut[ii][cs]=-24.9793 working row 12, i=3, j=3, bi=0, bj=0 ii=12, ut[ii][ii]=-101.25, ut[ii][cs]=7 working row 17, i=3, j=4, bi=0, bj=0 ii=17, ut[ii][ii]=-96.8728, ut[ii][cs]=87.9453 working row 8, i=4, j=2, bi=0, bj=0 ii=8, ut[ii][ii]=-88.2953, ut[ii][cs]=9.7605 working row 13, i=4, j=3, bi=0, bj=0 ii=13, ut[ii][ii]=-84.4877, ut[ii][cs]=42.5907 working row 18, i=4, j=4, bi=0, bj=0 ii=18, ut[ii][ii]=-80.9177, 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.00 0.00 8.25 0.00 0.00 i=1, 0.00 0.00 -15.00 0.00 0.00 i=2, 0.00 0.00 4.50 0.00 0.00 i=3, 0.00 0.00 3.00 0.00 0.00 i=4, 0.00 0.00 -0.75 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, -0.19 1.50 -0.00 -1.50 0.19 i=1, -0.62 5.00 -0.00 -5.00 0.62 i=2, 1.12 -9.00 0.00 9.00 -1.12 i=3, -0.38 3.00 -0.00 -3.00 0.38 i=4, 0.06 -0.50 0.00 0.50 -0.06 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, -0.75 12.00 -22.50 12.00 -0.75 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 -0.75 0.00 0.00 i=1, 0.00 0.00 -2.50 0.00 0.00 i=2, 0.00 0.00 4.50 0.00 0.00 i=3, 0.00 0.00 -1.50 0.00 0.00 i=4, 0.00 0.00 0.25 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.25 -2.00 0.00 2.00 -0.25 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row edge i=1, j=2 working row 5, i=1, j=2, bi=-1, bj=0 ii=5, ut[ii][ii]=-125.6, ut[ii][cs]=-71.3669 working row edge i=1, j=3 working row 10, i=1, j=3, bi=-1, bj=0 ii=10, ut[ii][ii]=-120.381, ut[ii][cs]=-28.0771 working row edge i=1, j=4 working row 15, i=1, j=4, bi=-1, bj=0 ii=15, ut[ii][ii]=-114.513, ut[ii][cs]=23.9413 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.00 0.00 -0.75 0.00 0.00 i=1, 0.00 0.00 3.00 0.00 0.00 i=2, 0.00 0.00 4.50 0.00 0.00 i=3, 0.00 0.00 -15.00 0.00 0.00 i=4, 0.00 0.00 8.25 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, -0.06 0.50 -0.00 -0.50 0.06 i=1, 0.38 -3.00 0.00 3.00 -0.38 i=2, -1.12 9.00 -0.00 -9.00 1.12 i=3, 0.62 -5.00 0.00 5.00 -0.62 i=4, 0.19 -1.50 0.00 1.50 -0.19 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, -0.75 12.00 -22.50 12.00 -0.75 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 -0.25 0.00 0.00 i=1, 0.00 0.00 1.50 0.00 0.00 i=2, 0.00 0.00 -4.50 0.00 0.00 i=3, 0.00 0.00 2.50 0.00 0.00 i=4, 0.00 0.00 0.75 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.25 -2.00 0.00 2.00 -0.25 i=4, 0.00 0.00 0.00 0.00 0.00 working row edge i=5, j=2 working row 9, i=5, j=2, bi=1, bj=0 ii=9, ut[ii][ii]=-60.1492, ut[ii][cs]=-28.9059 working row edge i=5, j=3 working row 14, i=5, j=3, bi=1, bj=0 ii=14, ut[ii][ii]=-54.708, ut[ii][cs]=-52.6363 working row edge i=5, j=4 working row 19, i=5, j=4, bi=1, bj=0 ii=19, ut[ii][ii]=-49.1876, ut[ii][cs]=-67.9943 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.00 -0.75 0.00 0.00 0.00 i=1, 0.00 12.00 0.00 0.00 0.00 i=2, 0.00 -22.50 0.00 0.00 0.00 i=3, 0.00 12.00 0.00 0.00 0.00 i=4, 0.00 -0.75 0.00 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, -0.19 -0.62 1.12 -0.38 0.06 i=1, 1.50 5.00 -9.00 3.00 -0.50 i=2, -0.00 -0.00 0.00 -0.00 0.00 i=3, -1.50 -5.00 9.00 -3.00 0.50 i=4, 0.19 0.62 -1.12 0.38 -0.06 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 8.25 -15.00 4.50 3.00 -0.75 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.25 0.00 0.00 0.00 i=1, 0.00 -2.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 2.00 0.00 0.00 0.00 i=4, 0.00 -0.25 0.00 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, -0.75 -2.50 4.50 -1.50 0.25 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row edge i=2, j=1 working row 1, i=2, j=1, bi=0, bj=-1 ii=1, ut[ii][ii]=-83.8603, ut[ii][cs]=66.5411 working row edge i=3, j=1 working row 2, i=3, j=1, bi=0, bj=-1 ii=2, ut[ii][ii]=-75.6133, ut[ii][cs]=-1.71905 working row edge i=4, j=1 working row 3, i=4, j=1, bi=0, bj=-1 ii=3, ut[ii][ii]=-64.1629, ut[ii][cs]=-23.7704 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.00 0.00 0.00 -0.75 0.00 i=1, 0.00 0.00 0.00 12.00 0.00 i=2, 0.00 0.00 0.00 -22.50 0.00 i=3, 0.00 0.00 0.00 12.00 0.00 i=4, 0.00 0.00 0.00 -0.75 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, -0.06 0.38 -1.12 0.62 0.19 i=1, 0.50 -3.00 9.00 -5.00 -1.50 i=2, -0.00 0.00 -0.00 0.00 0.00 i=3, -0.50 3.00 -9.00 5.00 1.50 i=4, 0.06 -0.38 1.12 -0.62 -0.19 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, -0.75 3.00 4.50 -15.00 8.25 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.25 0.00 i=1, 0.00 0.00 0.00 -2.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 2.00 0.00 i=4, 0.00 0.00 0.00 -0.25 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, -0.25 1.50 -4.50 2.50 0.75 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row edge i=2, j=5 working row 21, i=2, j=5, bi=0, bj=1 ii=21, ut[ii][ii]=-75.2208, ut[ii][cs]=-350.31 working row edge i=3, j=5 working row 22, i=3, j=5, bi=0, bj=1 ii=22, ut[ii][ii]=-68.7103, ut[ii][cs]=-396.662 working row edge i=4, j=5 working row 23, i=4, j=5, bi=0, bj=1 ii=23, ut[ii][ii]=-59.4455, 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.00 8.25 0.00 0.00 0.00 i=1, 0.00 -15.00 0.00 0.00 0.00 i=2, 0.00 4.50 0.00 0.00 0.00 i=3, 0.00 3.00 0.00 0.00 0.00 i=4, 0.00 -0.75 0.00 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.56 1.88 -3.38 1.12 -0.19 i=1, 1.88 6.25 -11.25 3.75 -0.62 i=2, -3.38 -11.25 20.25 -6.75 1.12 i=3, 1.12 3.75 -6.75 2.25 -0.38 i=4, -0.19 -0.62 1.12 -0.38 0.06 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 8.25 -15.00 4.50 3.00 -0.75 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 -0.75 0.00 0.00 0.00 i=1, 0.00 -2.50 0.00 0.00 0.00 i=2, 0.00 4.50 0.00 0.00 0.00 i=3, 0.00 -1.50 0.00 0.00 0.00 i=4, 0.00 0.25 0.00 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, -0.75 -2.50 4.50 -1.50 0.25 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row 0, i=1, j=1, bi=-1, bj=-1 ii=0, ut[ii][ii]=-81.6256, 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.00 0.00 0.00 8.25 0.00 i=1, 0.00 0.00 0.00 -15.00 0.00 i=2, 0.00 0.00 0.00 4.50 0.00 i=3, 0.00 0.00 0.00 3.00 0.00 i=4, 0.00 0.00 0.00 -0.75 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.19 -1.12 3.38 -1.88 -0.56 i=1, 0.62 -3.75 11.25 -6.25 -1.88 i=2, -1.12 6.75 -20.25 11.25 3.38 i=3, 0.38 -2.25 6.75 -3.75 -1.12 i=4, -0.06 0.38 -1.12 0.62 0.19 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, -0.75 3.00 4.50 -15.00 8.25 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 -0.75 0.00 i=1, 0.00 0.00 0.00 -2.50 0.00 i=2, 0.00 0.00 0.00 4.50 0.00 i=3, 0.00 0.00 0.00 -1.50 0.00 i=4, 0.00 0.00 0.00 0.25 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, -0.25 1.50 -4.50 2.50 0.75 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row 20, i=1, j=5, bi=-1, bj=1 ii=20, ut[ii][ii]=-83.6535, 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.00 -0.75 0.00 0.00 0.00 i=1, 0.00 3.00 0.00 0.00 0.00 i=2, 0.00 4.50 0.00 0.00 0.00 i=3, 0.00 -15.00 0.00 0.00 0.00 i=4, 0.00 8.25 0.00 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.19 0.62 -1.12 0.38 -0.06 i=1, -1.12 -3.75 6.75 -2.25 0.38 i=2, 3.38 11.25 -20.25 6.75 -1.12 i=3, -1.88 -6.25 11.25 -3.75 0.62 i=4, -0.56 -1.88 3.38 -1.12 0.19 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 8.25 -15.00 4.50 3.00 -0.75 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 -0.25 0.00 0.00 0.00 i=1, 0.00 1.50 0.00 0.00 0.00 i=2, 0.00 -4.50 0.00 0.00 0.00 i=3, 0.00 2.50 0.00 0.00 0.00 i=4, 0.00 0.75 0.00 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, -0.75 -2.50 4.50 -1.50 0.25 i=4, 0.00 0.00 0.00 0.00 0.00 working row 4, i=5, j=1, bi=1, bj=-1 ii=4, ut[ii][ii]=-52.3478, ut[ii][cs]=-51.7228 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.00 0.00 0.00 -0.75 0.00 i=1, 0.00 0.00 0.00 3.00 0.00 i=2, 0.00 0.00 0.00 4.50 0.00 i=3, 0.00 0.00 0.00 -15.00 0.00 i=4, 0.00 0.00 0.00 8.25 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.06 -0.38 1.12 -0.62 -0.19 i=1, -0.38 2.25 -6.75 3.75 1.12 i=2, 1.12 -6.75 20.25 -11.25 -3.38 i=3, -0.62 3.75 -11.25 6.25 1.88 i=4, -0.19 1.12 -3.38 1.88 0.56 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, -0.75 3.00 4.50 -15.00 8.25 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 -0.25 0.00 i=1, 0.00 0.00 0.00 1.50 0.00 i=2, 0.00 0.00 0.00 -4.50 0.00 i=3, 0.00 0.00 0.00 2.50 0.00 i=4, 0.00 0.00 0.00 0.75 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, -0.25 1.50 -4.50 2.50 0.75 i=4, 0.00 0.00 0.00 0.00 0.00 working row 24, i=5, j=5, bi=1, bj=1 ii=24, ut[ii][ii]=-26.9627, 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.087 -0.700 -3.425 0.700 -0.087 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.087 0.700 -3.230 -0.700 0.087 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.087 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.087 0.000 -28.077 i=11, j=0 .. 25 -0.700 -3.425 0.700 -0.087 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.087 0.000 -15.545 i=12, j=0 .. 25 0.087 -0.700 -3.000 0.700 -0.087 -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.087 0.700 -3.000 -0.700 0.087 7.000 i=13, j=0 .. 25 0.000 0.087 -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.087 0.700 -2.407 -0.700 42.591 i=14, j=0 .. 25 0.000 -0.087 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.087 -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.087 -0.700 -2.728 0.700 -0.087 -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.087 0.700 -2.407 -0.700 0.087 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] -1.41 1.67 3.70 5.15 6.44 1.00 3.85 5.96 7.78 9.74 2.74 5.59 8.00 10.41 13.26 4.04 7.11 10.04 13.26 17.22 5.11 8.63 12.30 16.56 21.85 check_soln values against PDE maxerr=5.58241e-13, rmserr=2.64364e-13, avgerr=2.09696e-13 exact solution u(x,y) -1.41 1.67 3.70 5.15 6.44 1.00 3.85 5.96 7.78 9.74 2.74 5.59 8.00 10.41 13.26 4.04 7.11 10.04 13.26 17.22 5.11 8.63 12.30 16.56 21.85 difference exact-computed solution u(x,y)-us[x,y] 0.00 0.00 0.00 0.00 0.00 -0.00 0.00 0.00 0.00 0.00 -0.00 -0.00 0.00 -0.00 0.00 -0.00 -0.00 -0.00 -0.00 0.00 -0.00 -0.00 -0.00 -0.00 -0.00 check against known analytic solution total error=5.55112e-14, average error=2.22045e-15, max error=7.10543e-15