The dynamic simulation of slider in hard disk drive is performed using Factored Implicit Finite Difference method. The modified Reynolds equation with Fukui and Kaneko model is employed as a governing equation. Equations of motion for the slider of three degrees of freedom are solved simultaneously with the modified Reynolds equation. The transient responses of the slider for disk step bumps and slider impulse forces are shown for various cases and are compared for the iteration algorithm and new algorithm.