function res = bvpfun(ya,yb)
% function res = bvpfun(ya,yb)
% Kartik's MATLAB code
% November 19, 2007
% User file that computes the residual in the boundary conditions

res = [ya(1)-17; yb(1)-43/3];