Abstract
The nonlinear steady wave equations based on the variational principle, were solved using the Newton-Raphson method, to find the numerical solutions for both surface and internal solitary waves. The ratio of kinetic energy to total energy for surface solitary waves increased, as the ratio of wave height to the still water depth increased. Conversely, the ratio of kinetic energy to total energy for internal solitary waves increased, when the ratio of wave height to the maximum wave height is average. The solitary waves propagation also was numerically simulated, by applying the time-dependent model, with the initial conditions obtained using the present method for the large-amplitude surface and internal solitary waves.