We investigate the soliton dynamics in tapered parabolic index fibers via symbolic computation for a variety of dispersion profiles to inspect how a specific dispersion profile controls the optical soliton. By means of AKNS procedure, Lax pair is constructed for nonlinear Schrodinger equation with variable coefficients. Using obtained Lax pair, multi-soliton solutions are generated via Darboux transformation technique. Using multi-soliton solutions, soliton dynamics in tapered parabolic index fiber with the hyperbolic, Gaussian, exponential, and linear profiles are discussed. Results obtained in this study will be of certain potential application on construction of the nonlinear optical devices by soliton control. Results obtained in this study will be of certain value to the studies on the propagation and application of the soliton in the tapered parabolic index fiber and dispersion-managed fiber system.