The kinetic modeling of falling film on the horizontal tube is derived using integral boundary layer equation with more realistic velocity distribution. In order to improve the accuracy, the surface tension, gravity force, inertia force and shear stress have been included. The obtained equations are solved numerically. The results showed that compared with Nusselt theory, the equation with surface tension is more accordant with experimental results at same Re. Film thickness decreases first and then increases, and the angle of minimum thickness at circumferential position is bigger than 90?. The increment of flow rate, decrement of shear stress coefficient and diameter can enhance the angle of minimum thickness and the change of parameters will not affect the distribution of liquid film. The angle of minimum thickness will increase when the flow rate increases, the shear stress coefficient decreases or the diameter decreases. Moreover, the surface tension plays an important role on the liquid film flow.