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