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

relay1987

" y; d7 S) n4 ]/ n
有A、B、C三个水库往甲、乙、丙、丁四个小区供应自来水，其中C水库由于条件限制不能给丁小区供水，这三个水库每天的供水量分别为50，60，50（单位：千吨），这四个小区每天的基本用水量为30，70，10，10（单位：千吨），在保证基本用水量的前提下，每个小区还需要一定量的额外用水，这四个小区每天额外需求量分别为50，70，20，40（单位：千吨），小区只需按照900元/千吨的价格缴纳费用，水库往这四个小区的引水管理费如下表，且除引水管理费之外，每千吨水还有450元的其他费用。
4 Y' V/ L1 E0 @) y! g% P! D  E2 r

	甲 # a+ h! c2 E' _8 q& s; f5 s0 R% u/ b	乙 0 v. q8 o. _2 l: {3 v2 W6 ?	丙 % R& [, b* v" i/ ^, j* d( c	丁 * |/ Z2 ~+ ^% n' |. y. H+ v, y
A 4 i: E: Z/ Q; y+ Z7 \	160	130 ( f6 [* X5 p: o- ?$ X  ^. J	220	170
B	140	130	190	150 ! v6 T% ~* z' b/ P
C 6 H' h# h) L6 J; E2 a	190	200	230 5 ]5 Q4 C* T3 m, w1 e	/

* P7 d5 |, |+ s7 K1 U


' n! ~( j9 N( x; @2 [6 V) q
         单位：元/千吨
针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

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

mrx

【1】求解：
1 _! T7 D' {( |* Z' \: P设水库A、B、C分别向甲、乙、丙、丁四个小区供水量为x11、x12、x13、x14、x21、x22、x23、x24、x31、x32、x33、x34，则有：
2 |5 Q. ^( k. ?$ m2 Lmax 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
6 R3 i7 E1 P6 X. ost
x34=0;
x11+x12+x13+x14<=50;
x21+x22+x23+x24<=60;
3 ?  e7 y( {7 b% l! L  |: ]x31+x32+x33+x34<=50;
x11+x21+x31>=30;
x11+x21+x31<=80;
( t9 X5 E) ^7 C# fx12+x22+x32>=70;
" A& C% C2 G' s2 g5 ex12+x22+x32<=140;
1 ^% W, y. R3 Z4 o/ P) {* fx13+x23+x33>=10;
x13+x23+x33<=30;
2 H4 a8 S2 c2 hx14+x24+x34>=10;
/ Z, ~$ m6 `  t6 d# }x14+x24+x34<=50;
end
% m' i: Q4 N: N  D0 oLINGO程序：
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;
5 }3 {" ^  o; N- {x31+x32+x33+x34<=50;
x11+x21+x31>=30;
& a. a3 ^" Z$ [0 v' d% `8 jx11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
# C- x) j' v" d) A7 }x14+x24+x34>=10;
# {( c- F0 r: q4 Mx14+x24+x34<=50;
end
) I- z& l# d" o4 h- P# U运行结果：
5 H( R/ l1 _" Z& Z/ V: Y) Y  Global optimal solution found.
2 R/ B, j2 G" F- _/ y; a  Objective value:                              47600.00
" j4 W1 \( A6 P+ {  Total solver iterations:                             7


Variable           Value        Reduced Cost
2 L& h) v7 E+ G3 a( `! g5 CX11        0.000000            30.00000
5 a- r% ^- c$ H  ^0 ]8 R: v7 IX12        50.00000            0.000000
X13        0.000000            50.00000
X14        0.000000            20.00000
X21        0.000000            10.00000
3 s/ e; Y( j- n( h/ _" Y) l0 }X22        50.00000            0.000000
X23        0.000000            20.00000
0 R: H% q3 [9 u/ K( FX24        10.00000            0.000000
  S: @5 Q8 R2 Y6 ?" p% UX31        40.00000            0.000000
X32        0.000000            10.00000
3 A! t2 c' y! j$ @' Y1 pX33        10.00000            0.000000
X34        0.000000            240.0000
: D) Z) F0 n' `+ Q% v
Row    Slack or Surplus      Dual Price
  d* B; G- F5 Q' e3 `0 K7 G5 d0 I# P1        47600.00            1.000000
2        0.000000            320.0000
3        0.000000            320.0000
( l# P2 X6 b# ~8 O% d4        0.000000            260.0000
* n" Q. u0 D* X8 |- K5        10.00000            0.000000
6        40.00000            0.000000
7        30.00000            0.000000
8        40.00000            0.000000
9        0.000000           -40.00000
10        20.00000            0.000000
% ]& M; N' V& |/ N11        0.000000           -20.00000
12        40.00000            0.000000
4 b7 Z7 q3 `: d# o. w
6 d4 G8 H) B% \4 R9 n【2】求解：
模型：
max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
st
x34=0;
! c4 T( y0 g; O& }2 e. qx11+x12+x13+x14<=100;
" {% s' T: i* N) f3 B. N" m) Ix21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
4 f, O  C6 l* d1 Q0 l# fx11+x21+x31>=30;
* _. f- W' L6 ~6 M% x) Gx11+x21+x31<=80;
x12+x22+x32>=70;
& o/ }" g( d" P+ [% rx12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
x14+x24+x34<=50;
6 A3 V5 N4 }; u% F+ `7 V  W4 y- {end
LINGO程序：
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<=100;
* O' x8 b# u3 D% }" v5 R2 f% v' }6 ix21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
x11+x21+x31>=30;
x11+x21+x31<=80;
7 Z- w+ i/ X* n! _$ Cx12+x22+x32>=70;
. v) p9 v* n- w8 k4 y9 A4 c2 K& tx12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
/ n- q5 Z% c2 k* x" Sx14+x24+x34>=10;
) A1 E" R& [7 Qx14+x24+x34<=50;
( s& W- y' t1 v! Zend
( `9 `! s8 c/ s+ z  H5 ~运行结果：
7 K. ^1 A: @- n, Q0 e  Global optimal solution found.
# }& {* c- O$ W+ J( f4 B  Objective value:                              88700.00
# e) z1 s# q. J) F5 l0 o  Total solver iterations:                             7
1 x% r6 j5 h8 z! u' O
" T$ N( O; H# S( W0 q" e) L
Variable           Value        Reduced Cost
X11        0.000000            20.00000
% Y# W8 T" K' B0 o; c8 sX12        100.0000            0.000000
X13        0.000000            40.00000
: H7 F3 b. {" rX14        0.000000            20.00000
X21        30.00000            0.000000
( j1 b5 _1 l3 }X22        40.00000            0.000000
X23        0.000000            10.00000
, w' l* {( {4 R, `) mX24        50.00000            0.000000
X31        50.00000            0.000000
+ J* g$ d# {( H, x/ OX32        0.000000            20.00000
X33        30.00000            0.000000
7 i8 O) E0 h$ }9 z( U  w  ~- u! IX34        0.000000            250.0000

Row    Slack or Surplus      Dual Price
1        88700.00            1.000000
2 Z+ |) @5 U  e2        0.000000            50.00000
0 U0 D/ u: s2 b, C5 R2 a- l$ @3 z/ f2 l3        0.000000            50.00000
4        20.00000            0.000000
5        50.00000            0.000000
6        0.000000            260.0000
2 d. v: F$ C! {% |+ A! t& t7        70.00000            0.000000
2 P9 `7 u& M0 o9 l" t, {4 ]/ |$ j+ G8        0.000000            270.0000
3 u8 |/ j6 L- k1 h1 D3 z% r9        20.00000            0.000000
10        0.000000            220.0000
; `5 L/ m( P- T1 q8 W$ `11        40.00000            0.000000
4 r% i1 w5 T7 ]12        0.000000            250.0000

alanloveai

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