%rackete on a target
% define boundary conditions u(a)=1 and u(b)=1 
% writen in the form g(y(a),y(b))=0, residual
function res = bvp_bc_def(ya,yb)
    res = [ ya(1) - 1	
	    yb(1) - 1 ];
return;