急求该题答案，我试了10次也不会运行程序

relay1987

' c- k2 g& s& q有A、B、C三个水库往甲、乙、丙、丁四个小区供应自来水，其中C水库由于条件限制不能给丁小区供水，这三个水库每天的供水量分别为50，60，50（单位：千吨），这四个小区每天的基本用水量为30，70，10，10（单位：千吨），在保证基本用水量的前提下，每个小区还需要一定量的额外用水，这四个小区每天额外需求量分别为50，70，20，40（单位：千吨），小区只需按照900元/千吨的价格缴纳费用，水库往这四个小区的引水管理费如下表，且除引水管理费之外，每千吨水还有450元的其他费用。

0 G* b" P& v+ S$ U( I2 P# s	甲	乙 5 J' {/ l' s; A- w8 ^' R* ?	丙 1 s# E/ z5 J* r0 ?	丁 " f& K7 W8 ?" Z7 w& U* j2 y
A	160 0 K0 F$ d: q0 S1 ^	130 5 \3 Z! v& ]/ x# E% O& w$ i	220	170 : l0 k/ z" @9 J- V7 q8 s
B % Q% B4 f! ]) U$ W  F	140	130 : V, _2 j) W! P) w* Z. e. M	190 # f+ @* I8 H+ r0 R* K' z	150 & ]( T/ w8 ]0 m$ t7 G
C , a& S" R4 h5 B* n$ d" J8 x/ J7 e	190	200	230	/ ' Y7 i$ R0 v) z* H4 H# K3 E

9 Q6 e8 u' f8 p( X
) K* M$ _" f: w5 a

8 [8 H6 a" i$ K" h+ z3 T6 O         单位：元/千吨
针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

希望大家快些帮帮我，我参加学校比赛，老师病了，我不想打扰他，希望大家帮帮我，后天比赛，这些软件我学的不好，所以不怎么会用！今天看书写了个程序，可是不会运行，请高手们救救我吧，写好后发我的邮箱就行了relay1987@163.com，谢谢了！急急急，期待您的妙手回春！

mrx

【1】求解：
设水库A、B、C分别向甲、乙、丙、丁四个小区供水量为x11、x12、x13、x14、x21、x22、x23、x24、x31、x32、x33、x34，则有：
max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
st
x34=0;
# _3 B1 C+ E4 o- z; S2 }x11+x12+x13+x14<=50;
! o5 q) s+ m5 \# \. E+ x% \x21+x22+x23+x24<=60;
x31+x32+x33+x34<=50;
* D! H/ w7 x$ fx11+x21+x31>=30;
x11+x21+x31<=80;
7 l8 F2 `$ s/ m8 |8 Jx12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
  J. {+ T) U4 s8 ^% O$ fx13+x23+x33<=30;
x14+x24+x34>=10;
$ ?' `, V: u( I* Dx14+x24+x34<=50;
. Q  n$ R* Q7 w1 F* lend
7 o, S8 s) {9 N: NLINGO程序：
max=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
x11+x12+x13+x14<=50;
x21+x22+x23+x24<=60;
6 X$ O( _5 N. B3 }: l0 H, L) qx31+x32+x33+x34<=50;
x11+x21+x31>=30;
x11+x21+x31<=80;
" M/ r/ C: i+ x: A! y9 mx12+x22+x32>=70;
: m6 p4 W0 G3 L5 |4 sx12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
x14+x24+x34<=50;
end
% k7 e$ W4 o' {- g运行结果：
  Global optimal solution found.
  Objective value:                              47600.00
  Total solver iterations:                             7
! @6 c9 z6 ^% p. Z

