An optimized finite difference Crank-Nicolson iterative scheme for the 2D Sobolev equation