%

function [pl,ql,pr,qr] = bc_Neumann(xl,ul,xr,ur,t)

pl = 0;
ql = 1;
pr = 0;
qr = 1;

end