Non-oscillatory, non-free parameter scheme is employed to split three-dimensional,time-dependent hydrodynamic model for numerical computation. Computations ofbackground fluid field and disturbance propagation are carried out for the coronalmass ejection occurring in May 1998. A similar density distribution to that of radial magnetic on solar source surface can be determined from the HAO K-coronameter polarization-brightness data, which can be taken as inner boundary densitydistribution. The propagation of disturbance shows a sharp discontinuity withoutspurious oscillation and dispersion. Also, the numerical test shows that the schemekeeps a stable iteration without any artificial diffusion being added. On the otherhand, the source code programmed with this hydrodynamic model can be used tomodel the specific coronal mass ejection when the configuration of magnetic field isnot decisive the solar wind flow.