Thank you @vv . Yes, Indeed. There are too many solutions and some of them are identical. I just learnt that the physical solutions are expected to sit on the two lines Im(x_j) ± Pi/2. Therefore, how to avoid physically meaningless solutions and identical solutions?

@Thomas Richard  it is done.