function [X,Y,U] = postprocU(x,u,m,n)

% [X,Y,U] = postprocU(x,y,u,m,n)
%            It creates coordinate matrices X, Y
%            and the corresponding solution matrix U,
%            each being m by n.

X = zeros(m,n);
Y = zeros(m,n);
U = zeros(m,n);

index = 1;

for i=1:m
    for j=1:n
        X(i,j) = x(1,index);
        Y(i,j) = x(2,index);
        U(i,j) = u(index);
        
        index = index + 1;
    end
end
