%BVALUE2: MATLAB script M-file that computes v0
%for the ballistics object.
%
%Author: Peter Howard
g = 9.81; b0 = 0; %b0 is initial guess at b
rhs = @(t,y,b) [y(2);-g-b*abs(y(2))*y(2)];
bc = @(yleft,yright,b) [yleft(1)-4.06;yright(1);yleft(2)];
initial_partition = linspace(0,.95,25);
initial_solution = [4.06;-1];
solinit = bvpinit(initial_partition,initial_solution,b0);
sol = bvp4c(rhs,bc,solinit);
fprintf('b = %6.4f\n',sol.parameters)