square_eq.m

％ A 2D homogeneous convection-diffusion case (u=exp(-ex*deta*x-ex*deta*y) with a square with ％ all

% Uniform node data for sqaure domain [0,1]*[0,1]

function [x,y]=square_eq(a,b,n)

     x=0:a/(n-1):a;
     y=(b/a)*x;