function res = BC_BVP1( y1, y2 )
  res(1,1) = y1(1);
  res(2,1) = y2(1) + 1;
end