function w = ODEFunfromShooting(x,y,Fun, dyFun,dypFun)
%
% function used in ODE solver for the nonlinear shooting method.
%
y1    = y(1:2);
fun   = feval(Fun,x,y1);
dyfun = feval(dyFun,x,y1);
dypfun= feval(dypFun,x,y1);
w     = [y(2);fun;y(4);dyfun*y(3)+dypfun*y(4)];
