function bc = bcs(ya,yb)
bc = [ya(1); yb(1)];