close all


dL = (diff(x_s).^2 + diff(z_s).^2).^.5;
L = diff(s);

n = 1049;

e = (dL(n,:)-L(n,:))./L(n,:);


T = 11675*e.*e-25086*e.*e.*e+19235*e.*e.*e.*e;

figure
plot(snap_dt*(1:length(e)),e);


figure
plot(e,T,'x')
hold all
plot(e,T_s(n,1:end)*4.45,'o');
plot(e,5e4*e,'+');


