v = 1000;
a = 30/180*pi;
[t,y] = ode45(@bullet,[0 45],[0 0 v*cos(a) v*sin(a)]);
plot(y(:,1),y(:,2));
grid