Numerical simulation of two-soliton and three-soliton molecules evolution in passively mode-locked fiber laser