This paper presents a digital modeling technique of the distributed system. The basic idea of the proposed technique is to discretize a continuous system with respect to the spatial coordinates using bilinear method. The response of the discretized system is analyzed by Laplace transform and z-transform. The computational results in torsional shaft and Timoshenko beam using the proposed technique are compared with the exact solutions and the results of finite element method.