1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84
|
#
# jython examples for jas.
# $Id$
#
from jas import Ring
from edu.jas import application
rs = """
# polynomial ring:
Rat(x1,x2,x3,y1,y2) G|3|
""";
ps = """
(
( y1 + y2 - 1 ),
( x1 - y1^2 - y1 - y2 ),
( x2 - y1 - y2^2 ),
( x3 - y1 y2 )
)
""";
r = Ring( rs );
print "Ring: " + str(r);
i = r.ideal( ps );
print "Ideal: " + str(i);
g = i.GB();
print "seq GB:", g;
rsi = """
# polynomial ring:
Rat(x1,x2,x3) G
""";
ri = Ring( rsi );
print "Ring: " + str(ri);
y = application.Ideal(g.pset).intersect(ri.ring);
len = y.list.size();
print "seq intersect y: ", y;
rs = """
# polynomial ring:
Rat(y1,y2,x1,x2,x3) G|2|
""";
ps = """
(
( y1 + y2 - 1 ),
( x1 - y1^2 - y1 - y2 ),
( x2 - y1 - y2^2 ),
( x3 - y1 y2 )
)
""";
r = Ring( rs );
print "Ring: " + str(r);
i = r.ideal( ps );
print "Ideal: " + str(i);
g = i.GB();
print "seq GB:", g;
rsb = """
# polynomial ring:
Rat(y1,y2) G
""";
rb = Ring( rsb );
print "Ring: " + str(rb);
print
y = application.Ideal(g.pset).intersect(rb.ring);
len = y.list.size();
print "seq intersect y: ", y;
print
|