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