UTSx := UTS(FRAC INT,x,0)

(1)

monx2 := monomial(1,2)$UTSx

(2)

sqrt(monx2)

(3)

The last line should of course be and not . Clearly, sqrt() shifts the coefficients of the result such that the leading term is always O(1).

I'm forwarding this fix proposed by Igor himself, lest it is lost in the archives.

The relevant message in axiom-math is

http://lists.gnu.org/archive/html/axiom-math/2006-08/msg00025.html

For simplicity, here is the fix:

  diff -ur axiom-20050901/src/algebra/sttaylor.spad.pamphlet axiom-20050901.new/src/algebra/sttaylor.spad.pamphlet
  --- axiom-20050901/src/algebra/sttaylor.spad.pamphlet 2005-01-30 07:03:10.000000000 -0500
  +++ axiom-20050901.new/src/algebra/sttaylor.spad.pamphlet     2006-08-21 15:49:25.000000000 -0400
  @@ -448,6 +448,9 @@

             RATPOWERS => co**rn * YS(powerrn(rn,(invCo :: A) * x,#1))
             error "** rational power of coefficient undefined"
  +        zeropref : ST A := new(ord::I::NNI,0::A)$LA :: ST A
  +        concat(zeropref, power)
  +

       if A has Field then
         mapdiv(x,y) == delay

If somebody could complete the documentation as outlined in

http://lists.gnu.org/archive/html/axiom-developer/2006-08/msg00486.html

that would be great.

Finally, I have to confess that it is very likely myself who introduced this bug, while fixing #63

Martin

Status: open => fix proposed

Here is the documented patch. (Documentation is partially due to Waldek Hebisch)

powern.patch

Martin

Fixed in OpenAxiom