sort output of roots or residue so that the repeated roots occur together

The output of the roots command come out in complex conjugate pairs.

e.g.

roots([1 8 38 88 121])

ans =

-2.0000 + 2.6458i

-2.0000 - 2.6458i

-2.0000 + 2.6458i

-2.0000 - 2.6458i

I want it to come out as:

-2.0000 + 2.6458i

-2.0000 + 2.6458i

-2.0000 - 2.6458i

-2.0000 - 2.6458i

so that I can algorithmically recognize repeated roots.

(The roots command used to do this. The current output is not an improvement>)

Please advise!

Dr. Vallorie Peridier, Temple University

Answer 1

[KALYAN ACHARJYA]

>> data=roots([1 8 38 88 121])

data =
  -2.0000 + 2.6458i
  -2.0000 - 2.6458i
  -2.0000 + 2.6458i
  -2.0000 - 2.6458i

>> result=sort([real(data) imag(data)],'descend')*[1;1i]
result=

  -2.0000 + 2.6458i
  -2.0000 + 2.6458i
  -2.0000 - 2.6458i
  -2.0000 - 2.6458i

>>