The Sagdeev pseudo-potential approach has been employed extensively in theoretical studies to determine large-amplitude (fully) nonlinear solutions in a variety of multi-species plasmas. Although these solutions are repeatedly considered as solitary waves (and even solitons), their temporal stability has never been proven. In this paper, a numerical study of the Vlasov-Poisson system is made to follow their temporal evolution in the presence of numerical noise and thereby test their longtime propagation stability. Considering the ion-acoustic regime, both constituents of the plasma, i.e. electrons and ions are treated following their distribution functions in these set of fully-kinetic simulations. The findings reveal that the stability of Sagdeev solution depends on a combination of two parameters, i.e. velocity and trapping parameter. It is shown that there exists a critical value of trapping parameter for both fast and slow solutions which separates the stable from unstable solutions. In case of stable solutions, it is shown that these nonlinear structures can propagate for long periods, which confirms their status as solitary waves. Stable solutions are reported for both Maxwellian and Kappa distribution functions. For unstable solutions, it is demonstrated that the instability causes the Sagdeev solution to decay by emitting ion-acoustic wave-packets on its propagation trail. The instability is shown to take place in a large range of velocity and even for Sagdeev solutions with velocity much higher than ion-sound speed. Besides, in order to validate our simulation code two precautionary measures are taken. Firstly, the well-known effect of the ion dynamics on a stationary electron hole solution is presented as a benchmarking test of the approach. Secondly, In order to verify the numerical accuracy of the simulations, the conservation of energy and entropy are presented. *