%V0VALUE: MATLAB script M-file that computes v0
%for the ballistics object.
%Author: Peter Howard
%
g = 9.81; b = .0645;
rhs = @(t,y) [y(2);-g-b*abs(y(2))*y(2)];
bc = @(yleft,yright) [yleft(1)-.39;yright(1)];
initial_partition = linspace(0,2.13,25);
initial_solution = [.39;1];
solinit = bvpinit(initial_partition,initial_solution);
sol = bvp4c(rhs,bc,solinit)
v0 = sol.y(2,1)