+ args = l;
+ }
+ // Find the symbol to factor in at this stage
+ if (!is_ex_of_type(args.op(0), symbol))
+ throw (std::runtime_error("sqrfree(): invalid factorization variable"));
+ const symbol x = ex_to_symbol(args.op(0));
+ // convert the argument from something in Q[X] to something in Z[X]
+ numeric lcm = lcm_of_coefficients_denominators(a);
+ ex tmp = multiply_lcm(a,lcm);
+ // find the factors
+ exvector factors = sqrfree_yun(tmp,x);
+ // construct the next list of symbols with the first element popped
+ lst newargs;
+ for (int i=1; i<args.nops(); ++i)
+ newargs.append(args.op(i));
+ // recurse down the factors in remaining vars
+ if (newargs.nops()>0) {
+ for (exvector::iterator i=factors.begin(); i!=factors.end(); ++i)
+ *i = sqrfree(*i, newargs);
+ }
+ // Done with recursion, now construct the final result
+ ex result = _ex1();
+ exvector::iterator it = factors.begin();
+ for (int p = 1; it!=factors.end(); ++it, ++p)
+ result *= power(*it, p);
+ // Yun's algorithm does not account for constant factors. (For
+ // univariate polynomials it works only in the monic case.) We can
+ // correct this by inserting what has been lost back into the result:
+ result = result * quo(tmp, result, x);
+ return result * lcm.inverse();