s="comment">% cornermod = true;</span>
0044 <span class="comment">% res = 50;</span>
0045 <span class="comment">% close all; warning off;</span>
0046 
0047 <span class="comment">% generate mesh</span>
0048 node = [0 0; lx 0; lx ly ; 0 ly];
0049 edge = [1 2 ; 2 3 ; 3 4; 4 1]; 
0050 hdata.hmax = res; 
0051 options.output = true;
0052 [p,t] = mesh2d(node,edge,hdata,options);
0053 
0054 <span class="comment">% generate connectivity</span>
0055 [e,te,et2,bnd] = connectivity(p,t);
0056 
0057 
0058 <span class="co