A Lax-Wendroff type semi-implicit numerical scheme is employed to numerically integrate the nonlinear one-dimensional unmagnetized plasma multi-fluid equations for ideal gases to obtain soliton solutions from an initial density hump profile. The time evolution of such solitons is studied and is...