function w = ODEFunfromLinearShooting(x,y,pFun, qFun,rFun)
%
% function used in ODE solver for the nonlinear shooting method.
%
pfun   = feval(pFun,x);
qfun   = feval(qFun,x);
rfun   = feval(rFun,x);
w     = [y(2);pfun*y(2)+qfun*y(1)+rfun;y(4);pfun*y(4)+qfun*y(3)];

