Movimiento de un proyectil (Matlab)

Definimos una función:

function[hmax,dmax]=trayectoria(v0,te)
% Velocidad inicial v0, ángulo te
g=9.81;%Gravedad
v0x=v0*cos(te*pi/180); %Componente de la Velocidad inicial
v0y=v0*sin(te*pi/180);%Componente de la Velocidad inicial
thmax=v0y/g; %Tiempo en que el Proyectil alcanza el punto más alto
hmax=v0^2/(2*g); %Altura máxima
ttot=2*thmax;%Tiempo total
dmax=v0x*ttot;%Distancia máxima
tplot=linspace(0,ttot,200);%Vector tiempo de 200 elementos
x=v0*tplot;
y=v0y*tplot-0.5*g*tplot.^2;%Posición del proyectil
plot(x,y) %Gráfico

Llamamos a la función:

[h,d]=trayectoria(230,45)

Y se muestra la trayectoria del proyectil: