For scientific engine analyses, authors proposed a simultaneous and approximate analyses for connecting rod big end bearings, for piston pin bearings and for crank shaft main bearings. However, to save computation time for engine bearing analysis, it is better to innovate new iteration methods to introduce solutions of partial differential equations in finite width bearing theory. In this paper, some kinds of Runge-Kutta methods are proposed by using variable step methods with inherent algorithm of engine bearings, and the characteristic of variable step's Runge-Kutta methods in engine bearing analyses are clarified.