axiom-developer
[Top][All Lists]

## Re: [Axiom-developer] groebner?

 From: Martin Rubey Subject: Re: [Axiom-developer] groebner? Date: 16 Feb 2006 20:24:34 +0100 User-agent: Gnus/5.09 (Gnus v5.9.0) Emacs/21.4

```"Raymond E. Rogers" <address@hidden> writes:

> Could someone tell me if this is wrong; or if it's right why?

It's wrong, I guess.

> list:=[p,Vr,Vt,e]
> eq1H:=((-Vr^3+Vr^2)*Vt+Vr^3-Vr^2)*p :: HDMP(list,POLY INT)

^^^^^^^^
||||||||

you are saying that the coefficients of your polynomial are polynomials. That's
probably confusing.

> (24) -> eq1H
>                3     2        3     2
>    (24)  ((- Vr  + Vr )Vt + Vr  - Vr )p

>From this output you can deduce that ((-Vr^3+Vr^2)*Vt+Vr^3-Vr^2) is interpreted
as coefficient and p as variable...

Thus the output of

> (25) -> groebner [eq1H]
> (25) ->
>    (25)  [p]

Is correct, in some sense. Although I agree that it is confusing. Try

eq1H := (((-Vr^3+Vr^2)*Vt+Vr^3-Vr^2)*p) :: HDMP(list,INT)

instead. You need the parentheses, because otherwise Axiom will coerce only p
to a HDMP, and eq1H will be of type POLY INT

Martin

```