Created attachment 666 [details] detailed description Value_i (integral value at the point i) is correct for i>=1, but for i=0 it is incorrect: Value_0 = q_0*dltau != 0. See attachment with detailed description.
Your suggestion will lead to a zero-value for a particle range at the lowest energy threshold (Value_0). It's not right from the perspective of radiophysics. The range should be a non-zero value because the particle has kinetic energy (= 990 eV). We approximately compute Value_0 in the way shown in that code.