clc;
clear all;
close all;
g=9.8;
teta=30;
v0=30;
trad=teta*pi/180;
xmax=v0^2*sin(2*trad)/g;
dx=xmax/50;
x=0:dx:xmax;
y=x*tan(trad)-g*x.^2/(2*v0^2*(cos(trad))^2);
plot(x,y);
