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

pl=[ul(1), ul(2)];
ql=[0,0];

pr=[ur(1), ur(2)]; 
qr=[0,0]; 
