function [t,y] = reech_stride(x,fe);

% [t,y] = reech_stride(x,fe);
% x time series
% fe sampling frequency (the sampling period should be about half of the
% time between two strides)
%t time vector
% y regularly resampled signal

v = [0 ; cumsum(x)];
stime = 0.5*[v(2:end)+v(1:end-1)];
t = 0:1/fe:v(end);
y = spline(stime,x,t);
t = t(:);
y = y(:);

