1.2.53 `sm1.saturation'
-----------------------

sm1.saturation(T)
     :  T = [I,J,V].   It returns saturation of I with respect to
     J^infty.    V is a list of variables.
Example:
      sm1.saturation([[x2^2,x2*x4, x2, x4^2], [x2,x4], [x2,x4]]);

