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
|
# {i,j,k | k=0; 1<=i<=10; 1<=j<=10} |
# {i,j,k | k=0; 28<=i<=31; 1<=j<=10}
D 2
5 5
0 0 0 1 0
1 1 0 0 -1
1 -1 0 0 10
1 0 1 0 -1
1 0 -1 0 10
5 5
0 0 0 1 0
1 1 0 0 -28
1 -1 0 0 31
1 0 1 0 -1
1 0 -1 0 10
# {i,j,k | k=0; i+j>=25; i<=25; j<=25} |
# {i,j,k | k=0; 2<=j<=9; 30<=i<=32}
D 2
4 5
0 0 0 1 0
1 1 1 0 -25
1 -1 0 0 25
1 0 -1 0 25
5 5
0 0 0 1 0
1 0 1 0 -2
1 0 -1 0 9
1 1 0 0 -30
1 -1 0 0 32
F 2
DomSimplify[readDom["{i,j,k | k=0; 1<=i<=10; 1<=j<=10} |
{i,j,k | k=0; 28<=i<=31; 1<=j<=10}"],
readDom["{i,j,k | k=0; i+j>=25; i<=25; j<=25} |
{i,j,k | k=0; 2<=j<=9; 30<=i<=32}"]]
(*answer =
domain[3, {i, j, k}, {pol[2, 4, 0, 2,
{{1, 1, 0, 0, -28}, {1, -1, 0, 0, 31}},
{{0, 0, 1, 0, 0}, {0, 0, 0, 1, 0}, {1, 28, 0, 0, 1}, {1, 31, 0, 0, 1}}]}]
{i,j,k|i-28>=0;-i+31>=0}
*)
|
