BC for geometry with six edges, mixed Neumann and Dirichlet, x-dependent
0001 function bc=ac6bcfx(p,u,lam) 0002 % BC for geometry with six edges, mixed Neumann and Dirichlet, x-dependent 0003 sf=10^4; % stiff spring factor for Dirichlet part 0004 qd = mat2str(sf);gd = [mat2str(sf*lam) '*x']; % pseudo-Dirichlet: 0005 qn = '0'; gn = '0'; % Neumann part 0006 % ordering in geometry from rec: bottom, right, top, left 0007 bc = gnbcs(p.neq,qn,gn,qn,gn,qn,gn,qn,gn,qd,gd,qn,gn); 0008