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