r:=rule(j*j==1)   
      
   Change the definition of * so every   
   x*x it reach it Return 1   
      
   so 3*3 is 1 ok   
      
   Possibly could be ok too if it is possible specific one set as in   
      
   versors:=[[sin(x),cos(x)] for x in [0.1, 0.2,0.3, 0.4,0.5, 0.6, 0.7, 0.8]   
      
   r:=rule(sqrt(a^2+b^2)==1|[a,b] in versors)   
      
   Or define a generic abstract set of Versors   
      
   r:=rule(sqrt(a^2+b^2)==1|[a,b] in Versors)   
      
   That has the rule   
   ----   
      
   But how to define a rule just for a symbol j   
   I tried rule(%j*%j==1) but seems to me in that rule not only j*j=1 but even   
   x*x=1   
      
   --- SoupGate-Win32 v1.05   
    * Origin: you cannot sedate... all the things you hate (1:229/2)