clc
clear all
for i=1:100
    [a(i) b(i) c(i)] = main();
    a=a*9.81
    b=b*9.81;
    c=c*9.81;
    plot(i,a(i),'--rs','LineWidth',2,...
                'MarkerEdgeColor','k',...
                'MarkerFaceColor','g',...
                'MarkerSize',10)
    ylim([-9.8 9.8]);
    xlim([0 100]);
    hold on
    pause(0.01)
end