The Saint-Venant equation is a mathematical model which could be used to study water flow in an open channel, river, etc. The Kurganov-Petrova (KP) method, which is a second-order scheme, is used to solve the Saint-Venant equations with good stability. The water flow of a river between two hydropower stations in Norway has been simulated in this study using MATLAB and OpenModelica. The KP scheme has been used to discretize the Saint-Venant equations in the spatial domain, yielding a collection of Ordinary Differential Equations (ODEs). These are then integrated with time using the variable steplength solvers in MATLAB: ode23t, ode23s, ode45, and fixed step-length solvers: The Euler method, the second and fourth order Runge Kutta method (RK2 and RK4). In OpenModelica built-in, variable step-length DASSL solver has been used. From the simulation, it was observed that all solvers produce more or less similar results. Volumetric flowrate calculation indicated numerical oscillation with variable step-length solvers in MATLAB. The results indicated that it is reasonable to match the order of space and time discretization.