Código simulación futbol-3D

function mov_parabolico_3d(xo,zo,vo, alpha)

close

%%declaración del vector de tiempo de simulacion
%%dividimos la unidad de tiempo en fs

yo=0.1;
fs=10;
t=0:1/fs:10;

%Velocidades por cada eje

vox=vo*cos(20*pi/180);
voy=vo*sin(20*pi/180);
voz=2*vo*sin(alpha*pi/180);

%%declaracion de las ecuaciones de modelado del sistema

Sx=vox*t+xo; %%Vector Posicion en X con respecto al tiempo

Sy=-(1/2)*9.81*t.^2+voy*t+yo; %%Vector posicion en Y con respecto al tiempo

Sz=voz*t+zo; %%Vector Posicion en X con respecto al tiempo

Symax=find(Sy<=0);

%Dibujar arco futbol

plot3([-100 -100 100 100],[300 300 300 300],[0 20 20 0] ,'-r','LineWidth',3)
hold on
patch([-300 -300 300 300],[0 300 300 0],[0 0 0 0],'g')
hold on
plot3([-200 -200 200 200],[300 200 200 300],[0 0 0 0] ,'-w','LineWidth',2)
%% configuración grafica para la animacion

for i=1:Symax(1)
    plot3(Sz(i),Sx(i),Sy(i),'ok');
    axis([-300,300,0,300,0,100]);
    xlabel('Sz(t)');
    ylabel('Sx(t)');
    zlabel('Sy(t)');
    pause(1/fs);
    grid on;
    hold on;
end

end