In this paper, an analytical and numerical computation of multi-solitons in Korteweg-de Vries (KdV) equation is presented. The KdV equation, which is classic of all model equations of nonlinear waves in the soliton phenomena, is described. In the analytical computation, the multi-solitons in KdV equation are computed symbolically using computer symbolic manipulator<span style="white-space:nowrap;">—</span>Wolfram Mathematica via Hirota method because of the lengthy algebraic computation in the method. For the numerical computation, Crank-Nicolson implicit scheme is used to obtain numerical algorithm for the KdV equation. The simulations of solitons in MATLAB as well as results concerning collision or interactions between solitons are presented. Comparing the analytical and numerical solutions, it is observed that the results are identically equal with little ripples in solitons after a collision in the numerical simulations;however there is no significant effect to cause a change in their properties. This supports the existence of solitons solutions and the theoretical assertion that solitons indeed collide with one another and come out without change of properties or identities.