clear
v0=[10:10:60 100];
tb=[20 18.5 18 18 18 17.5 17.3];
td=[38 21.5 19.5 19.5 19 19 19];
plot(v0,tb,'k-',v0,td,'k--')
plot(v0,tb,'k',v0,td,'k')
title('time to failure vs. impact velocity')
xlabel('v0 (m/s)'), ylabel('time (us)'), grid
legend('brittle','ductile');

%v0i=min(v0):5:max(v0);
%tbi=spline(v0,tb,v0i);
%tdi=spline(v0,td,v0i);
%hold on, plot(v0i,[tb td],'k-'), hold off
