The recursive formulas of modified Bessel functions give the relativistic expressions for energy and momentum. Modified Bessel functions are solutions to a continuous time, one-dimensional discrete jump process. The jump process is analyzed from two inertial frames with a relative constant velocity; the average distance of a particle along the chain corresponds to the distance between two observers in the two inertial frames. The recursion relations of modified Bessel functions are compared to the 'k calculus' which uses the radial Doppler effect to derive relativistic kinematics. The Doppler effect predicts that the frequency is a decreasing function of the velocity, and the Planck frequency, which increases with velocity, does not transform like the frequency of a clock. The Lorentz transformation can be interpreted as energy and momentum conservation relations through the addition formula for hyperbolic cosine and sine, respectively. The addition formula for the hyperbolic tangent gives the well-known relativistic formula for the addition of velocities. In the non-relativistic and ultra-relativistic limits the distributions of the particle's position are Gaussian and Poisson, respectively.