Dynamics of the river plume is a very complicated non-linear problem with the free boundary changing in time and space. Mixing with the ambient water through the boundary makes the problem more complicated. In this paper we reduced 3-dimensional problem into 1-dimensional one by using the integral analysis method. Basic equations have been integrated over the lateral and vertical variations. For these integrations we adopted the well-established assumption that the flow-axis component of plume velocity and the density difference of the plume with the ambient water have Gaussian distributions in directions which are perpendicular to the flow-axis of the plume. We also used the result of our previous study on the lateral spreading velocity of the plume derived under the same assumption. And entrainment was included as a mixing process. The resultant 1-dimensional equations were solved by Runge-Kutta numerical method. Consequently, comparatively easy method of numerical analysis is presented for the 3-dimensional river plume. The method can also be used for the analysis of the thermal plume of cooling water of power plants.