This version keeps your exact normalised residual Rnorm(B,N), solves Rnorm(B,N)=1, and also prints the scaled quantities N(B(n)-1) for odd and even n separately so the limiting constants are easier to spot.