clc;
clear all;
close all;
g=9.8;
teta=5:5:70;
nsize=size(teta');
v0=80;
trad=teta*pi/180;
for i=1:nsize
    xmax=v0^2*sin(2*trad(i))/g;
    dx=xmax/50;
    x=0:dx:xmax;
    y=x*tan(trad(i))-g*x.^2/(2*v0^2*(cos(trad(i)))^2);
plot(x,y);
title('kutxit gasrola');
xlabel('x');
ylabel('y');
grid on;
hold on;
end


