idx = find(pos < 1E-8, 1, 'first');
endTime = t(idx);
plot(t(1:idx), pos(1:idx), '--','DisplayName','Verlet solution');
legend('show')
xlim([0 endTime])
ylim([0 height])
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!