open('bang_bang_double_integrator_sim')
sim('bang_bang_double_integrator_sim')

%%
figure
subplot(3,1,1)
plot(sim_u,'LineWidth',2)
axis([0 3.5 -1.1 1.1])
xlabel('')
ylabel('u')
title('')
grid on
subplot(3,1,2)
plot(sim_v,'LineWidth',2)
xlabel('')
ylabel('v')
title('')
grid on
subplot(3,1,3)
plot(sim_y,'LineWidth',2)
xlabel('Time [s]')
ylabel('y')
title('')
grid on

%%

% print -depsc2 ../figures/bang_bang_simulation.eps
% !epstopdf ../figures/bang_bang_simulation.eps