Dynamic simulation in robotic systems can be considered as a useful tool not only for the design of both mechanical and control systems, but also for planning the tasks of robotic systems. Usually, the dynamic model suffers from discontinuities in some parts of it, such as the use of Coulomb friction model and the contact problem. These discontinuities could lead to stiff differential equations in the simulation process. In this paper, we present an algorithm that solves the discontinuity problem of the Coulomb friction model without applying any normalization. It consists of the application of an external switch that divides the integration interval into subintervals, the calculation of the friction force in the stick phase, and further improvements that enhance its stability. This algorithm can be implemented directly in the available commercial integration routines with event-detecting capability. Results are shown by a simulation process of a simple 1-DoF oscillator and a 3-DoF parallel robot prototype considering Coulomb friction in its joints. Both simulations show that the stiffness problem has been solved. This algorithm is presented in the form of a flowchart that can be extended to solve other types of discontinuity.