Original array with boundary values and initial guesses of zero: i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 0.00000 0.00000 0.00000 11.00000 j = 2 2.00000 0.00000 0.00000 0.00000 10.00000 j = 1 1.00000 0.00000 0.00000 0.00000 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 * * * Calculations and results for iteration 1 using successive over relaxation * * * Calculating u(1,1) = -omega * [ AS(1,1) * u(1,0) + AW(1,1) * u(0,1) + AE(1,1) * uOld(2,1) + AN(1,1) * uOld(1,2) ] / AP(1,1) + (1 - omega) * uOld(1,1) u(1,1) = (1.17100) [ -(1.00000)(2.00000) - (1.00000)(1.00000) - (1.00000)(0.00000) - (1.00000)(0.00000) ] / (-4.00000) + (-0.17100)(0.00000) = 0.87825 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 0.00000 0.00000 0.00000 11.00000 j = 2 2.00000 0.00000 0.00000 0.00000 10.00000 j = 1 1.00000 0.87825 0.00000 0.00000 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,1) = -omega * [ AS(2,1) * u(2,0) + AW(2,1) * u(1,1) + AE(2,1) * uOld(3,1) + AN(2,1) * uOld(2,2) ] / AP(2,1) + (1 - omega) * uOld(2,1) u(2,1) = (1.17100) [ -(1.00000)(4.00000) - (1.00000)(0.87825) - (1.00000)(0.00000) - (1.00000)(0.00000) ] / (-4.00000) + (-0.17100)(0.00000) = 1.42811 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 0.00000 0.00000 0.00000 11.00000 j = 2 2.00000 0.00000 0.00000 0.00000 10.00000 j = 1 1.00000 0.87825 1.42811 0.00000 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,1) = -omega * [ AS(3,1) * u(3,0) + AW(3,1) * u(2,1) + AE(3,1) * uOld(4,1) + AN(3,1) * uOld(3,2) ] / AP(3,1) + (1 - omega) * uOld(3,1) u(3,1) = (1.17100) [ -(1.00000)(6.00000) - (1.00000)(1.42811) - (1.00000)(9.00000) - (1.00000)(0.00000) ] / (-4.00000) + (-0.17100)(0.00000) = 4.80933 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 0.00000 0.00000 0.00000 11.00000 j = 2 2.00000 0.00000 0.00000 0.00000 10.00000 j = 1 1.00000 0.87825 1.42811 4.80933 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,2) = -omega * [ AS(1,2) * u(1,1) + AW(1,2) * u(0,2) + AE(1,2) * uOld(2,2) + AN(1,2) * uOld(1,3) ] / AP(1,2) + (1 - omega) * uOld(1,2) u(1,2) = (1.17100) [ -(1.00000)(0.87825) - (1.00000)(2.00000) - (1.00000)(0.00000) - (1.00000)(0.00000) ] / (-4.00000) + (-0.17100)(0.00000) = 0.84261 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 0.00000 0.00000 0.00000 11.00000 j = 2 2.00000 0.84261 0.00000 0.00000 10.00000 j = 1 1.00000 0.87825 1.42811 4.80933 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,2) = -omega * [ AS(2,2) * u(2,1) + AW(2,2) * u(1,2) + AE(2,2) * uOld(3,2) + AN(2,2) * uOld(2,3) ] / AP(2,2) + (1 - omega) * uOld(2,2) u(2,2) = (1.17100) [ -(1.00000)(1.42811) - (1.00000)(0.84261) - (1.00000)(0.00000) - (1.00000)(0.00000) ] / (-4.00000) + (-0.17100)(0.00000) = 0.66475 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 0.00000 0.00000 0.00000 11.00000 j = 2 2.00000 0.84261 0.66475 0.00000 10.00000 j = 1 1.00000 0.87825 1.42811 4.80933 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,2) = -omega * [ AS(3,2) * u(3,1) + AW(3,2) * u(2,2) + AE(3,2) * uOld(4,2) + AN(3,2) * uOld(3,3) ] / AP(3,2) + (1 - omega) * uOld(3,2) u(3,2) = (1.17100) [ -(1.00000)(4.80933) - (1.00000)(0.66475) - (1.00000)(10.00000) - (1.00000)(0.00000) ] / (-4.00000) + (-0.17100)(0.00000) = 4.53004 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 0.00000 0.00000 0.00000 11.00000 j = 2 2.00000 0.84261 0.66475 4.53004 10.00000 j = 1 1.00000 0.87825 1.42811 4.80933 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,3) = -omega * [ AS(1,3) * u(1,2) + AW(1,3) * u(0,3) + AE(1,3) * uOld(2,3) + AN(1,3) * uOld(1,4) ] / AP(1,3) + (1 - omega) * uOld(1,3) u(1,3) = (1.17100) [ -(1.00000)(0.84261) - (1.00000)(3.00000) - (1.00000)(0.00000) - (1.00000)(6.00000) ] / (-4.00000) + (-0.17100)(0.00000) = 2.88142 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 2.88142 0.00000 0.00000 11.00000 j = 2 2.00000 0.84261 0.66475 4.53004 10.00000 j = 1 1.00000 0.87825 1.42811 4.80933 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,3) = -omega * [ AS(2,3) * u(2,2) + AW(2,3) * u(1,3) + AE(2,3) * uOld(3,3) + AN(2,3) * uOld(2,4) ] / AP(2,3) + (1 - omega) * uOld(2,3) u(2,3) = (1.17100) [ -(1.00000)(0.66475) - (1.00000)(2.88142) - (1.00000)(0.00000) - (1.00000)(8.00000) ] / (-4.00000) + (-0.17100)(0.00000) = 3.38014 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 2.88142 3.38014 0.00000 11.00000 j = 2 2.00000 0.84261 0.66475 4.53004 10.00000 j = 1 1.00000 0.87825 1.42811 4.80933 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,3) = -omega * [ AS(3,3) * u(3,2) + AW(3,3) * u(2,3) + AE(3,3) * uOld(4,3) + AN(3,3) * uOld(3,4) ] / AP(3,3) + (1 - omega) * uOld(3,3) u(3,3) = (1.17100) [ -(1.00000)(4.53004) - (1.00000)(3.38014) - (1.00000)(11.00000) - (1.00000)(10.00000) ] / (-4.00000) + (-0.17100)(0.00000) = 8.46346 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 2.88142 3.38014 8.46346 11.00000 j = 2 2.00000 0.84261 0.66475 4.53004 10.00000 j = 1 1.00000 0.87825 1.42811 4.80933 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 At end of 1 iterations, the maximum residual is 1.88047 and the maximum relative increment is 1.00000 Residual array i = 1 i = 2 i = 3 j = 3 0.42426 1.62226 1.23591 j = 2 0.76350 1.88047 1.45435 j = 1 0.43943 1.15997 0.43021 Relative increment array i = 1 i = 2 i = 3 j = 3 1.00000 1.00000 1.00000 j = 2 1.00000 1.00000 1.00000 j = 1 1.00000 1.00000 1.00000 * * * Calculations and results for iteration 2 using successive over relaxation * * * Calculating u(1,1) = -omega * [ AS(1,1) * u(1,0) + AW(1,1) * u(0,1) + AE(1,1) * uOld(2,1) + AN(1,1) * uOld(1,2) ] / AP(1,1) + (1 - omega) * uOld(1,1) u(1,1) = (1.17100) [ -(1.00000)(2.00000) - (1.00000)(1.00000) - (1.00000)(1.42811) - (1.00000)(0.84261) ] / (-4.00000) + (-0.17100)(0.87825) = 1.39282 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 2.88142 3.38014 8.46346 11.00000 j = 2 2.00000 0.84261 0.66475 4.53004 10.00000 j = 1 1.00000 1.39282 1.42811 4.80933 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,1) = -omega * [ AS(2,1) * u(2,0) + AW(2,1) * u(1,1) + AE(2,1) * uOld(3,1) + AN(2,1) * uOld(2,2) ] / AP(2,1) + (1 - omega) * uOld(2,1) u(2,1) = (1.17100) [ -(1.00000)(4.00000) - (1.00000)(1.39282) - (1.00000)(4.80933) - (1.00000)(0.66475) ] / (-4.00000) + (-0.17100)(1.42811) = 2.93708 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 2.88142 3.38014 8.46346 11.00000 j = 2 2.00000 0.84261 0.66475 4.53004 10.00000 j = 1 1.00000 1.39282 2.93708 4.80933 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,1) = -omega * [ AS(3,1) * u(3,0) + AW(3,1) * u(2,1) + AE(3,1) * uOld(4,1) + AN(3,1) * uOld(3,2) ] / AP(3,1) + (1 - omega) * uOld(3,1) u(3,1) = (1.17100) [ -(1.00000)(6.00000) - (1.00000)(2.93708) - (1.00000)(9.00000) - (1.00000)(4.53004) ] / (-4.00000) + (-0.17100)(4.80933) = 5.75485 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 2.88142 3.38014 8.46346 11.00000 j = 2 2.00000 0.84261 0.66475 4.53004 10.00000 j = 1 1.00000 1.39282 2.93708 5.75485 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,2) = -omega * [ AS(1,2) * u(1,1) + AW(1,2) * u(0,2) + AE(1,2) * uOld(2,2) + AN(1,2) * uOld(1,3) ] / AP(1,2) + (1 - omega) * uOld(1,2) u(1,2) = (1.17100) [ -(1.00000)(1.39282) - (1.00000)(2.00000) - (1.00000)(0.66475) - (1.00000)(2.88142) ] / (-4.00000) + (-0.17100)(0.84261) = 1.88731 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 2.88142 3.38014 8.46346 11.00000 j = 2 2.00000 1.88731 0.66475 4.53004 10.00000 j = 1 1.00000 1.39282 2.93708 5.75485 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,2) = -omega * [ AS(2,2) * u(2,1) + AW(2,2) * u(1,2) + AE(2,2) * uOld(3,2) + AN(2,2) * uOld(2,3) ] / AP(2,2) + (1 - omega) * uOld(2,2) u(2,2) = (1.17100) [ -(1.00000)(2.93708) - (1.00000)(1.88731) - (1.00000)(4.53004) - (1.00000)(3.38014) ] / (-4.00000) + (-0.17100)(0.66475) = 3.61437 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 2.88142 3.38014 8.46346 11.00000 j = 2 2.00000 1.88731 3.61437 4.53004 10.00000 j = 1 1.00000 1.39282 2.93708 5.75485 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,2) = -omega * [ AS(3,2) * u(3,1) + AW(3,2) * u(2,2) + AE(3,2) * uOld(4,2) + AN(3,2) * uOld(3,3) ] / AP(3,2) + (1 - omega) * uOld(3,2) u(3,2) = (1.17100) [ -(1.00000)(5.75485) - (1.00000)(3.61437) - (1.00000)(10.00000) - (1.00000)(8.46346) ] / (-4.00000) + (-0.17100)(4.53004) = 7.37338 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 2.88142 3.38014 8.46346 11.00000 j = 2 2.00000 1.88731 3.61437 7.37338 10.00000 j = 1 1.00000 1.39282 2.93708 5.75485 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,3) = -omega * [ AS(1,3) * u(1,2) + AW(1,3) * u(0,3) + AE(1,3) * uOld(2,3) + AN(1,3) * uOld(1,4) ] / AP(1,3) + (1 - omega) * uOld(1,3) u(1,3) = (1.17100) [ -(1.00000)(1.88731) - (1.00000)(3.00000) - (1.00000)(3.38014) - (1.00000)(6.00000) ] / (-4.00000) + (-0.17100)(2.88142) = 3.68407 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 3.68407 3.38014 8.46346 11.00000 j = 2 2.00000 1.88731 3.61437 7.37338 10.00000 j = 1 1.00000 1.39282 2.93708 5.75485 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,3) = -omega * [ AS(2,3) * u(2,2) + AW(2,3) * u(1,3) + AE(2,3) * uOld(3,3) + AN(2,3) * uOld(2,4) ] / AP(2,3) + (1 - omega) * uOld(2,3) u(2,3) = (1.17100) [ -(1.00000)(3.61437) - (1.00000)(3.68407) - (1.00000)(8.46346) - (1.00000)(8.00000) ] / (-4.00000) + (-0.17100)(3.38014) = 6.37829 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 3.68407 6.37829 8.46346 11.00000 j = 2 2.00000 1.88731 3.61437 7.37338 10.00000 j = 1 1.00000 1.39282 2.93708 5.75485 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,3) = -omega * [ AS(3,3) * u(3,2) + AW(3,3) * u(2,3) + AE(3,3) * uOld(4,3) + AN(3,3) * uOld(3,4) ] / AP(3,3) + (1 - omega) * uOld(3,3) u(3,3) = (1.17100) [ -(1.00000)(7.37338) - (1.00000)(6.37829) - (1.00000)(11.00000) - (1.00000)(10.00000) ] / (-4.00000) + (-0.17100)(8.46346) = 8.72630 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 3.68407 6.37829 8.72630 11.00000 j = 2 2.00000 1.88731 3.61437 7.37338 10.00000 j = 1 1.00000 1.39282 2.93708 5.75485 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 At end of 2 iterations, the maximum residual is 1.02964 and the maximum relative increment is 0.81608 Residual array i = 1 i = 2 i = 3 j = 3 0.63233 0.37211 0.03838 j = 2 0.78551 1.02964 0.34950 j = 1 0.56327 0.75343 0.57276 Relative increment array i = 1 i = 2 i = 3 j = 3 0.21787 0.47006 0.03012 j = 2 0.55354 0.81608 0.38562 j = 1 0.36945 0.51377 0.16430 * * * Calculations and results for iteration 3 using successive over relaxation * * * Calculating u(1,1) = -omega * [ AS(1,1) * u(1,0) + AW(1,1) * u(0,1) + AE(1,1) * uOld(2,1) + AN(1,1) * uOld(1,2) ] / AP(1,1) + (1 - omega) * uOld(1,1) u(1,1) = (1.17100) [ -(1.00000)(2.00000) - (1.00000)(1.00000) - (1.00000)(2.93708) - (1.00000)(1.88731) ] / (-4.00000) + (-0.17100)(1.39282) = 2.05242 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 3.68407 6.37829 8.72630 11.00000 j = 2 2.00000 1.88731 3.61437 7.37338 10.00000 j = 1 1.00000 2.05242 2.93708 5.75485 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,1) = -omega * [ AS(2,1) * u(2,0) + AW(2,1) * u(1,1) + AE(2,1) * uOld(3,1) + AN(2,1) * uOld(2,2) ] / AP(2,1) + (1 - omega) * uOld(2,1) u(2,1) = (1.17100) [ -(1.00000)(4.00000) - (1.00000)(2.05242) - (1.00000)(5.75485) - (1.00000)(3.61437) ] / (-4.00000) + (-0.17100)(2.93708) = 4.01244 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 3.68407 6.37829 8.72630 11.00000 j = 2 2.00000 1.88731 3.61437 7.37338 10.00000 j = 1 1.00000 2.05242 4.01244 5.75485 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,1) = -omega * [ AS(3,1) * u(3,0) + AW(3,1) * u(2,1) + AE(3,1) * uOld(4,1) + AN(3,1) * uOld(3,2) ] / AP(3,1) + (1 - omega) * uOld(3,1) u(3,1) = (1.17100) [ -(1.00000)(6.00000) - (1.00000)(4.01244) - (1.00000)(9.00000) - (1.00000)(7.37338) ] / (-4.00000) + (-0.17100)(5.75485) = 6.74037 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 3.68407 6.37829 8.72630 11.00000 j = 2 2.00000 1.88731 3.61437 7.37338 10.00000 j = 1 1.00000 2.05242 4.01244 6.74037 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,2) = -omega * [ AS(1,2) * u(1,1) + AW(1,2) * u(0,2) + AE(1,2) * uOld(2,2) + AN(1,2) * uOld(1,3) ] / AP(1,2) + (1 - omega) * uOld(1,2) u(1,2) = (1.17100) [ -(1.00000)(2.05242) - (1.00000)(2.00000) - (1.00000)(3.61437) - (1.00000)(3.68407) ] / (-4.00000) + (-0.17100)(1.88731) = 3.00023 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 3.68407 6.37829 8.72630 11.00000 j = 2 2.00000 3.00023 3.61437 7.37338 10.00000 j = 1 1.00000 2.05242 4.01244 6.74037 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,2) = -omega * [ AS(2,2) * u(2,1) + AW(2,2) * u(1,2) + AE(2,2) * uOld(3,2) + AN(2,2) * uOld(2,3) ] / AP(2,2) + (1 - omega) * uOld(2,2) u(2,2) = (1.17100) [ -(1.00000)(4.01244) - (1.00000)(3.00023) - (1.00000)(7.37338) - (1.00000)(6.37829) ] / (-4.00000) + (-0.17100)(3.61437) = 5.46071 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 3.68407 6.37829 8.72630 11.00000 j = 2 2.00000 3.00023 5.46071 7.37338 10.00000 j = 1 1.00000 2.05242 4.01244 6.74037 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,2) = -omega * [ AS(3,2) * u(3,1) + AW(3,2) * u(2,2) + AE(3,2) * uOld(4,2) + AN(3,2) * uOld(3,3) ] / AP(3,2) + (1 - omega) * uOld(3,2) u(3,2) = (1.17100) [ -(1.00000)(6.74037) - (1.00000)(5.46071) - (1.00000)(10.00000) - (1.00000)(8.72630) ] / (-4.00000) + (-0.17100)(7.37338) = 7.79314 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 3.68407 6.37829 8.72630 11.00000 j = 2 2.00000 3.00023 5.46071 7.79314 10.00000 j = 1 1.00000 2.05242 4.01244 6.74037 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,3) = -omega * [ AS(1,3) * u(1,2) + AW(1,3) * u(0,3) + AE(1,3) * uOld(2,3) + AN(1,3) * uOld(1,4) ] / AP(1,3) + (1 - omega) * uOld(1,3) u(1,3) = (1.17100) [ -(1.00000)(3.00023) - (1.00000)(3.00000) - (1.00000)(6.37829) - (1.00000)(6.00000) ] / (-4.00000) + (-0.17100)(3.68407) = 4.75034 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.75034 6.37829 8.72630 11.00000 j = 2 2.00000 3.00023 5.46071 7.79314 10.00000 j = 1 1.00000 2.05242 4.01244 6.74037 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,3) = -omega * [ AS(2,3) * u(2,2) + AW(2,3) * u(1,3) + AE(2,3) * uOld(3,3) + AN(2,3) * uOld(2,4) ] / AP(2,3) + (1 - omega) * uOld(2,3) u(2,3) = (1.17100) [ -(1.00000)(5.46071) - (1.00000)(4.75034) - (1.00000)(8.72630) - (1.00000)(8.00000) ] / (-4.00000) + (-0.17100)(6.37829) = 6.79522 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.75034 6.79522 8.72630 11.00000 j = 2 2.00000 3.00023 5.46071 7.79314 10.00000 j = 1 1.00000 2.05242 4.01244 6.74037 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,3) = -omega * [ AS(3,3) * u(3,2) + AW(3,3) * u(2,3) + AE(3,3) * uOld(4,3) + AN(3,3) * uOld(3,4) ] / AP(3,3) + (1 - omega) * uOld(3,3) u(3,3) = (1.17100) [ -(1.00000)(7.79314) - (1.00000)(6.79522) - (1.00000)(11.00000) - (1.00000)(10.00000) ] / (-4.00000) + (-0.17100)(8.72630) = 8.92630 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.75034 6.79522 8.92630 11.00000 j = 2 2.00000 3.00023 5.46071 7.79314 10.00000 j = 1 1.00000 2.05242 4.01244 6.74037 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 At end of 3 iterations, the maximum residual is 0.56563 and the maximum relative increment is 0.37095 Residual array i = 1 i = 2 i = 3 j = 3 0.05147 0.01089 0.02920 j = 2 0.56563 0.06045 0.01130 j = 1 0.45075 0.55093 0.03897 Relative increment array i = 1 i = 2 i = 3 j = 3 0.22446 0.06136 0.02241 j = 2 0.37095 0.33811 0.05386 j = 1 0.32137 0.26801 0.14621 * * Calculations and results for iteration 4 using successive over relaxation * * Calculating u(1,1) = -omega * [ AS(1,1) * u(1,0) + AW(1,1) * u(0,1) + AE(1,1) * uOld(2,1) + AN(1,1) * uOld(1,2) ] / AP(1,1) + (1 - omega) * uOld(1,1) u(1,1) = (1.17100) [ -(1.00000)(2.00000) - (1.00000)(1.00000) - (1.00000)(4.01244) - (1.00000)(3.00023) ] / (-4.00000) + (-0.17100)(2.05242) = 2.58025 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.75034 6.79522 8.92630 11.00000 j = 2 2.00000 3.00023 5.46071 7.79314 10.00000 j = 1 1.00000 2.58025 4.01244 6.74037 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,1) = -omega * [ AS(2,1) * u(2,0) + AW(2,1) * u(1,1) + AE(2,1) * uOld(3,1) + AN(2,1) * uOld(2,2) ] / AP(2,1) + (1 - omega) * uOld(2,1) u(2,1) = (1.17100) [ -(1.00000)(4.00000) - (1.00000)(2.58025) - (1.00000)(6.74037) - (1.00000)(5.46071) ] / (-4.00000) + (-0.17100)(4.01244) = 4.81211 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.75034 6.79522 8.92630 11.00000 j = 2 2.00000 3.00023 5.46071 7.79314 10.00000 j = 1 1.00000 2.58025 4.81211 6.74037 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,1) = -omega * [ AS(3,1) * u(3,0) + AW(3,1) * u(2,1) + AE(3,1) * uOld(4,1) + AN(3,1) * uOld(3,2) ] / AP(3,1) + (1 - omega) * uOld(3,1) u(3,1) = (1.17100) [ -(1.00000)(6.00000) - (1.00000)(4.81211) - (1.00000)(9.00000) - (1.00000)(7.79314) ] / (-4.00000) + (-0.17100)(6.74037) = 6.92883 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.75034 6.79522 8.92630 11.00000 j = 2 2.00000 3.00023 5.46071 7.79314 10.00000 j = 1 1.00000 2.58025 4.81211 6.92883 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,2) = -omega * [ AS(1,2) * u(1,1) + AW(1,2) * u(0,2) + AE(1,2) * uOld(2,2) + AN(1,2) * uOld(1,3) ] / AP(1,2) + (1 - omega) * uOld(1,2) u(1,2) = (1.17100) [ -(1.00000)(2.58025) - (1.00000)(2.00000) - (1.00000)(5.46071) - (1.00000)(4.75034) ] / (-4.00000) + (-0.17100)(3.00023) = 3.81711 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.75034 6.79522 8.92630 11.00000 j = 2 2.00000 3.81711 5.46071 7.79314 10.00000 j = 1 1.00000 2.58025 4.81211 6.92883 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,2) = -omega * [ AS(2,2) * u(2,1) + AW(2,2) * u(1,2) + AE(2,2) * uOld(3,2) + AN(2,2) * uOld(2,3) ] / AP(2,2) + (1 - omega) * uOld(2,2) u(2,2) = (1.17100) [ -(1.00000)(4.81211) - (1.00000)(3.81711) - (1.00000)(7.79314) - (1.00000)(6.79522) ] / (-4.00000) + (-0.17100)(5.46071) = 5.86316 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.75034 6.79522 8.92630 11.00000 j = 2 2.00000 3.81711 5.86316 7.79314 10.00000 j = 1 1.00000 2.58025 4.81211 6.92883 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,2) = -omega * [ AS(3,2) * u(3,1) + AW(3,2) * u(2,2) + AE(3,2) * uOld(4,2) + AN(3,2) * uOld(3,3) ] / AP(3,2) + (1 - omega) * uOld(3,2) u(3,2) = (1.17100) [ -(1.00000)(6.92883) - (1.00000)(5.86316) - (1.00000)(10.00000) - (1.00000)(8.92630) ] / (-4.00000) + (-0.17100)(7.79314) = 7.95290 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.75034 6.79522 8.92630 11.00000 j = 2 2.00000 3.81711 5.86316 7.95290 10.00000 j = 1 1.00000 2.58025 4.81211 6.92883 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,3) = -omega * [ AS(1,3) * u(1,2) + AW(1,3) * u(0,3) + AE(1,3) * uOld(2,3) + AN(1,3) * uOld(1,4) ] / AP(1,3) + (1 - omega) * uOld(1,3) u(1,3) = (1.17100) [ -(1.00000)(3.81711) - (1.00000)(3.00000) - (1.00000)(6.79522) - (1.00000)(6.00000) ] / (-4.00000) + (-0.17100)(4.75034) = 4.92920 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.92920 6.79522 8.92630 11.00000 j = 2 2.00000 3.81711 5.86316 7.95290 10.00000 j = 1 1.00000 2.58025 4.81211 6.92883 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,3) = -omega * [ AS(2,3) * u(2,2) + AW(2,3) * u(1,3) + AE(2,3) * uOld(3,3) + AN(2,3) * uOld(2,4) ] / AP(2,3) + (1 - omega) * uOld(2,3) u(2,3) = (1.17100) [ -(1.00000)(5.86316) - (1.00000)(4.92920) - (1.00000)(8.92630) - (1.00000)(8.00000) ] / (-4.00000) + (-0.17100)(6.79522) = 6.95266 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.92920 6.95266 8.92630 11.00000 j = 2 2.00000 3.81711 5.86316 7.95290 10.00000 j = 1 1.00000 2.58025 4.81211 6.92883 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,3) = -omega * [ AS(3,3) * u(3,2) + AW(3,3) * u(2,3) + AE(3,3) * uOld(4,3) + AN(3,3) * uOld(3,4) ] / AP(3,3) + (1 - omega) * uOld(3,3) u(3,3) = (1.17100) [ -(1.00000)(7.95290) - (1.00000)(6.95266) - (1.00000)(11.00000) - (1.00000)(10.00000) ] / (-4.00000) + (-0.17100)(8.92630) = 8.98496 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.92920 6.95266 8.98496 11.00000 j = 2 2.00000 3.81711 5.86316 7.95290 10.00000 j = 1 1.00000 2.58025 4.81211 6.92883 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 At end of 4 iterations, the maximum residual is 0.32706 and the maximum relative increment is 0.21400 Residual array i = 1 i = 2 i = 3 j = 3 0.01324 0.00833 0.00857 j = 2 0.02604 0.02053 0.00866 j = 1 0.32706 0.03096 0.01242 Relative increment array i = 1 i = 2 i = 3 j = 3 0.03629 0.02264 0.00653 j = 2 0.21400 0.06864 0.02009 j = 1 0.20457 0.16618 0.02720 * * * Calculations and results for iteration 5 using successive over relaxation * * * Calculating u(1,1) = -omega * [ AS(1,1) * u(1,0) + AW(1,1) * u(0,1) + AE(1,1) * uOld(2,1) + AN(1,1) * uOld(1,2) ] / AP(1,1) + (1 - omega) * uOld(1,1) u(1,1) = (1.17100) [ -(1.00000)(2.00000) - (1.00000)(1.00000) - (1.00000)(4.81211) - (1.00000)(3.81711) ] / (-4.00000) + (-0.17100)(2.58025) = 2.96323 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.92920 6.95266 8.98496 11.00000 j = 2 2.00000 3.81711 5.86316 7.95290 10.00000 j = 1 1.00000 2.96323 4.81211 6.92883 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,1) = -omega * [ AS(2,1) * u(2,0) + AW(2,1) * u(1,1) + AE(2,1) * uOld(3,1) + AN(2,1) * uOld(2,2) ] / AP(2,1) + (1 - omega) * uOld(2,1) u(2,1) = (1.17100) [ -(1.00000)(4.00000) - (1.00000)(2.96323) - (1.00000)(6.92883) - (1.00000)(5.86316) ] / (-4.00000) + (-0.17100)(4.81211) = 4.96047 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.92920 6.95266 8.98496 11.00000 j = 2 2.00000 3.81711 5.86316 7.95290 10.00000 j = 1 1.00000 2.96323 4.96047 6.92883 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,1) = -omega * [ AS(3,1) * u(3,0) + AW(3,1) * u(2,1) + AE(3,1) * uOld(4,1) + AN(3,1) * uOld(3,2) ] / AP(3,1) + (1 - omega) * uOld(3,1) u(3,1) = (1.17100) [ -(1.00000)(6.00000) - (1.00000)(4.96047) - (1.00000)(9.00000) - (1.00000)(7.95290) ] / (-4.00000) + (-0.17100)(6.92883) = 6.98681 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.92920 6.95266 8.98496 11.00000 j = 2 2.00000 3.81711 5.86316 7.95290 10.00000 j = 1 1.00000 2.96323 4.96047 6.98681 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,2) = -omega * [ AS(1,2) * u(1,1) + AW(1,2) * u(0,2) + AE(1,2) * uOld(2,2) + AN(1,2) * uOld(1,3) ] / AP(1,2) + (1 - omega) * uOld(1,2) u(1,2) = (1.17100) [ -(1.00000)(2.96323) - (1.00000)(2.00000) - (1.00000)(5.86316) - (1.00000)(4.92920) ] / (-4.00000) + (-0.17100)(3.81711) = 3.95973 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.92920 6.95266 8.98496 11.00000 j = 2 2.00000 3.95973 5.86316 7.95290 10.00000 j = 1 1.00000 2.96323 4.96047 6.98681 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,2) = -omega * [ AS(2,2) * u(2,1) + AW(2,2) * u(1,2) + AE(2,2) * uOld(3,2) + AN(2,2) * uOld(2,3) ] / AP(2,2) + (1 - omega) * uOld(2,2) u(2,2) = (1.17100) [ -(1.00000)(4.96047) - (1.00000)(3.95973) - (1.00000)(7.95290) - (1.00000)(6.95266) ] / (-4.00000) + (-0.17100)(5.86316) = 5.97239 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.92920 6.95266 8.98496 11.00000 j = 2 2.00000 3.95973 5.97239 7.95290 10.00000 j = 1 1.00000 2.96323 4.96047 6.98681 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,2) = -omega * [ AS(3,2) * u(3,1) + AW(3,2) * u(2,2) + AE(3,2) * uOld(4,2) + AN(3,2) * uOld(3,3) ] / AP(3,2) + (1 - omega) * uOld(3,2) u(3,2) = (1.17100) [ -(1.00000)(6.98681) - (1.00000)(5.97239) - (1.00000)(10.00000) - (1.00000)(8.98496) ] / (-4.00000) + (-0.17100)(7.95290) = 7.99171 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.92920 6.95266 8.98496 11.00000 j = 2 2.00000 3.95973 5.97239 7.99171 10.00000 j = 1 1.00000 2.96323 4.96047 6.98681 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,3) = -omega * [ AS(1,3) * u(1,2) + AW(1,3) * u(0,3) + AE(1,3) * uOld(2,3) + AN(1,3) * uOld(1,4) ] / AP(1,3) + (1 - omega) * uOld(1,3) u(1,3) = (1.17100) [ -(1.00000)(3.95973) - (1.00000)(3.00000) - (1.00000)(6.95266) - (1.00000)(6.00000) ] / (-4.00000) + (-0.17100)(4.92920) = 4.98646 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.98646 6.95266 8.98496 11.00000 j = 2 2.00000 3.95973 5.97239 7.99171 10.00000 j = 1 1.00000 2.96323 4.96047 6.98681 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,3) = -omega * [ AS(2,3) * u(2,2) + AW(2,3) * u(1,3) + AE(2,3) * uOld(3,3) + AN(2,3) * uOld(2,4) ] / AP(2,3) + (1 - omega) * uOld(2,3) u(2,3) = (1.17100) [ -(1.00000)(5.97239) - (1.00000)(4.98646) - (1.00000)(8.98496) - (1.00000)(8.00000) ] / (-4.00000) + (-0.17100)(6.95266) = 6.99164 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.98646 6.99164 8.98496 11.00000 j = 2 2.00000 3.95973 5.97239 7.99171 10.00000 j = 1 1.00000 2.96323 4.96047 6.98681 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,3) = -omega * [ AS(3,3) * u(3,2) + AW(3,3) * u(2,3) + AE(3,3) * uOld(4,3) + AN(3,3) * uOld(3,4) ] / AP(3,3) + (1 - omega) * uOld(3,3) u(3,3) = (1.17100) [ -(1.00000)(7.99171) - (1.00000)(6.99164) - (1.00000)(11.00000) - (1.00000)(10.00000) ] / (-4.00000) + (-0.17100)(8.98496) = 8.99770 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.98646 6.99164 8.99770 11.00000 j = 2 2.00000 3.95973 5.97239 7.99171 10.00000 j = 1 1.00000 2.96323 4.96047 6.98681 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 At end of 5 iterations, the maximum residual is 0.02079 and the maximum relative increment is 0.12924 Residual array i = 1 i = 2 i = 3 j = 3 0.00139 0.00251 0.00186 j = 2 0.02079 0.00350 0.00248 j = 1 0.01682 0.02013 0.00123 Relative increment array i = 1 i = 2 i = 3 j = 3 0.01148 0.00558 0.00142 j = 2 0.03602 0.01829 0.00486 j = 1 0.12924 0.02991 0.00830 * * * Calculations and results for iteration 6 using successive over relaxation * * * Calculating u(1,1) = -omega * [ AS(1,1) * u(1,0) + AW(1,1) * u(0,1) + AE(1,1) * uOld(2,1) + AN(1,1) * uOld(1,2) ] / AP(1,1) + (1 - omega) * uOld(1,1) u(1,1) = (1.17100) [ -(1.00000)(2.00000) - (1.00000)(1.00000) - (1.00000)(4.96047) - (1.00000)(3.95973) ] / (-4.00000) + (-0.17100)(2.96323) = 2.98293 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.98646 6.99164 8.99770 11.00000 j = 2 2.00000 3.95973 5.97239 7.99171 10.00000 j = 1 1.00000 2.98293 4.96047 6.98681 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,1) = -omega * [ AS(2,1) * u(2,0) + AW(2,1) * u(1,1) + AE(2,1) * uOld(3,1) + AN(2,1) * uOld(2,2) ] / AP(2,1) + (1 - omega) * uOld(2,1) u(2,1) = (1.17100) [ -(1.00000)(4.00000) - (1.00000)(2.98293) - (1.00000)(6.98681) - (1.00000)(5.97239) ] / (-4.00000) + (-0.17100)(4.96047) = 4.98982 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.98646 6.99164 8.99770 11.00000 j = 2 2.00000 3.95973 5.97239 7.99171 10.00000 j = 1 1.00000 2.98293 4.98982 6.98681 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,1) = -omega * [ AS(3,1) * u(3,0) + AW(3,1) * u(2,1) + AE(3,1) * uOld(4,1) + AN(3,1) * uOld(3,2) ] / AP(3,1) + (1 - omega) * uOld(3,1) u(3,1) = (1.17100) [ -(1.00000)(6.00000) - (1.00000)(4.98982) - (1.00000)(9.00000) - (1.00000)(7.99171) ] / (-4.00000) + (-0.17100)(6.98681) = 6.99685 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.98646 6.99164 8.99770 11.00000 j = 2 2.00000 3.95973 5.97239 7.99171 10.00000 j = 1 1.00000 2.98293 4.98982 6.99685 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,2) = -omega * [ AS(1,2) * u(1,1) + AW(1,2) * u(0,2) + AE(1,2) * uOld(2,2) + AN(1,2) * uOld(1,3) ] / AP(1,2) + (1 - omega) * uOld(1,2) u(1,2) = (1.17100) [ -(1.00000)(2.98293) - (1.00000)(2.00000) - (1.00000)(5.97239) - (1.00000)(4.98646) ] / (-4.00000) + (-0.17100)(3.95973) = 3.98984 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.98646 6.99164 8.99770 11.00000 j = 2 2.00000 3.98984 5.97239 7.99171 10.00000 j = 1 1.00000 2.98293 4.98982 6.99685 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,2) = -omega * [ AS(2,2) * u(2,1) + AW(2,2) * u(1,2) + AE(2,2) * uOld(3,2) + AN(2,2) * uOld(2,3) ] / AP(2,2) + (1 - omega) * uOld(2,2) u(2,2) = (1.17100) [ -(1.00000)(4.98982) - (1.00000)(3.98984) - (1.00000)(7.99171) - (1.00000)(6.99164) ] / (-4.00000) + (-0.17100)(5.97239) = 5.99389 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.98646 6.99164 8.99770 11.00000 j = 2 2.00000 3.98984 5.99389 7.99171 10.00000 j = 1 1.00000 2.98293 4.98982 6.99685 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,2) = -omega * [ AS(3,2) * u(3,1) + AW(3,2) * u(2,2) + AE(3,2) * uOld(4,2) + AN(3,2) * uOld(3,3) ] / AP(3,2) + (1 - omega) * uOld(3,2) u(3,2) = (1.17100) [ -(1.00000)(6.99685) - (1.00000)(5.99389) - (1.00000)(10.00000) - (1.00000)(8.99770) ] / (-4.00000) + (-0.17100)(7.99171) = 7.99803 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.98646 6.99164 8.99770 11.00000 j = 2 2.00000 3.98984 5.99389 7.99803 10.00000 j = 1 1.00000 2.98293 4.98982 6.99685 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,3) = -omega * [ AS(1,3) * u(1,2) + AW(1,3) * u(0,3) + AE(1,3) * uOld(2,3) + AN(1,3) * uOld(1,4) ] / AP(1,3) + (1 - omega) * uOld(1,3) u(1,3) = (1.17100) [ -(1.00000)(3.98984) - (1.00000)(3.00000) - (1.00000)(6.99164) - (1.00000)(6.00000) ] / (-4.00000) + (-0.17100)(4.98646) = 4.99690 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99690 6.99164 8.99770 11.00000 j = 2 2.00000 3.98984 5.99389 7.99803 10.00000 j = 1 1.00000 2.98293 4.98982 6.99685 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,3) = -omega * [ AS(2,3) * u(2,2) + AW(2,3) * u(1,3) + AE(2,3) * uOld(3,3) + AN(2,3) * uOld(2,4) ] / AP(2,3) + (1 - omega) * uOld(2,3) u(2,3) = (1.17100) [ -(1.00000)(5.99389) - (1.00000)(4.99690) - (1.00000)(8.99770) - (1.00000)(8.00000) ] / (-4.00000) + (-0.17100)(6.99164) = 6.99806 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99690 6.99806 8.99770 11.00000 j = 2 2.00000 3.98984 5.99389 7.99803 10.00000 j = 1 1.00000 2.98293 4.98982 6.99685 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,3) = -omega * [ AS(3,3) * u(3,2) + AW(3,3) * u(2,3) + AE(3,3) * uOld(4,3) + AN(3,3) * uOld(3,4) ] / AP(3,3) + (1 - omega) * uOld(3,3) u(3,3) = (1.17100) [ -(1.00000)(7.99803) - (1.00000)(6.99806) - (1.00000)(11.00000) - (1.00000)(10.00000) ] / (-4.00000) + (-0.17100)(8.99770) = 8.99925 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99690 6.99806 8.99925 11.00000 j = 2 2.00000 3.98984 5.99389 7.99803 10.00000 j = 1 1.00000 2.98293 4.98982 6.99685 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 At end of 6 iterations, the maximum residual is 0.01199 and the maximum relative increment is 0.00755 Residual array i = 1 i = 2 i = 3 j = 3 0.00008 0.00055 0.00023 j = 2 0.00359 0.00005 0.00054 j = 1 0.01199 0.00360 0.00012 Relative increment array i = 1 i = 2 i = 3 j = 3 0.00209 0.00092 0.00017 j = 2 0.00755 0.00359 0.00079 j = 1 0.00660 0.00588 0.00143 * * * Calculations and results for iteration 7 using successive over relaxation * * * Calculating u(1,1) = -omega * [ AS(1,1) * u(1,0) + AW(1,1) * u(0,1) + AE(1,1) * uOld(2,1) + AN(1,1) * uOld(1,2) ] / AP(1,1) + (1 - omega) * uOld(1,1) u(1,1) = (1.17100) [ -(1.00000)(2.00000) - (1.00000)(1.00000) - (1.00000)(4.98982) - (1.00000)(3.98984) ] / (-4.00000) + (-0.17100)(2.98293) = 2.99696 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99690 6.99806 8.99925 11.00000 j = 2 2.00000 3.98984 5.99389 7.99803 10.00000 j = 1 1.00000 2.99696 4.98982 6.99685 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,1) = -omega * [ AS(2,1) * u(2,0) + AW(2,1) * u(1,1) + AE(2,1) * uOld(3,1) + AN(2,1) * uOld(2,2) ] / AP(2,1) + (1 - omega) * uOld(2,1) u(2,1) = (1.17100) [ -(1.00000)(4.00000) - (1.00000)(2.99696) - (1.00000)(6.99685) - (1.00000)(5.99389) ] / (-4.00000) + (-0.17100)(4.98982) = 4.99814 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99690 6.99806 8.99925 11.00000 j = 2 2.00000 3.98984 5.99389 7.99803 10.00000 j = 1 1.00000 2.99696 4.99814 6.99685 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,1) = -omega * [ AS(3,1) * u(3,0) + AW(3,1) * u(2,1) + AE(3,1) * uOld(4,1) + AN(3,1) * uOld(3,2) ] / AP(3,1) + (1 - omega) * uOld(3,1) u(3,1) = (1.17100) [ -(1.00000)(6.00000) - (1.00000)(4.99814) - (1.00000)(9.00000) - (1.00000)(7.99803) ] / (-4.00000) + (-0.17100)(6.99685) = 6.99942 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99690 6.99806 8.99925 11.00000 j = 2 2.00000 3.98984 5.99389 7.99803 10.00000 j = 1 1.00000 2.99696 4.99814 6.99942 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,2) = -omega * [ AS(1,2) * u(1,1) + AW(1,2) * u(0,2) + AE(1,2) * uOld(2,2) + AN(1,2) * uOld(1,3) ] / AP(1,2) + (1 - omega) * uOld(1,2) u(1,2) = (1.17100) [ -(1.00000)(2.99696) - (1.00000)(2.00000) - (1.00000)(5.99389) - (1.00000)(4.99690) ] / (-4.00000) + (-0.17100)(3.98984) = 3.99815 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99690 6.99806 8.99925 11.00000 j = 2 2.00000 3.99815 5.99389 7.99803 10.00000 j = 1 1.00000 2.99696 4.99814 6.99942 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,2) = -omega * [ AS(2,2) * u(2,1) + AW(2,2) * u(1,2) + AE(2,2) * uOld(3,2) + AN(2,2) * uOld(2,3) ] / AP(2,2) + (1 - omega) * uOld(2,2) u(2,2) = (1.17100) [ -(1.00000)(4.99814) - (1.00000)(3.99815) - (1.00000)(7.99803) - (1.00000)(6.99806) ] / (-4.00000) + (-0.17100)(5.99389) = 5.99881 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99690 6.99806 8.99925 11.00000 j = 2 2.00000 3.99815 5.99881 7.99803 10.00000 j = 1 1.00000 2.99696 4.99814 6.99942 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,2) = -omega * [ AS(3,2) * u(3,1) + AW(3,2) * u(2,2) + AE(3,2) * uOld(4,2) + AN(3,2) * uOld(3,3) ] / AP(3,2) + (1 - omega) * uOld(3,2) u(3,2) = (1.17100) [ -(1.00000)(6.99942) - (1.00000)(5.99881) - (1.00000)(10.00000) - (1.00000)(8.99925) ] / (-4.00000) + (-0.17100)(7.99803) = 7.99960 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99690 6.99806 8.99925 11.00000 j = 2 2.00000 3.99815 5.99881 7.99960 10.00000 j = 1 1.00000 2.99696 4.99814 6.99942 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,3) = -omega * [ AS(1,3) * u(1,2) + AW(1,3) * u(0,3) + AE(1,3) * uOld(2,3) + AN(1,3) * uOld(1,4) ] / AP(1,3) + (1 - omega) * uOld(1,3) u(1,3) = (1.17100) [ -(1.00000)(3.99815) - (1.00000)(3.00000) - (1.00000)(6.99806) - (1.00000)(6.00000) ] / (-4.00000) + (-0.17100)(4.99690) = 4.99942 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99942 6.99806 8.99925 11.00000 j = 2 2.00000 3.99815 5.99881 7.99960 10.00000 j = 1 1.00000 2.99696 4.99814 6.99942 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,3) = -omega * [ AS(2,3) * u(2,2) + AW(2,3) * u(1,3) + AE(2,3) * uOld(3,3) + AN(2,3) * uOld(2,4) ] / AP(2,3) + (1 - omega) * uOld(2,3) u(2,3) = (1.17100) [ -(1.00000)(5.99881) - (1.00000)(4.99942) - (1.00000)(8.99925) - (1.00000)(8.00000) ] / (-4.00000) + (-0.17100)(6.99806) = 6.99960 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99942 6.99960 8.99925 11.00000 j = 2 2.00000 3.99815 5.99881 7.99960 10.00000 j = 1 1.00000 2.99696 4.99814 6.99942 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,3) = -omega * [ AS(3,3) * u(3,2) + AW(3,3) * u(2,3) + AE(3,3) * uOld(4,3) + AN(3,3) * uOld(3,4) ] / AP(3,3) + (1 - omega) * uOld(3,3) u(3,3) = (1.17100) [ -(1.00000)(7.99960) - (1.00000)(6.99960) - (1.00000)(11.00000) - (1.00000)(10.00000) ] / (-4.00000) + (-0.17100)(8.99925) = 8.99989 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99942 6.99960 8.99989 11.00000 j = 2 2.00000 3.99815 5.99881 7.99960 10.00000 j = 1 1.00000 2.99696 4.99814 6.99942 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 At end of 7 iterations, the maximum residual is 0.00211 and the maximum relative increment is 0.00468 Residual array i = 1 i = 2 i = 3 j = 3 0.00002 0.00006 0.00009 j = 2 0.00065 0.00006 0.00007 j = 1 0.00211 0.00066 0.00002 Relative increment array i = 1 i = 2 i = 3 j = 3 0.00051 0.00022 0.00007 j = 2 0.00208 0.00082 0.00020 j = 1 0.00468 0.00167 0.00037 * * * Calculations and results for iteration 8 using successive over relaxation * * * Calculating u(1,1) = -omega * [ AS(1,1) * u(1,0) + AW(1,1) * u(0,1) + AE(1,1) * uOld(2,1) + AN(1,1) * uOld(1,2) ] / AP(1,1) + (1 - omega) * uOld(1,1) u(1,1) = (1.17100) [ -(1.00000)(2.00000) - (1.00000)(1.00000) - (1.00000)(4.99814) - (1.00000)(3.99815) ] / (-4.00000) + (-0.17100)(2.99696) = 2.99943 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99942 6.99960 8.99989 11.00000 j = 2 2.00000 3.99815 5.99881 7.99960 10.00000 j = 1 1.00000 2.99943 4.99814 6.99942 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,1) = -omega * [ AS(2,1) * u(2,0) + AW(2,1) * u(1,1) + AE(2,1) * uOld(3,1) + AN(2,1) * uOld(2,2) ] / AP(2,1) + (1 - omega) * uOld(2,1) u(2,1) = (1.17100) [ -(1.00000)(4.00000) - (1.00000)(2.99943) - (1.00000)(6.99942) - (1.00000)(5.99881) ] / (-4.00000) + (-0.17100)(4.99814) = 4.99964 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99942 6.99960 8.99989 11.00000 j = 2 2.00000 3.99815 5.99881 7.99960 10.00000 j = 1 1.00000 2.99943 4.99964 6.99942 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,1) = -omega * [ AS(3,1) * u(3,0) + AW(3,1) * u(2,1) + AE(3,1) * uOld(4,1) + AN(3,1) * uOld(3,2) ] / AP(3,1) + (1 - omega) * uOld(3,1) u(3,1) = (1.17100) [ -(1.00000)(6.00000) - (1.00000)(4.99964) - (1.00000)(9.00000) - (1.00000)(7.99960) ] / (-4.00000) + (-0.17100)(6.99942) = 6.99988 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99942 6.99960 8.99989 11.00000 j = 2 2.00000 3.99815 5.99881 7.99960 10.00000 j = 1 1.00000 2.99943 4.99964 6.99988 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,2) = -omega * [ AS(1,2) * u(1,1) + AW(1,2) * u(0,2) + AE(1,2) * uOld(2,2) + AN(1,2) * uOld(1,3) ] / AP(1,2) + (1 - omega) * uOld(1,2) u(1,2) = (1.17100) [ -(1.00000)(2.99943) - (1.00000)(2.00000) - (1.00000)(5.99881) - (1.00000)(4.99942) ] / (-4.00000) + (-0.17100)(3.99815) = 3.99963 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99942 6.99960 8.99989 11.00000 j = 2 2.00000 3.99963 5.99881 7.99960 10.00000 j = 1 1.00000 2.99943 4.99964 6.99988 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,2) = -omega * [ AS(2,2) * u(2,1) + AW(2,2) * u(1,2) + AE(2,2) * uOld(3,2) + AN(2,2) * uOld(2,3) ] / AP(2,2) + (1 - omega) * uOld(2,2) u(2,2) = (1.17100) [ -(1.00000)(4.99964) - (1.00000)(3.99963) - (1.00000)(7.99960) - (1.00000)(6.99960) ] / (-4.00000) + (-0.17100)(5.99881) = 5.99975 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99942 6.99960 8.99989 11.00000 j = 2 2.00000 3.99963 5.99975 7.99960 10.00000 j = 1 1.00000 2.99943 4.99964 6.99988 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,2) = -omega * [ AS(3,2) * u(3,1) + AW(3,2) * u(2,2) + AE(3,2) * uOld(4,2) + AN(3,2) * uOld(3,3) ] / AP(3,2) + (1 - omega) * uOld(3,2) u(3,2) = (1.17100) [ -(1.00000)(6.99988) - (1.00000)(5.99975) - (1.00000)(10.00000) - (1.00000)(8.99989) ] / (-4.00000) + (-0.17100)(7.99960) = 7.99993 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99942 6.99960 8.99989 11.00000 j = 2 2.00000 3.99963 5.99975 7.99993 10.00000 j = 1 1.00000 2.99943 4.99964 6.99988 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,3) = -omega * [ AS(1,3) * u(1,2) + AW(1,3) * u(0,3) + AE(1,3) * uOld(2,3) + AN(1,3) * uOld(1,4) ] / AP(1,3) + (1 - omega) * uOld(1,3) u(1,3) = (1.17100) [ -(1.00000)(3.99963) - (1.00000)(3.00000) - (1.00000)(6.99960) - (1.00000)(6.00000) ] / (-4.00000) + (-0.17100)(4.99942) = 4.99987 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99987 6.99960 8.99989 11.00000 j = 2 2.00000 3.99963 5.99975 7.99993 10.00000 j = 1 1.00000 2.99943 4.99964 6.99988 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,3) = -omega * [ AS(2,3) * u(2,2) + AW(2,3) * u(1,3) + AE(2,3) * uOld(3,3) + AN(2,3) * uOld(2,4) ] / AP(2,3) + (1 - omega) * uOld(2,3) u(2,3) = (1.17100) [ -(1.00000)(5.99975) - (1.00000)(4.99987) - (1.00000)(8.99989) - (1.00000)(8.00000) ] / (-4.00000) + (-0.17100)(6.99960) = 6.99993 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99987 6.99993 8.99989 11.00000 j = 2 2.00000 3.99963 5.99975 7.99993 10.00000 j = 1 1.00000 2.99943 4.99964 6.99988 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,3) = -omega * [ AS(3,3) * u(3,2) + AW(3,3) * u(2,3) + AE(3,3) * uOld(4,3) + AN(3,3) * uOld(3,4) ] / AP(3,3) + (1 - omega) * uOld(3,3) u(3,3) = (1.17100) [ -(1.00000)(7.99993) - (1.00000)(6.99993) - (1.00000)(11.00000) - (1.00000)(10.00000) ] / (-4.00000) + (-0.17100)(8.99989) = 8.99998 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99987 6.99993 8.99998 11.00000 j = 2 2.00000 3.99963 5.99975 7.99993 10.00000 j = 1 1.00000 2.99943 4.99964 6.99988 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 At end of 8 iterations, the maximum residual is 0.00038 and the maximum relative increment is 0.00082 Residual array i = 1 i = 2 i = 3 j = 3 0.00002 0.00003 0.00001 j = 2 0.00013 0.00003 0.00003 j = 1 0.00038 0.00013 0.00002 Relative increment array i = 1 i = 2 i = 3 j = 3 0.00009 0.00005 0.00001 j = 2 0.00037 0.00016 0.00004 j = 1 0.00082 0.00030 0.00007 * * * Calculations and results for iteration 9 using successive over relaxation * * * Calculating u(1,1) = -omega * [ AS(1,1) * u(1,0) + AW(1,1) * u(0,1) + AE(1,1) * uOld(2,1) + AN(1,1) * uOld(1,2) ] / AP(1,1) + (1 - omega) * uOld(1,1) u(1,1) = (1.17100) [ -(1.00000)(2.00000) - (1.00000)(1.00000) - (1.00000)(4.99964) - (1.00000)(3.99963) ] / (-4.00000) + (-0.17100)(2.99943) = 2.99988 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99987 6.99993 8.99998 11.00000 j = 2 2.00000 3.99963 5.99975 7.99993 10.00000 j = 1 1.00000 2.99988 4.99964 6.99988 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,1) = -omega * [ AS(2,1) * u(2,0) + AW(2,1) * u(1,1) + AE(2,1) * uOld(3,1) + AN(2,1) * uOld(2,2) ] / AP(2,1) + (1 - omega) * uOld(2,1) u(2,1) = (1.17100) [ -(1.00000)(4.00000) - (1.00000)(2.99988) - (1.00000)(6.99988) - (1.00000)(5.99975) ] / (-4.00000) + (-0.17100)(4.99964) = 4.99992 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99987 6.99993 8.99998 11.00000 j = 2 2.00000 3.99963 5.99975 7.99993 10.00000 j = 1 1.00000 2.99988 4.99992 6.99988 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,1) = -omega * [ AS(3,1) * u(3,0) + AW(3,1) * u(2,1) + AE(3,1) * uOld(4,1) + AN(3,1) * uOld(3,2) ] / AP(3,1) + (1 - omega) * uOld(3,1) u(3,1) = (1.17100) [ -(1.00000)(6.00000) - (1.00000)(4.99992) - (1.00000)(9.00000) - (1.00000)(7.99993) ] / (-4.00000) + (-0.17100)(6.99988) = 6.99998 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99987 6.99993 8.99998 11.00000 j = 2 2.00000 3.99963 5.99975 7.99993 10.00000 j = 1 1.00000 2.99988 4.99992 6.99998 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,2) = -omega * [ AS(1,2) * u(1,1) + AW(1,2) * u(0,2) + AE(1,2) * uOld(2,2) + AN(1,2) * uOld(1,3) ] / AP(1,2) + (1 - omega) * uOld(1,2) u(1,2) = (1.17100) [ -(1.00000)(2.99988) - (1.00000)(2.00000) - (1.00000)(5.99975) - (1.00000)(4.99987) ] / (-4.00000) + (-0.17100)(3.99963) = 3.99992 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99987 6.99993 8.99998 11.00000 j = 2 2.00000 3.99992 5.99975 7.99993 10.00000 j = 1 1.00000 2.99988 4.99992 6.99998 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,2) = -omega * [ AS(2,2) * u(2,1) + AW(2,2) * u(1,2) + AE(2,2) * uOld(3,2) + AN(2,2) * uOld(2,3) ] / AP(2,2) + (1 - omega) * uOld(2,2) u(2,2) = (1.17100) [ -(1.00000)(4.99992) - (1.00000)(3.99992) - (1.00000)(7.99993) - (1.00000)(6.99993) ] / (-4.00000) + (-0.17100)(5.99975) = 5.99995 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99987 6.99993 8.99998 11.00000 j = 2 2.00000 3.99992 5.99995 7.99993 10.00000 j = 1 1.00000 2.99988 4.99992 6.99998 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,2) = -omega * [ AS(3,2) * u(3,1) + AW(3,2) * u(2,2) + AE(3,2) * uOld(4,2) + AN(3,2) * uOld(3,3) ] / AP(3,2) + (1 - omega) * uOld(3,2) u(3,2) = (1.17100) [ -(1.00000)(6.99998) - (1.00000)(5.99995) - (1.00000)(10.00000) - (1.00000)(8.99998) ] / (-4.00000) + (-0.17100)(7.99993) = 7.99998 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99987 6.99993 8.99998 11.00000 j = 2 2.00000 3.99992 5.99995 7.99998 10.00000 j = 1 1.00000 2.99988 4.99992 6.99998 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(1,3) = -omega * [ AS(1,3) * u(1,2) + AW(1,3) * u(0,3) + AE(1,3) * uOld(2,3) + AN(1,3) * uOld(1,4) ] / AP(1,3) + (1 - omega) * uOld(1,3) u(1,3) = (1.17100) [ -(1.00000)(3.99992) - (1.00000)(3.00000) - (1.00000)(6.99993) - (1.00000)(6.00000) ] / (-4.00000) + (-0.17100)(4.99987) = 4.99998 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99998 6.99993 8.99998 11.00000 j = 2 2.00000 3.99992 5.99995 7.99998 10.00000 j = 1 1.00000 2.99988 4.99992 6.99998 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(2,3) = -omega * [ AS(2,3) * u(2,2) + AW(2,3) * u(1,3) + AE(2,3) * uOld(3,3) + AN(2,3) * uOld(2,4) ] / AP(2,3) + (1 - omega) * uOld(2,3) u(2,3) = (1.17100) [ -(1.00000)(5.99995) - (1.00000)(4.99998) - (1.00000)(8.99998) - (1.00000)(8.00000) ] / (-4.00000) + (-0.17100)(6.99993) = 6.99998 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99998 6.99998 8.99998 11.00000 j = 2 2.00000 3.99992 5.99995 7.99998 10.00000 j = 1 1.00000 2.99988 4.99992 6.99998 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 Calculating u(3,3) = -omega * [ AS(3,3) * u(3,2) + AW(3,3) * u(2,3) + AE(3,3) * uOld(4,3) + AN(3,3) * uOld(3,4) ] / AP(3,3) + (1 - omega) * uOld(3,3) u(3,3) = (1.17100) [ -(1.00000)(7.99998) - (1.00000)(6.99998) - (1.00000)(11.00000) - (1.00000)(10.00000) ] / (-4.00000) + (-0.17100)(8.99998) = 9.00000 i = 0 i = 1 i = 2 i = 3 i = 4 j = 4 4.00000 6.00000 8.00000 10.00000 12.00000 j = 3 3.00000 4.99998 6.99998 9.00000 11.00000 j = 2 2.00000 3.99992 5.99995 7.99998 10.00000 j = 1 1.00000 2.99988 4.99992 6.99998 9.00000 j = 0 0.00000 2.00000 4.00000 6.00000 8.00000 At end of 9 iterations, the maximum residual is 0.00008 and the maximum relative increment is 0.00015 Residual array i = 1 i = 2 i = 3 j = 3 0.00000 0.00000 0.00000 j = 2 0.00003 0.00000 0.00000 j = 1 0.00008 0.00003 0.00000 Relative increment array i = 1 i = 2 i = 3 j = 3 0.00002 0.00001 0.00000 j = 2 0.00007 0.00003 0.00001 j = 1 0.00015 0.00006 0.00001