function xy = MyG(z0alpha0)
% function xy = MyG(z0alpha0)
x_B = 2.00;
y_B = 3.05;
z0 = z0alpha0(1);
alpha0 = z0alpha0(2);
s0 = 30; % (x_B/z0)/cos(alpha0); % approximation
G0 = zeros(5, 1);
G0(1) = 0;    % x0
G0(2) = 1.75; % y0
G0(3) = z0;   % z0
G0(4) = s0*cos(alpha0); % vx0, alpha0 from ground plane
G0(5) = s0*sin(alpha0); % vy0
[t, G] = ode45(@DG, [0 1], G0);
% [t, G] = eulexpl(@DG, [0 1], G0, 500);
plot(G(:, 1), G(:, 2))
hold on
xy = zeros(2,1);
xy(1) = G(end, 1) - x_B;
xy(2) = G(end, 2) - y_B;
end