This paper is comparing numerical schemes for a differential equation with convection and fourth-order diffusion. Our model equation is $h_t+(h^2-h^3)_x=-(h^3h_{xxx})_x$, which arises in the context of thin film flow driven the competing effects of an induced surface tension gradient and gravity. These films arise in thin coating flows and are of great technical and scientific interest. Here we focus on the several numerical methods to apply the model equation and the comparison and analysis of the numerical results. The convection terms are treated with well known WENO methods and the diffusion term is treated implicitly. The diffusion and convection schemes are combined using a fractional step-splitting method.