con2d.m

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

% scaled nonsingular general solution of convection-diffuison problems

 function  Con2d=Con2d(scale,m,Lr,detax,detay)
 
 Con2d=exp(-scale+detax+detay+Lr)*besseli(m,Lr,1);