init_rms: L-2 norm of the residual at the start of iterations
den = -2.0/dx^2 - 2.0/dy^2
k=0
for k = 1:max_iter