% Convection - Diffusion - 2D; Dirichle BC on [0,1]*[0,1]
% compute an exact solution sol=x^3*y^2*(1-x)*(1-y)
%
function sol=cd_sol(Node)

nnode=size(Node,2);

for k=1:nnode,
    xcoor = Node(1,k);
    ycoor = Node(2,k);
    sol(k,1)=(xcoor^3)*(ycoor^2)*(1-xcoor)*(1-ycoor);
end

%
