BC for geometry with six edges, mixed DBC and NBC but fixed coeff
0001 function bc=ac6bcf(p,u,lam) 0002 % BC for geometry with six edges, mixed DBC and NBC but fixed coeff 0003 sf=10^4; % stiff spring factor... 0004 qd=sf; gd = sf*lam; % Dirichlet part 0005 qn = zeros(p.neq);gn = zeros(p.neq,1); % Neumann part 0006 % ordering in geometry from rec: bottom, right, top, left 0007 bc = gnbc(p.neq,qn,gn,qn,gn,qn,gn,qn,gn,qd,gd,qn,gn);