We discuss algorithms applicable to the numerical solution of second-order ordinary differential equations by finite differences. We make particular reference to the solution of the dissipative particle dynamics fluid model, and present extensive results comparing one of the algorithms discussed with the standard method of solution. These results show the successful modeling of phase separation and surface tension in a,binary immiscible fluid mixture.