I am solving the KdV equation using Hirota's method. I used the perturbation series to solve for the one-, two-, three-, and N-soliton solutions but I am having some challenges in the substitution of the D-operator in the series to solve for the unknown functions.