theta = 12*pi/180;
omega = 350;
v0 = 72;

q = [ 0 ; 0 ; v0*cos(theta) ; v0*sin(theta); omega ];

allQ = [];

dt = 0.01; % Matt's rule of thumb

while q(2) >= 0 % ball above ground
    
    dq = golf(q);
    q = q + dq*dt; % Euler's method
    
    allQ = [allQ, q];
    
end

plot(allQ(1,:), allQ(2,:))