#! /usr/bin/octave

load 20090405.dat
[time, x, y, z, dt, dx, dy, dz, ds, hr, dist] = forerunner(data_20090405);
N = size(time);
m = 190/2.2;
vx(1)=0;
vy(1)=0;
vz(1)=0;
dv_x(1) = 0;
dv_y(1) = 0;
dv_z(1) = 0;
ax(1) = 0;
ay(1) = 0;
az(1) = 0;
ax(2) = 0;
ay(2) = 0;
az(2) = 0;
W(1)=0;
W(2)=0'
for n = 2:N(2)
vx(n) = dx(n)/dt(n);
vy(n) = dy(n)/dt(n);
vz(n) = dz(n)/dt(n);
dv_x(n) = vx(n) - vx(n-1);
dv_y(n) = vy(n) - vy(n-1);
dv_z(n) = vz(n) - vz(n-1);
endfor
for n = 3:N(2)
ax(n) = dv_x(n)/dt(n);
ay(n) = dv_y(n)/dt(n);
az(n) = dv_z(n)/dt(n);
W(n) = m*(ax(n)*dx(n) + ay(n)*dy(n) + 9.8*dz(n));
endfor
__gnuplot_set__ term png
__gnuplot_set__ output '2090405_work.png';
plot(W);
cumulative = sum(W)
greatest = max(W)