function f = shoot_chemical ( alpha )
%
NSTEP = 100;
[x,y] = rk3 ( 'chemical_ode', [ 0.0, 1.0 ], [ 0; alpha ], NSTEP );
f = y(1,NSTEP+1);


