rectangle with NBC, 1 component (scalar)
0001 function [geo, bc]=recnbc1(lx,ly) 0002 % rectangle with NBC, 1 component (scalar) 0003 geo=rec(lx,ly); 0004 bc=[[ 1 1 1 1 ]; 0005 [0 0 0 0]; 0006 [1 1 1 1]; 0007 [1 1 1 1]; 0008 [48 48 48 48]; 0009 [48 48 48 48]; 0010 [48 48 48 48]; 0011 [48 48 48 48]; 0012 [49 49 49 49]; 0013 [48 48 48 48]]; 0014 end