recnbc2

 rectangle with Neumann BC, 2 components

0001 function [geo, bc]=recnbc2(lx,ly) 
0002 % rectangle with Neumann BC, 2 components
0003  geo=rec(lx,ly);
0004  bc=[[  2     2     2     2]; % Neumann BC, exported from gui
0005      [0     0     0     0];
0006      [1     1     1     1];
0007      [1     1     1     1];
0008      [1     1     1     1];
0009      [1     1     1     1];
0010      [1     1     1     1];
0011      [1     1     1     1];
0012     [48    48    48    48];
0013     [48    48    48    48];
0014     [48    48    48    48];
0015     [48    48    48    48];
0016     [48    48    48    48];
0017     [48    48    48    48];
0018     [48    48    48    48];
0019     [48    48    48    48];
0020     [48    48    48    48];
0021     [48    48    48    48];
0022     [48    48    48    48];
0023     [48    48    48    48];
0024     [49    49    49    49];
0025     [48    48    48    48];
0026     [48    48    48    48];
0027     [49    49    49    49];
0028     [48    48    48    48];
0029     [48    48    48    48]];
0030 end