convergence for 2D case is O(h) not O(h^2)