The file "hrv2631.m" is a matlab function to calculate human response to vibration. I understand that was 'working' version.
The file "sample-acc.txt" is an ascii file (ex Linux) which contains vertical (z) acceleration-time history at 200 samples per sec. The accelerations are in g units, hence will require conversion to m/sec^2 before being passed to hrv2631.
My translation of hrv2631.m into an Rlab2 function produces the following results using the data in "sample-acc.txt" (the Rlab2 version takes an additional parameter denoting units of the signal and converts 'g' values to m/sec^2 if that is what is supplied);

RMS = 0.714
VDV = 1.8
TVDV15 = 6.63e+4