Changes http://wiki.axiom-developer.org/338EvalMissesMainVariable/diff
--
??changed:
\begin{axiom}
-g := operator 'g
+g := operator 'g;
eval(g(1783), g, i+->i)
??changed:
smprep(lop, lexp, lfunc, p) ==
- (v := mainVariable p) case "failed" => p::%
- symbolIfCan(k := v::K) case SY => p::%
- g := (op := operator k)
- (arg := [eval(a,lop,lexp,lfunc) for a in argument k]$List(%))
- q := map(eval(#1::%, lop, lexp, lfunc),
- univariate(p, k))$SparseUnivariatePolynomialFunctions2(MP, %)
- (n := position(name op, lop)) < minIndex lop => q g
- a:% := 0
- f := eval((lfunc.n) arg, lop, lexp, lfunc)
- e := lexp.n
- while q ^= 0 repeat
- m := degree q
- qr := divide(m, e)
- t1 := f ** (qr.quotient)::N
- t2 := g ** (qr.remainder)::N
- a := a + leadingCoefficient(q) * t1 * t2
- q := reductum q
- a
+ ..(v := mainVariable p) case "failed" => p::%
+ ..symbolIfCan(k := v::K) case SY => p::%
+ ..g := (op := operator k)
+ .....(arg := [eval(a,lop,lexp,lfunc) for a in argument k]$List(%))
+ ..q := map(eval(#1::%, lop, lexp, lfunc),
+ ...........univariate(p, k))$SparseUnivariatePolynomialFunctions2(MP, %)
+ ..(n := position(name op, lop)) < minIndex lop => q g
+ ..a:% := 0
+ ..f := eval((lfunc.n) arg, lop, lexp, lfunc)
+ ..e := lexp.n
+ ..while q ^= 0 repeat
+ ....m := degree q
+ ....qr := divide(m, e)
+ ....t1 := f ** (qr.quotient)::N
+ ....t2 := g ** (qr.remainder)::N
+ ....a := a + leadingCoefficient(q) * t1 * t2
+ ....q := reductum q
+ ..a
--
forwarded from http://wiki.axiom-developer.org/[EMAIL PROTECTED]
