This function will calculate all the pth roots of unity mod r and return the result as an expression sequence.

The order, p, of the root must be prime.

Note that there will always be at least the root 1.

The command with(numtheory,rootsunity) allows the use of the abbreviated form of this command.