function f = chemical_shoot ( alpha )
%
%  function f = chemical_shoot ( alpha )
%
%  CHEMICAL_SHOOT uses the shooting method to determine an approximate
%  solution to a particular boundary value problem.
%
%  Modified:
%
%    11 February 2000
%
%  Author:
%
%    John Burkardt
%
%  Parameters:
%
%    Input, real ALPHA, the value of Y' at the left endpoint.
%
%    Output, real F, the error, that is, the difference between the
%    compute value of the solution at the right endpoint, and the
%    desired value, determined by the boundary condition:
%
%      Y(RIGHT) = 0
%
%    so
%
%      F = Y(RIGHT).
%
nstep = 100;
[x,y] = rk3 ( 'chemical_ode', [ 0.0, 1.0 ], [ 0.0; alpha ], nstep );
f = y(1,nstep+1);


