% zsave = [-10:0.1:10];
% sf = [];
% for z=zsave
%     sf = [sf shootingfunction(z)];
% end
% plot(zsave,sf)
fzero(@shootingfunction,1.0)
%fzero(@shootingfunction,3.0)
%fzero(@shootingfunction,4.0)