Maple calculations: Difference between revisions

From Polymath Wiki
Jump to navigationJump to search
New page: : with(simplex); : A3 := 8*a + 12*b + 6*c + d; : A4 := 16*a+32*b+24*c+8*d+e; : A5 := 32*a+80*b+80*c+40*d+10*e+f; : A6 := 64*a+192*b+240*c+160*d+60*e+12*f+g; : A7 := 128*a+448*b+672*c+560*d...
 
No edit summary
Line 1: Line 1:
: with(simplex);
The computations below were done in Maple 12.
: A3 := 8*a + 12*b + 6*c + d;
 
: A4 := 16*a+32*b+24*c+8*d+e;
:with(Optimization);
: A5 := 32*a+80*b+80*c+40*d+10*e+f;
:A3 := 8*a + 12*b + 6*c + d;
: A6 := 64*a+192*b+240*c+160*d+60*e+12*f+g;
:A4 := 16*a+32*b+24*c+8*d+e;
: A7 := 128*a+448*b+672*c+560*d+280*e+84*f+14*g+h;
:A5 := 32*a+80*b+80*c+40*d+10*e+f;
: X1 := [2*a+b <= 2];
:A6 := 64*a+192*b+240*c+160*d+60*e+12*f+g;
: X2 := [4*a+2*b+c <= 4, 2*a+c <= 2, 2*b+c <= 2, c <= 1];
:A7 := 128*a+448*b+672*c+560*d+280*e+84*f+14*g+h;
: X3 := [op(X2), op(subs([a=b,b=c,c=d],X2)), 8*a+6*b+6*c+2*d<=11, 4*a+4*b+3*c+d<=6, 7*a+3*b+3*c+d <= 7, 4*a+6*c+2*d<=7, 6*a+3*c+d<=5, 2*a+d <= 2];
:X1 := {2*a+b <= 2};
: X4 := [op(X3), op(subs([a=b,b=c,c=d,d=e],X3)), op(subs([b=c,c=e],X2)),A4 <= 43];
:X2 := {4*a+2*b+c <= 4, 2*a+c <= 2, 2*b+c <= 2, c <= 1};
: X5 := [op(X4), op(subs([a=b,b=c,c=d,d=e,e=f],X4)), 2*a+f <= 2, A5 <= 124];
:X3 := {op(X2), op(subs([a=b,b=c,c=d],X2)), 8*a+6*b+6*c+2*d<=11, 4*a+4*b+3*c+d<=6, 7*a+3*b+3*c+d <= 7, 4*a+6*c+2*d<=7, 6*a+3*c+d<=5, 2*a+d <= 2};
: X6 := [op(X5), op(subs([a=b,b=c,c=d,d=e,e=f,f=g],X5)), op(subs([b=d,c=g],X2)), op(subs([b=c,c=e,d=g],X3))];
:X4 := {op(X3), op(subs([a=b,b=c,c=d,d=e],X3)), op(subs([b=c,c=e],X2)),A4 <= 43};
: evalf(subs(maximize(A6,X6,NONNEGATIVE),A6));
:X5 := {op(X4), op(subs([a=b,b=c,c=d,d=e,e=f],X4)), 2*a+f <= 2, A5 <= 124};
364.6984127
:X6 := {op(X5), op(subs([a=b,b=c,c=d,d=e,e=f,f=g],X5)), op(subs([b=d,c=g],X2)), op(subs([b=c,c=e,d=g],X3))};
: evalf(subs(maximize(A6,X6,NONNEGATIVE),[64*a,192*b,240*c,160*d,60*e,12*f,g]));
:Y6 := subs([a=a/64,b=b/192,c=c/240,d=d/160,e=e/60,f=f/12],X6);
[29.81627297, 98.26771654, 128.5039370, 77.27034121, 21.73228346, 8.409448819, 0.]
:LPSolve(a+b+c+d+e+f+g,Y6,assume=nonnegint,depthlimit=100,maximize);
: evalf(subs(maximize(A6,[op(X6),g=1],NONNEGATIVE),A6));
[364, [a = 31, b = 94, c = 126, d = 83, e = 23, f = 7, g = 0]]
355
:LPSolve(a+b+c+d+e+f+g,{op(Y6),g=1},assume=nonnegint,depthlimit=100,maximize);
: X6 := [op(X6), A6 <= 364];
[355, [a = 24, b = 96, c = 120, d = 80, e = 30, f = 4, g = 1]
: X7 := [op(X6), op(subs([a=b,b=c,c=d,d=e,e=f,f=g,g=h],X6)), 2*a+h <= 2];
:X6 := {op(X6), A6 <= 364};
: evalf(subs(maximize(A7,X7,NONNEGATIVE),A7));
:X7 := {op(X6), op(subs([a=b,b=c,c=d,d=e,e=f,f=g,g=h],X6)), 2*a+h <= 2};
1092
:Y7 := subs([a=a/128,b=b/448,c=c/672,d=d/560,e=e/280,f=f/84,g=g/14],X7);
:LPSolve(a+b+c+d+e+f+g+h,Y7,assume=nonnegint,depthlimit=100,maximize);
 
[1092, [a = 63, b = 224, c = 336, d = 280, e = 140, f = 42, g = 7, h = 0]]
 
:LPSolve(a+b+c+d+e+f+g+h,{op(Y7),h=1},assume=nonnegint,depthlimit=100,maximize);
 
[1077, [a = 64, b = 224, c = 336, d = 280, e = 140, f = 28, g = 4, h = 1]]

Revision as of 17:06, 11 March 2009

The computations below were done in Maple 12.

with(Optimization);
A3 := 8*a + 12*b + 6*c + d;
A4 := 16*a+32*b+24*c+8*d+e;
A5 := 32*a+80*b+80*c+40*d+10*e+f;
A6 := 64*a+192*b+240*c+160*d+60*e+12*f+g;
A7 := 128*a+448*b+672*c+560*d+280*e+84*f+14*g+h;
X1 := {2*a+b <= 2};
X2 := {4*a+2*b+c <= 4, 2*a+c <= 2, 2*b+c <= 2, c <= 1};
X3 := {op(X2), op(subs([a=b,b=c,c=d],X2)), 8*a+6*b+6*c+2*d<=11, 4*a+4*b+3*c+d<=6, 7*a+3*b+3*c+d <= 7, 4*a+6*c+2*d<=7, 6*a+3*c+d<=5, 2*a+d <= 2};
X4 := {op(X3), op(subs([a=b,b=c,c=d,d=e],X3)), op(subs([b=c,c=e],X2)),A4 <= 43};
X5 := {op(X4), op(subs([a=b,b=c,c=d,d=e,e=f],X4)), 2*a+f <= 2, A5 <= 124};
X6 := {op(X5), op(subs([a=b,b=c,c=d,d=e,e=f,f=g],X5)), op(subs([b=d,c=g],X2)), op(subs([b=c,c=e,d=g],X3))};
Y6 := subs([a=a/64,b=b/192,c=c/240,d=d/160,e=e/60,f=f/12],X6);
LPSolve(a+b+c+d+e+f+g,Y6,assume=nonnegint,depthlimit=100,maximize);
[364, [a = 31, b = 94, c = 126, d = 83, e = 23, f = 7, g = 0]]
LPSolve(a+b+c+d+e+f+g,{op(Y6),g=1},assume=nonnegint,depthlimit=100,maximize);
[355, [a = 24, b = 96, c = 120, d = 80, e = 30, f = 4, g = 1]
X6 := {op(X6), A6 <= 364};
X7 := {op(X6), op(subs([a=b,b=c,c=d,d=e,e=f,f=g,g=h],X6)), 2*a+h <= 2};
Y7 := subs([a=a/128,b=b/448,c=c/672,d=d/560,e=e/280,f=f/84,g=g/14],X7);
LPSolve(a+b+c+d+e+f+g+h,Y7,assume=nonnegint,depthlimit=100,maximize);
[1092, [a = 63, b = 224, c = 336, d = 280, e = 140, f = 42, g = 7, h = 0]]
LPSolve(a+b+c+d+e+f+g+h,{op(Y7),h=1},assume=nonnegint,depthlimit=100,maximize);
[1077, [a = 64, b = 224, c = 336, d = 280, e = 140, f = 28, g = 4, h = 1]]