[0e1846] | 1 | proc mptest (lv) |
---|
| 2 | { |
---|
| 3 | print(lv); |
---|
[feaddd] | 4 | write("MPfile:w /tmp/t.mp", lv); |
---|
| 5 | return(lv == read("MPfile: /tmp/t.mp")); |
---|
[0e1846] | 6 | }; |
---|
| 7 | |
---|
| 8 | proc mpt(lv) |
---|
| 9 | { |
---|
| 10 | lv; |
---|
[feaddd] | 11 | write("MPfile:w /tmp/t.mp", lv); |
---|
| 12 | def r = read("MPfile: /tmp/t.mp"); |
---|
[0e1846] | 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 | |
---|
[feaddd] | 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; |
---|
[0e1846] | 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 | |
---|