The attached worksheet develops a procedure for extrapolating boundary data from a square to its interior. Specifically, let's consider the square [−1,1]×[−1,1] and the continuous functions u1, u2, u3, u4 defined on its edges. The procedure constructs a continuous function u(x,y) in the interior of the square which matches the boundary data.
The function u(x,y) is necessarily discontinuous at a corner if the boundary data of the edges meeting at that corner are inconsistent. However, if the boundary data are consistent at all corners, then u(x,y) is continuous everywhere including the boundary.
Here is what the extrapolating function looks like: