Martin Rubey wrote:
>
> Rob A Nendorf <[email protected]> writes:
>
> > I'm brand new to axiom, and need to use the WEIER package to calculate
> > the distinguished polynomial of a given power series (I need to use
> > the weierstrass function). I'm completely lost as to how to implement
> > the package.
>
> I just tried to use WEIER and was unable to do so, too. Anybody in the
> know?
>
AFAICS WEIER is buggy. The following should work:
T := TaylorSeries(Fraction(Integer))
pFI := Polynomial(Fraction(Integer))
ts := (x+y)::pFI::mT
wP := WeierstrassPreparation(Fraction(Integer))
weierstrass(y, ts)$wP
but does not. After applying the patch below it "almost" works,
that is the result is computed and contains correct coefficients,
but trying to display it leads to error due to missing entries
in stream. ATM I am investigating if this is bug in display
routines or WEIER produced bad series.
On infinite series WEIER seem to work OK:
ts1 := 1 + ts;
ts2 := recip(ts1);
ts3 := ts2 - 1;
weierstrass(y, ts3)$wP
(12) [x + O(11)]
Type: List(TaylorSeries(Fraction(Integer)))
The patch follows:
--- ../ax-build24/src/algebra/WEIER.spad 2010-10-10 17:12:17.000000000
-0400
+++ WEIER.spad 2010-10-26 11:58:05.000000000 -0400
@@ -124,7 +124,7 @@
else
p:=tp(v,sts) pretend (ST SMPS)
b:StS:=rest(((first p pretend StS)),a::NNI)
- c:=retractIfCan first b
+ c:=retractIfCan(first b)@Union(R, "failed")
c case "failed"=>_
error "the coefficient of the lowest degree of the variable should _
be a constant"
--
Waldek Hebisch
[email protected]
--
You received this message because you are subscribed to the Google Groups
"FriCAS - computer algebra system" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/fricas-devel?hl=en.
