With the increase in dimension, we only have to add another for loop in the algorithm in the previous section, and the new formula to accommodate the existence of the additional dimension. As a result, using the formula that was derived in section, we attain the result that we desired. The result was similar to that of the 1-D solution, but the problem domain was extended, hence increasing the need to parallelize the algorithm.
The resulting images have been put together as an animation. Please see it at: http://cs.erin.utoronto.ca/ e0gagq6g/492/poisson.mpeg