Procedure from library rinvar.lib (see rinvar_lib).
rinvar.lib
ReynoldsImage(RO, f); list RO, poly f
compute the Reynolds image of the polynomial f where RO represents the Reynolds operator
poly