This paper describes the generation and propagation of internal solitary wave in a two-layer fluid system by numerical analysis. Characteristics of interfacial soliton such as wave type, wave height, wave celerity are investigated numerically with respect to an extent of initial disturbance, fluid thicknesses of the two fluids and etc. The difference between the internal wave propagation on sloping beach and flat bottom was also examined. Laboratory experiments were conducted in the wave flume and compared with the results of numerical computation for verification.