1 | proc mptest (lv) |
---|
2 | { |
---|
3 | print(lv); |
---|
4 | write("MPfile:w /tmp/t.mp", lv); |
---|
5 | return(lv == read("MPfile: /tmp/t.mp")); |
---|
6 | }; |
---|
7 | |
---|
8 | proc mpt(lv) |
---|
9 | { |
---|
10 | lv; |
---|
11 | write("MPfile:w /tmp/t.mp", lv); |
---|
12 | def r = read("MPfile: /tmp/t.mp"); |
---|
13 | r; |
---|
14 | }; |
---|
15 | |
---|
16 | int i; |
---|
17 | int i1 = 1; |
---|
18 | int i2 = -100092; |
---|
19 | |
---|
20 | intvec iv; |
---|
21 | intvec iv1 = 1,2,3; |
---|
22 | intvec iv2 = -1,2,-3,4,-5; |
---|
23 | |
---|
24 | intmat im; |
---|
25 | intmat im1[2][3]=1,3,5,7,8; |
---|
26 | intmat im2[3][3]= -1,2,-3,4,-5; |
---|
27 | |
---|
28 | string s; |
---|
29 | string s1 = "Hello World"; |
---|
30 | |
---|
31 | list l; |
---|
32 | list l1 = i,iv,im,i1,iv1,im1; |
---|
33 | list l2 = l1, iv2, l, im2; |
---|
34 | |
---|
35 | ring r; |
---|
36 | |
---|
37 | poly p; |
---|
38 | number n; |
---|
39 | number n1 = 2; |
---|
40 | number n2 = -7; |
---|
41 | |
---|
42 | poly p1 = x + y; |
---|
43 | poly p2 = xyz - 2x3y4z5 + 3xy -4yz + 5z; |
---|
44 | poly p3 = p2^4; |
---|
45 | |
---|
46 | vector v; |
---|
47 | vector v1 = [p1, p2]; |
---|
48 | vector v2 = [p, p1,p2, p3] + p1*gen(5); |
---|
49 | |
---|
50 | ideal id; |
---|
51 | ideal id1 = p, p1, p2, p3; |
---|
52 | ideal id2 = p, p1+p2, p2+3, p1+p2+p3, p1, p2, p3; |
---|
53 | |
---|
54 | module mv; |
---|
55 | module mv1 = v, v1, v2; |
---|
56 | module mv2 = v1+v2, p1*v2, p2*(v1-v2), v, v1, v2; |
---|
57 | |
---|
58 | matrix m; |
---|
59 | matrix m1[2][3] = p, p1, p2, id1; |
---|
60 | matrix m2[4][4] = p, p1, p2, p3, id1, id2; |
---|
61 | |
---|
62 | |
---|
63 | ideal j = std(id1); |
---|
64 | |
---|
65 | qring q = j; |
---|
66 | poly p11 = x + y; |
---|
67 | poly p12 = xyz - 2x3y4z5 + 3xy -4yz + 5z; |
---|
68 | poly p13 = p12^4; |
---|
69 | |
---|
70 | ring rtm = 32003, (x,y,z,t), dp; |
---|
71 | ideal ex = t3y7z4 + 9x8 + 5t2x4y2 + 2t2xy2z3, |
---|
72 | 2t2x5y4+7txy6+9y8 + 2t2x2yz3, 3t4x2y3z2 + 9z8, 3x2y9 + |
---|
73 | y9 + 5x4; |
---|
74 | |
---|
75 | ring rt0 = 0, (x,y,z,t), dp; |
---|
76 | ideal ex = t3y7z4 + 9x8 + 5t2x4y2 + 2t2xy2z3, |
---|
77 | 2t2x5y4+7txy6+9y8 + 2t2x2yz3, 3t4x2y3z2 + 9z8, 3x2y9 + |
---|
78 | y9 + 5x4; |
---|
79 | |
---|
80 | |
---|
81 | ring rr = 32003,(a,b),dp; |
---|
82 | map f = r, a,b,a+b; |
---|
83 | map g = rr,a2,b2; |
---|
84 | map phi = g(f); |
---|
85 | |
---|
86 | |
---|
87 | ring r0 = 0, x, lp; |
---|
88 | number n; |
---|
89 | number n1 = 29734481274863241234589; |
---|
90 | number n2 = n1/(n1-6); |
---|
91 | |
---|
92 | poly p; |
---|
93 | poly p1 = n1*x + n2*x^4 + 7/3*x^2*x^5 - 6; |
---|
94 | poly p2 = x*x*x - 2*x^3*x^4*x^5 + 3*x*x -4*x*x + 5*x^2 + p1; |
---|
95 | poly p3 = p2*p2; |
---|
96 | |
---|
97 | vector v; |
---|
98 | vector v1 = [p1, p2]; |
---|
99 | vector v2 = [p, p1,p2, p3] + p1*gen(5); |
---|
100 | |
---|
101 | ideal id; |
---|
102 | ideal id1 = p, p1, p2, p3; |
---|
103 | ideal id2 = p, p1+p2, p2+3, p1+p2+p3, p1, p2, p3; |
---|
104 | |
---|
105 | module mv; |
---|
106 | module mv1 = v, v1, v2; |
---|
107 | module mv2 = v1+v2, p1*v2, p2*(v1-v2), v, v1, v2; |
---|
108 | |
---|
109 | matrix m; |
---|
110 | matrix m1[2][3] = p, p1, p2, id1; |
---|
111 | matrix m2[4][4] = p, p1, p2, p3, id1, id2; |
---|
112 | |
---|
113 | |
---|
114 | ring r1 = 0, x(1..10), ls; |
---|
115 | |
---|
116 | number n; |
---|
117 | number n1 = 29734481274863241234589; |
---|
118 | number n2 = n1/(n1-6); |
---|
119 | |
---|
120 | poly p; |
---|
121 | poly p1 = n1*x(1) + n2*x(5)^4 + 7/3*x(3)^2*x(4)^5 - 6; |
---|
122 | poly p2 = x(2)*x(5)*x(6) - 2*x(7)^3*x(8)^4*x(9)^5 + 3*x(1)*x(5) -4*x(5)*x(9) + |
---|
123 | 5*x(9)^2 + p1; |
---|
124 | poly p3 = p2*p2; |
---|
125 | |
---|
126 | vector v; |
---|
127 | vector v1 = [p1, p2]; |
---|
128 | vector v2 = [p, p1,p2, p3] + p1*gen(5); |
---|
129 | |
---|
130 | ideal id; |
---|
131 | ideal id1 = p, p1, p2, p3; |
---|
132 | ideal id2 = p, p1+p2, p2+3, p1+p2+p3, p1, p2, p3; |
---|
133 | |
---|
134 | module mv; |
---|
135 | module mv1 = v, v1, v2; |
---|
136 | module mv2 = v1+v2, p1*v2, p2*(v1-v2), v, v1, v2; |
---|
137 | |
---|
138 | matrix m; |
---|
139 | matrix m1[2][3] = p, p1, p2, id1; |
---|
140 | matrix m2[4][4] = p, p1, p2, p3, id1, id2; |
---|
141 | |
---|
142 | ring r2=10,(x(1..6)),(lp(2),dp(4)); |
---|
143 | number n; |
---|
144 | number n1 = 29734481274863241234589; |
---|
145 | number n2 = n1/(n1-6); |
---|
146 | |
---|
147 | poly p; |
---|
148 | poly p1 = n1*x(1) + n2*x(5)^4 + 7/3*x(3)^2*x(4)^5 - 6; |
---|
149 | poly p2 = x(2)*x(5)*x(6) - 2*x(1)^3*x(2)^4*x(3)^5 + 3*x(1)*x(5) -4*x(5)*x(3) + |
---|
150 | 5*x(3)^2 + p1; |
---|
151 | poly p3 = p2*p2; |
---|
152 | |
---|
153 | vector v; |
---|
154 | vector v1 = [p1, p2]; |
---|
155 | vector v2 = [p, p1,p2, p3] + p1*gen(5); |
---|
156 | |
---|
157 | ideal id; |
---|
158 | ideal id1 = p, p1, p2, p3; |
---|
159 | ideal id2 = p, p1+p2, p2+3, p1+p2+p3, p1, p2, p3; |
---|
160 | |
---|
161 | module mv; |
---|
162 | module mv1 = v, v1, v2; |
---|
163 | module mv2 = v1+v2, p1*v2, p2*(v1-v2), v, v1, v2; |
---|
164 | |
---|
165 | matrix m; |
---|
166 | matrix m1[2][3] = p, p1, p2, id1; |
---|
167 | matrix m2[4][4] = p, p1, p2, p3, id1, id2; |
---|
168 | |
---|
169 | ring r3=(7,a, b, c),(x,y,z),dp; |
---|
170 | number n; |
---|
171 | number n1 = 2*a^2*b*c -3*b*c + -1*a + 4; |
---|
172 | number n2 = n1 / (n1 - 5c)^2; |
---|
173 | |
---|
174 | poly p; |
---|
175 | poly p1 = x + y + a; |
---|
176 | poly p2 = 3*a^2*xyz*n1 - 2*n2*x3y4z5 + 3xy -4yz + n1*z; |
---|
177 | poly p3 = p2*p1; |
---|
178 | |
---|
179 | vector v; |
---|
180 | vector v1 = [p1, p2]; |
---|
181 | vector v2 = [p, p1,p2, p3] + p1*gen(5); |
---|
182 | |
---|
183 | ideal id; |
---|
184 | ideal id1 = p, p1, p2, p3; |
---|
185 | ideal id2 = p, p1+p2, p2+3, p1+p2+p3, p1, p2, p3; |
---|
186 | |
---|
187 | module mv; |
---|
188 | module mv1 = v, v1, v2; |
---|
189 | module mv2 = v1+v2, p1*v2, p2*(v1-v2), v, v1, v2; |
---|
190 | |
---|
191 | matrix m; |
---|
192 | matrix m1[2][3] = p, p1, p2, id1; |
---|
193 | matrix m2[4][4] = p, p1, p2, p3, id1, id2; |
---|
194 | |
---|
195 | ring r4=(7,a),(x,y,z),dp; |
---|
196 | minpoly=a^2+a+3; |
---|
197 | poly p; |
---|
198 | number n; |
---|
199 | number n1 = 2*a^2 + 3; |
---|
200 | number n2 = n1 / (n1 -5)^3; |
---|
201 | |
---|
202 | poly p1 = x + y + a; |
---|
203 | poly p2 = 3*a^2*xyz*n1 - 2*n2*x3y4z5 + 3xy -4yz + n1*z; |
---|
204 | poly p3 = p2^4; |
---|
205 | |
---|
206 | vector v; |
---|
207 | vector v1 = [p1, p2]; |
---|
208 | vector v2 = [p, p1,p2, p3] + p1*gen(5); |
---|
209 | |
---|
210 | ideal id; |
---|
211 | ideal id1 = p, p1, p2, p3; |
---|
212 | ideal id2 = p, p1+p2, p2+3, p1+p2+p3, p1, p2, p3; |
---|
213 | |
---|
214 | module mv; |
---|
215 | module mv1 = v, v1, v2; |
---|
216 | module mv2 = v1+v2, p1*v2, p2*(v1-v2), v, v1, v2; |
---|
217 | |
---|
218 | matrix m; |
---|
219 | matrix m1[2][3] = p, p1, p2, id1; |
---|
220 | matrix m2[4][4] = p, p1, p2, p3, id1, id2; |
---|
221 | |
---|
222 | ring r5=(49,a),(x,y,z),dp; |
---|
223 | poly p; |
---|
224 | number n; |
---|
225 | number n1 = 2*a^2 + 3; |
---|
226 | number n2 = n1 / (n1 -5)^3; |
---|
227 | |
---|
228 | poly p1 = x + y + a; |
---|
229 | poly p2 = 3*a^2*xyz*n1 - 2*n2*x3y4z5 + 3xy -4yz + n1*z; |
---|
230 | poly p3 = p2^4; |
---|
231 | |
---|
232 | vector v; |
---|
233 | vector v1 = [p1, p2]; |
---|
234 | vector v2 = [p, p1,p2, p3] + p1*gen(5); |
---|
235 | |
---|
236 | ideal id; |
---|
237 | ideal id1 = p, p1, p2, p3; |
---|
238 | ideal id2 = p, p1+p2, p2+3, p1+p2+p3, p1, p2, p3; |
---|
239 | |
---|
240 | module mv; |
---|
241 | module mv1 = v, v1, v2; |
---|
242 | module mv2 = v1+v2, p1*v2, p2*(v1-v2), v, v1, v2; |
---|
243 | |
---|
244 | matrix m; |
---|
245 | matrix m1[2][3] = p, p1, p2, id1; |
---|
246 | matrix m2[4][4] = p, p1, p2, p3, id1, id2; |
---|
247 | |
---|
248 | ring r6 = 32003, (x,y,z), wp(1,2,3); |
---|
249 | poly p = (xyz - 2x3y4z5 + 3xy -4yz + 5z + 2)^2; |
---|
250 | |
---|
251 | ring r7 = 32003, (x,y,z), Ws(1,2,3); |
---|
252 | poly p = (xyz - 2x3y4z5 + 3xy -4yz + 5z + 2)^2; |
---|
253 | |
---|
254 | ring r8 = 32003, (a,b,c,x,y,z), (dp(3), wp(1,2,3)); |
---|
255 | poly p = (xyzabc - 2x3y4z5a2b4c3 + 3xya2 -4yzb3c2 + 5zc + 2)^2; |
---|
256 | |
---|
257 | ring r9 = 32003, (a,b,c,x,y,z), (a(1,2,3,4,5),Dp(3), ds(3)); |
---|
258 | poly p = (xyzabc - 2x3y4z5a2b4c3 + 3xya2 -4yzb3c2 + 5zc + 2)^2; |
---|
259 | |
---|
260 | ring r10 = 32003, (x,y,z), M(1, 0, 0, 0, 1, 0, 0, 0, 1); |
---|
261 | poly p = (xyz - 2x3y4z5 + 3xy -4yz + 5z + 2)^2; |
---|
262 | |
---|
263 | ring r11 = 0, x(1..10), (lp(2), M(1, 2, 3, 1, 1, 1, 1, 0, 0), ds(2), ws(1,2,3)); |
---|
264 | poly p1 = x(1) + 123399456085/(123399456085 + 1)*x(5)^4 + 7/3*x(3)^2*x(4)^5 - 6; |
---|
265 | poly p2 = x(2)*x(5)*x(6) - 2*x(7)^3*x(8)^4*x(9)^5 + 3*x(1)*x(5) -4*x(5)*x(9) + |
---|
266 | 5*x(9)^2 + p1; |
---|
267 | |
---|
268 | ring r12 = 0, (x,y,z), (C,lp); |
---|
269 | vector v1 = [x+y2,z3+xy]; |
---|
270 | vector v2 = [x,x,x]; |
---|
271 | |
---|
272 | ring r13 = 0, (x,y,z), (lp,c); |
---|
273 | vector v1 = [x+y2,z3+xy]; |
---|
274 | vector v2 = [x,x,x]; |
---|
275 | |
---|
276 | ring r14 = 0, (x,y,z), (c,lp); |
---|
277 | vector v1 = [x+y2,z3+xy]; |
---|
278 | vector v2 = [x,x,x]; |
---|
279 | |
---|