% A 2D homogeneous convection-diffusion case (u=exp(-ex*deta*x-ex*deta*y) with a square with % all Dirichlet boundary, note that reaction coefficient is not zero % by indirect BKM