clc; clear all; close all;
g = -9.8; 
v = 5; 
theta = 30 * pi/180; % launch angle, radians
x=0;
y=0;
vx=v*cos(theta);
vy=v*sin(theta);
dt=0.01;
plot(x,y,'o')
axis([-1 5 0 0.35])
%hold on
while(1)
    x=x+vx*dt;
    y=y+vy*dt;
    vy=vy+g*dt;
    plot(x,y,'o')
    axis([-1 5 0 0.35])
    pause(0.1)
    if(y<0)
      vy=vy-g*dt;
      vy=-vy;  
    end
    if(x>5)
    vx=-vx;
    end
    if(x<0)
    vx=-vx;
    end
end