A mathematical model of a gravitational flow of a liquid film above a fixed bed was obtained based on the fundamental Navier-Stokes formulation. It results a pair of non-linear partial differential equations denominated Shallow Waters. This PDE system was discretized with the MacCormack’s method and later implemented in FORTRAN. The results are compared with the Linear Stability Analysis solution and with other analytical studies.