I am trying to use the reduce() command on a rational pollynomial. I
first clear the denominator by multiplying by the denominator but when
I use reduce() I get an error. When I print out the polynomial it is
no longer rational but SAGE doesn't like it. I would appriciate any
advice. Below is an example script with ouput. Thank you.
David Stahl
Untitled
system:sage
{{{id=0|
R1 = PolynomialRing(RationalField(),9,
["x0","x1","x2","y0","y1","y2","a0","a1","a2"], "lex")
x0,x1,x2,y0,y1,y2,a0,a1,a2=R1.gens()
X=[
x0*y0 - x2*y1 + x1*y2,
a2*y0 + a1*y1 - a0*y2,
a1*x0 - a0*x1 - a2*x2
]
I=R1.ideal(X)
I2 = I.groebner_basis()
print X[2].reduce(I2)
Y=X[2]/a1
Z=a1*Y
print Z
print Z.reduce(I2)
///
0
-x1*a0 + x0*a1 - x2*a2
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/home/notebook/sage_notebook/worksheets/admin/61/code/3.py",
line 17, in <module>
exec compile(ur'print Z.reduce(I2)' + '\n', '', 'single')
File "/usr/local/sage/data/extcode/sage/", line 1, in <module>
TypeError: reduce() takes exactly 1 argument (2 given)
}}}
{{{id=1|
}}}
--~--~---------~--~----~------------~-------~--~----~
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/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---