function bc = fbc(y0,y1)
bc = [y0(1)+1; y1(1)-2];