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

relay1987

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

	甲 ! {. V6 `# v: \- a+ K% C	乙 . X, L$ S4 _% J. u6 I  k3 K5 O( L: p	丙	丁
A & t% C; V# W% A8 p  }) D. K2 F	160   W( q7 t" g6 w0 L  j( E4 L% J	130	220 2 y0 |6 {7 Z/ h! z: M	170
B ' t$ {$ }9 g0 U" u" d9 o6 z" ]6 {	140	130	190	150 . O" O2 Y5 D% I
C 1 \" s! Z6 t3 O, ?	190	200	230 ! u9 d# `* q3 F	/

: f0 [+ S0 c) P, Q. b* H
# B( V9 d: k& A1 l* r$ w3 d
! Y, X  Z2 U5 m5 F
         单位：元/千吨
针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
% L. u/ g- M6 A& Z2 l" `1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

alanloveai

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

mrx

【1】求解：
6 d$ @  f  Q, `) D设水库A、B、C分别向甲、乙、丙、丁四个小区供水量为x11、x12、x13、x14、x21、x22、x23、x24、x31、x32、x33、x34，则有：
* t1 P  y# S$ p& w6 S8 D- w1 Emax 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
st
x34=0;
2 D1 @4 {( n/ |# c5 h1 _x11+x12+x13+x14<=50;
x21+x22+x23+x24<=60;
( s# }1 i6 z8 ^6 Bx31+x32+x33+x34<=50;
x11+x21+x31>=30;
) n% I4 D5 U. Px11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
6 j7 H# Y8 x) Cx13+x23+x33<=30;
x14+x24+x34>=10;
+ r3 [/ o6 N, G( I0 ^x14+x24+x34<=50;
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;
! v1 ~3 V' r1 V: b- ~x11+x12+x13+x14<=50;
x21+x22+x23+x24<=60;
3 R$ w& e: z0 X5 u* \; A" d$ J7 Zx31+x32+x33+x34<=50;
x11+x21+x31>=30;
+ S: k' W" d& w% v& s; mx11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
+ T! {  z6 b2 F9 R  ax14+x24+x34>=10;
) m  |& [  q8 ^3 fx14+x24+x34<=50;
# n/ x! |8 p6 K7 h* q: ]end
运行结果：
  Global optimal solution found.
: |- P- M: Q# T- a' W  Objective value:                              47600.00
  Total solver iterations:                             7


Variable           Value        Reduced Cost
3 K% D6 O( v( z" k1 M& N2 d) wX11        0.000000            30.00000
  s+ W6 ~1 N6 s) `X12        50.00000            0.000000
X13        0.000000            50.00000
X14        0.000000            20.00000
, W' R. i/ K7 m& s% H# FX21        0.000000            10.00000
/ }3 |9 q; A* F4 {$ U$ eX22        50.00000            0.000000
  G- M, a- y2 N; T: n9 nX23        0.000000            20.00000
  y4 f. G6 _! MX24        10.00000            0.000000
X31        40.00000            0.000000
2 x9 D+ z3 r- P# Z: V7 |X32        0.000000            10.00000
X33        10.00000            0.000000
8 g1 R  d0 k% H+ MX34        0.000000            240.0000

Row    Slack or Surplus      Dual Price
. U+ X8 n' }2 R. }' B1        47600.00            1.000000
/ I9 \3 E: e  H2 t! g# ^0 _3 W2        0.000000            320.0000
+ Q( Q0 t: {5 h  Y+ @9 j* G  {$ J3        0.000000            320.0000
3 T8 V: |! ^8 L' ]+ {% n4 n4        0.000000            260.0000
5        10.00000            0.000000
: H% w1 |% U# O2 g7 V6        40.00000            0.000000
7        30.00000            0.000000
% R* f4 P+ n+ {/ a! Z8        40.00000            0.000000
9        0.000000           -40.00000
10        20.00000            0.000000
11        0.000000           -20.00000
1 J% m1 x& t$ ^( @0 ?. C  q12        40.00000            0.000000

/ K# N+ r' j2 A- d8 W) G+ x【2】求解：
模型：
max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
# f$ U& b7 O8 w+ fst
x34=0;
x11+x12+x13+x14<=100;
6 F* z  B) P) ^7 A6 xx21+x22+x23+x24<=120;
) O0 i$ w, u* u4 V  y& P2 ax31+x32+x33+x34<=100;
' Z/ ~1 ]. l- @/ o8 U# v* X9 |x11+x21+x31>=30;
x11+x21+x31<=80;
6 |  n0 t/ a* ]9 j* Rx12+x22+x32>=70;
) v# q! R' c% t0 o% N; {3 e5 lx12+x22+x32<=140;
x13+x23+x33>=10;
& V% ~/ E( J4 Q. o: vx13+x23+x33<=30;
7 d9 q6 D: j; i! cx14+x24+x34>=10;
x14+x24+x34<=50;
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;
x21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
0 O0 A  t$ [/ G. o  O$ [! Gx11+x21+x31>=30;
) v  i% w( ?- O$ X! m$ |$ ?x11+x21+x31<=80;
x12+x22+x32>=70;
+ b2 x) X2 a7 v! g0 W- ^9 J4 @x12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
: n: b/ W* `$ T2 Z1 Gx14+x24+x34<=50;
end
运行结果：
  Global optimal solution found.
+ k* A1 {0 Z# u4 }& s& h+ W  Objective value:                              88700.00
  Total solver iterations:                             7
& \: P0 |+ W5 h, m" k  E

; z: I( Z, N% ~$ f  E4 ^6 O# \Variable           Value        Reduced Cost
) |! B5 a* T; y3 F0 }X11        0.000000            20.00000
X12        100.0000            0.000000
X13        0.000000            40.00000
; i; F/ ~/ g5 g4 x/ B7 A. m# gX14        0.000000            20.00000
6 r. o3 [! V5 G0 m4 UX21        30.00000            0.000000
X22        40.00000            0.000000
X23        0.000000            10.00000
7 d4 [9 q/ m* c7 C, T7 k! ?X24        50.00000            0.000000
: S& M. ?4 e9 BX31        50.00000            0.000000
X32        0.000000            20.00000
X33        30.00000            0.000000
7 e7 G: X" t& bX34        0.000000            250.0000
* D( [8 |( B  S( H/ G' n
Row    Slack or Surplus      Dual Price
: q& x6 h) x  q6 @5 X7 l) H2 r1        88700.00            1.000000
/ f3 d' C" U# t* [' j8 p2        0.000000            50.00000
9 \5 v" k3 r+ p$ X# X3        0.000000            50.00000
: _5 b5 X' T6 Q. s$ @7 \3 V7 q$ M5 b4        20.00000            0.000000
5        50.00000            0.000000
6        0.000000            260.0000
7        70.00000            0.000000
8        0.000000            270.0000
9        20.00000            0.000000
5 C6 Y. I) h- \* t9 o10        0.000000            220.0000
; q; W6 g) U2 O4 s  Z9 C11        40.00000            0.000000
12        0.000000            250.0000

relay1987

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