% e( b- l9 R) k" b6 a" h9 K, }Variable           Value        Reduced Cost
X11        0.000000            30.00000
7 _- j/ G( ]" jX12        50.00000            0.000000
1 L$ f3 b/ M( @3 _2 ^8 B) X) \X13        0.000000            50.00000
X14        0.000000            20.00000
X21        0.000000            10.00000
X22        50.00000            0.000000
3 X( j, i6 F& A. c6 o& qX23        0.000000            20.00000
% Y6 n1 |7 v! Z3 zX24        10.00000            0.000000
2 N8 }. y' u, [$ G, V/ LX31        40.00000            0.000000
2 C% l5 [( I6 @4 ]0 Y% aX32        0.000000            10.00000
1 R6 D; ~+ o9 [) Z# t5 TX33        10.00000            0.000000
X34        0.000000            240.0000

7 v2 c( Y$ Z, O) F1 wRow    Slack or Surplus      Dual Price
1        47600.00            1.000000
2        0.000000            320.0000
3        0.000000            320.0000
4        0.000000            260.0000
( h$ I0 s- K. U: Z4 @9 ]5        10.00000            0.000000
8 h. z! s* a2 u6        40.00000            0.000000
. c$ R8 \: j$ O( ]7 i$ H' H& m  i; ?: z! n7        30.00000            0.000000
8        40.00000            0.000000
( T$ h1 |: J% {9        0.000000           -40.00000
& `3 h! A% K  G9 q+ z+ f1 O  W10        20.00000            0.000000
11        0.000000           -20.00000
12        40.00000            0.000000

* v" A2 L& _: ~# R/ {【2】求解：
" s0 v: ]6 T- S$ y  O& @模型：
: S2 J* X. H# Ymax 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
# k& i7 g8 D. S; ust
x34=0;
) t3 W6 f, _1 d$ B# {8 @  gx11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
: ]/ O4 B4 t4 U7 Ix31+x32+x33+x34<=100;
x11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
' ], O2 x* c0 E( r) O) y/ O8 nx12+x22+x32<=140;
x13+x23+x33>=10;
* q6 Z! _' W+ I' U- }x13+x23+x33<=30;
x14+x24+x34>=10;
, s4 n6 v+ H) y7 ]5 `% fx14+x24+x34<=50;
end
# A( H5 Y3 H0 v4 z) vLINGO程序：
( q/ a3 M+ `# a7 b& O5 Pmax=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
, b6 t: [7 M) T+ W9 K. n5 Jx11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
x11+x21+x31>=30;
x11+x21+x31<=80;
: m- T6 ^: I" n3 U" x6 |x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
4 r  e0 c: @/ x  I2 Fx13+x23+x33<=30;
x14+x24+x34>=10;
x14+x24+x34<=50;
# u1 ~' K7 b& F& o' M. Rend
6 C; h- W5 s$ O! H运行结果：
- H3 m% o/ }1 Q# s8 O  Global optimal solution found.
& K6 W7 J" q! N  Objective value:                              88700.00
  Total solver iterations:                             7


0 i) h8 e& d- j* d# @Variable           Value        Reduced Cost
X11        0.000000            20.00000
X12        100.0000            0.000000
& C$ ^' t  R% |3 [8 E, f& rX13        0.000000            40.00000
, O6 S0 G, c- _2 m9 X$ {4 gX14        0.000000            20.00000
X21        30.00000            0.000000
X22        40.00000            0.000000
% M4 O0 ]: i5 z1 b$ j+ ]X23        0.000000            10.00000
3 l/ T: O9 n; _0 M% q2 f9 c! \4 hX24        50.00000            0.000000
X31        50.00000            0.000000
X32        0.000000            20.00000
X33        30.00000            0.000000
X34        0.000000            250.0000

: @+ M/ C+ G" |4 T+ r4 d: WRow    Slack or Surplus      Dual Price
% A9 b6 t& a8 ], h5 P& W1        88700.00            1.000000
2        0.000000            50.00000
3        0.000000            50.00000
2 G* L1 ~- Q( n  P* q4        20.00000            0.000000
( M( u: [% {! k% l5        50.00000            0.000000
6 T0 F0 s5 n) ^# U2 ]/ ?- y6        0.000000            260.0000
) ]; F( r. g& M7        70.00000            0.000000
3 o8 m: J. d% @* ^' {( w/ A2 f6 S5 F8        0.000000            270.0000
9        20.00000            0.000000
10        0.000000            220.0000
11        40.00000            0.000000
$ H) H3 P9 {% ?  s12        0.000000            250.0000

alanloveai

强············