function y = bcfun(ya,yb,BV)
y =[ya(1)-BV.alpha;yb(1)-BV.beta];
