%Sample script file to solve a Boundary Value Problem.  In this case, we
%are using the bisection method.

%The following was used to sample the space- I need to bracket the
%solution.

%t=linspace(0,5,10);
%t=linspace(-1,1,10);
%
%for j=1:10
%    h(j)=BVPSample(t(j));
%end

%Now I know the solution is bracked: [-0.6, -0.2]
xc=bisect(@BVPSample,-0.6,-0.2,10^(-6));


%Test the solution:
%[t,y]=ode45(@BVPeqn,[0,1],[1;xc]);
%y(end,1)
%
%plot(t,y(:,1));