pdenu23_eq.c running xmin=-1, xmax=1.2, nx=11 ymin=-1.1, ymax=1.3, ny=10 zmin=-1.2, zmax=1.4, nz=9 u(xg[1],yg[1],zg[1])=-12.069 c(xg[1],yg[1],zg[1])=73.84 internal cells zeroed matrix initialized initial matrix, left side, upper 331.412 -932.608 275.439 -117.383 69.001 -43.298 36.844 593.812 -555.946 -192.780 70.387 -39.272 23.939 -20.080 -233.368 506.623 -367.260 -159.044 73.162 -40.775 32.843 143.525 -240.382 307.236 -216.555 -180.181 76.109 -55.559 -105.947 165.489 -165.180 283.185 -263.037 -152.156 89.830 90.284 -135.977 122.394 -148.410 248.567 -185.827 -220.404 -70.547 104.417 -89.875 97.472 -125.618 246.865 -148.606 initial matrix, right side, upper 0.000 0.000 0.000 0.000 0.000 0.000 4664.242 0.000 0.000 0.000 0.000 0.000 0.000 3641.459 0.000 0.000 0.000 0.000 0.000 0.000 3201.142 -0.077 0.000 0.000 0.000 0.000 0.000 2933.681 0.000 -0.077 0.000 0.000 0.000 0.000 2752.386 0.000 0.000 -0.077 0.000 0.000 0.000 2838.124 0.000 0.000 0.000 -0.077 0.000 0.000 3016.147 initial matrix, left side, lower 0.000 0.000 -19.912 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -19.912 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -19.912 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -19.912 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -19.912 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 initial matrix, right side, lower 96.938 -38.223 20.090 -15.791 11.240 -4.152 -7271.281 -334.622 91.969 -32.567 22.618 -15.330 5.510 -8033.692 -76.166 -310.081 79.761 -41.464 25.727 -8.843 -9467.845 43.284 -61.246 -340.333 108.388 -49.425 15.056 -11038.563 -27.279 33.053 -50.580 -380.070 135.654 -27.564 -12891.364 35.344 -40.659 56.036 -111.216 -337.180 114.169 -14142.121 -123.230 137.923 -178.958 327.293 -511.089 -37.924 -15834.210 exact solution xg[1]=-0.900000,yg[1]=-1.000000, zg[1]=-1.000000, u(x,y,z)=-12.069 xg[1]=-0.900000,yg[1]=-1.000000, zg[7]=1.200000, u(x,y,z)=22.339 xg[1]=-0.900000,yg[8]=1.200000, zg[1]=-1.000000, u(x,y,z)=-15.545 xg[1]=-0.900000,yg[8]=1.200000, zg[7]=1.200000, u(x,y,z)=23.703 xg[9]=1.000000,yg[1]=-1.000000, zg[1]=-1.000000, u(x,y,z)=-13 xg[9]=1.000000,yg[1]=-1.000000, zg[7]=1.200000, u(x,y,z)=26.424 xg[9]=1.000000,yg[8]=1.200000, zg[1]=-1.000000, u(x,y,z)=14.456 xg[9]=1.000000,yg[8]=1.200000, zg[7]=1.200000, u(x,y,z)=58.72 computed solution xg[1]=-0.900000,yg[1]=-1.000000, zg[1]=-1.000000, us=-12.069 xg[1]=-0.900000,yg[1]=-1.000000, zg[7]=1.200000, us=22.339 xg[1]=-0.900000,yg[8]=1.200000, zg[1]=-1.000000, us=-15.545 xg[1]=-0.900000,yg[8]=1.200000, zg[7]=1.200000, us=23.703 xg[9]=1.000000,yg[1]=-1.000000, zg[1]=-1.000000, us=-13 xg[9]=1.000000,yg[1]=-1.000000, zg[7]=1.200000, us=26.424 xg[9]=1.000000,yg[8]=1.200000, zg[1]=-1.000000, us=14.456 xg[9]=1.000000,yg[8]=1.200000, zg[7]=1.200000, us=58.72 difference exact-computed solution xg[1]=-0.900000,yg[1]=-1.000000, zg[1]=-1.000000, diff(x,y,z)=-2.71321e-11 xg[1]=-0.900000,yg[1]=-1.000000, zg[7]=1.200000, diff(x,y,z)=-4.83666e-11 xg[1]=-0.900000,yg[8]=1.200000, zg[1]=-1.000000, diff(x,y,z)=4.80682e-12 xg[1]=-0.900000,yg[8]=1.200000, zg[7]=1.200000, diff(x,y,z)=1.08287e-11 xg[9]=1.000000,yg[1]=-1.000000, zg[1]=-1.000000, diff(x,y,z)=-8.9976e-11 xg[9]=1.000000,yg[1]=-1.000000, zg[7]=1.200000, diff(x,y,z)=-1.01124e-10 xg[9]=1.000000,yg[8]=1.200000, zg[1]=-1.000000, diff(x,y,z)=2.16005e-12 xg[9]=1.000000,yg[8]=1.200000, zg[7]=1.200000, diff(x,y,z)=-5.80371e-11 total error=3.85759e-08, average error=7.65394e-11, max error=6.40306e-10