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


