眼科病床的合理安排

xlfd1988

2009高教社杯全国大学生数学建模竞赛

承  诺  书

9 N5 J; z, [; t; \$ X4 _, H我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
' B* A% F0 q. r0 X  }6 T: H' _0 X我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。

/ Y. ]. X6 P! A2 k8 g8 o- |我们参赛选择的题号是（从A/B/C/D中选择一项填写）：          B            
8 N0 p& K1 {" a        我们的参赛报名号为（如果赛区设置报名号的话）：    B甲2425                    
& }; f1 n- z7 q* P- b$ ^, x+ s所属学校（请填写完整的全名）：                  山东科技大学               
& o( \. f/ c' X参赛队员 (打印并签名) ：1.                    阎  冬                          
4 C$ }2 Q& E/ q; |$ u                       2.                    冯宗瑞                       
                       3.                    杨  超                     
指导教师或指导教师组负责人  (打印并签名)：          郭秀荣   王娟            
! U" `  |" ^9 J" S6 l& C, D
6 O' X3 D: |; i1 Z                                         日期：  2009  年  9  月 14 日
/ l- ]- ^+ O5 t+ y/ i0 R4 q

- ^2 s% X$ s( f$ X/ |; j1 v


赛区评阅编号（由赛区组委会评阅前进行编号）：
7 {1 I7 }( J1 k  u

, H% E. H& R$ K8 o8 E5 T2008高教社杯全国大学生数学建模竞赛

; I( d6 O. R3 M4 d编 号 专 用 页


- J- H! V# |" ]& e7 _' s
赛区评阅编号（由赛区组委会评阅前进行编号）：

7 u5 ?# {2 C7 [* H, k+ ~6 t

1 ?+ j  [0 H# A# r" x- U* R赛区评阅记录（可供赛区评阅时使用）：
5 k  x( t/ a/ p" `# {
评
( J9 o6 W2 ~9 _3 ^6 t0 A; }$ {阅
0 i4 E# |0 ]& J8 X人
                                                                               


评
分
) j3 U5 M; I0 T) b* ]% y" ]                                                                               

0 _. H" ^( E5 u备
$ c0 ?, P3 o0 n, k注
                                                                               
) [! J7 z( e% N5 j1 ]



7 N5 F/ ^6 q- l8 ~$ t/ {5 j/ \全国统一编号（由赛区组委会送交全国前编号）：




. C2 |& n1 F" w8 l/ |, y9 M
全国评阅编号（由全国组委会评阅前进行编号）：
  ]+ A) _$ c% N7 d眼科病床的合理安排的研究

- M0 C- t7 d+ @6 ^# O, {, R# j# s1 \【摘要】：合理分析床位利用情况对于提高医院经济效益，改善病房管理，挖掘潜力增强服务能力等有十分重要的意义。我们通过大量的分析给出了影响医院床位安排合理性的因素，并利用层次分析法给出了各个因素对于床位安排合理度的权值。利用excel对附表中的数据进行处理给出每个因素对某天的病床安排对某位病人的合理度评分，结合权值构造有权综合合理度模型、线性规划模型，并利用lingo编程求解，可以根据任一天的预出院人数给出该天哪些病人可以被安排住院。
    对于问题三，我们首先利用Excel处理数据，得出自8月1号至9月10号的每种类型的病人等待住院的时间 、已住院人数 以及等待住院人数 。利用matlab画出 与 以及 与 的散点图和拟合曲线，构造多元完全二次多项式回归模型，并给出交互式图形及回归曲线方程的估计参数。借此可预测任何时刻某类门诊病人的等待入院时间区间。
, Q! \$ j1 ]9 z考虑手术对病床安排合理度的影响因素，将问题二的模型进行改进，构造了新的线性规划模型。并利用lingo编程求解。同时对影响住院时间的因素进行分析。
对于问题五，先将每类型病人的入院到出院作为一个新的系统（这个系统的选取很关键），由于该系统满足排队论的条件，所以可利用马尔科夫排队模型给出每类病人的平均系统内逗留时间，再结合每类病人的平均等待入院时间，构造了非线性规划和线性规划两个模型，并编程运行求解，最终给出最优的床位分配比例。
, A5 L" ~2 x: f) Y* X! ~, x

) U. x; h8 ~4 H, `9 Q/ Q

: v7 Y) G) n6 t8 y6 R6 h4 H9 A
- p; }% {5 b/ ]! \! E
% q+ a7 m2 b1 h! k* E) u1 H

【关键词】：层次分析法，有权综合合理度模型，线性规划模型，回归模型，排队模型
非线性规划模型，lingo，matlab


5 N. P6 q4 F9 U. T5 r, j


# G2 V* o5 G: s8 j. m0 n
9 q6 i1 n; W  I0 o4 r
1 ]2 q- m% S# m! U7 H( x, T  D



4 T) f( y8 d, \+ A4 L2 M# s

' R2 x+ i: \' b! q6 [. m1.问题重述
. f' i& f1 _- I' O（略）
2.问题假设
1.假设该医院床位数固定，不加临时床位;
& k0 p, {$ K# p
2.假设医院手术条件充分，不考虑手术条件限制;
! u& n- D* U  x7 o; y: |5 m
1 O- c$ m4 c9 O1 I1 h6 }" h# B' D3.假设除外伤外其余眼科疾病不考虑急症;

4.假设病人的门诊及出院相互独立.
6 o$ g9 g! N7 K2 O9 ~1 Z& A! |3．符号说明
$ X0 N; m4 d( O% R- i
表示病人类型在第 日的病床安排对第 个人的合理性评分


表示等待住院时间在第 日的病床安排对第 个人的合理性评分

表示病床周转率在第 日的病床安排对第 个人的合理性评分

3 I- C5 I/ K7 }8 \- ~
表示第 日病床安排对第 个待住院者的合理度分值

' |4 P  @% w9 [7 O$ |2 }
, R& m# G$ x7 ?+ {表示视网膜病人已住院人数

表示视网膜病人待住院人数

" \& ?2 |/ R+ x$ ~: c( ~表示视网膜病人等待住院时间长度(以天为单位)

表示白内障病人已住院人数
# M( q- y& A; f, K, \, ^, E
表示白内障病人待住院人数

表示白内障病人等待住院时间长度(以天为单位)
8 a& c. ^6 i' ^' b6 s8 ~. |$ q
1 g3 \& M; H" i8 T$ ^9 j第 天预出院人数


第 类病人在新系统内的平均逗留时间
, B3 |5 ?, V4 u& ]4 Y3 v

第 类病人的日平均住院率
- S! g, O! Y9 c
+ B" M. n: B1 ?# }
系统对第 类病人的日平均服务率
' q% I; z1 o0 `$ h5 H; ]

. m. J% J0 f3 n% I9 C  E) S第 类病人的平均等待时间


各类病床分配比例


/ G0 w8 h0 V3 K3 J4．问题的分析及模型的建立及求解
* I! @  U( z5 y- W4.1问题一
* c* P; _0 a3 t2 n9 d病床安排评价指标是医院信息管理最活跃的要素之一。合理分析床位利用情况对于提高医院经济效益，改善病房管理，挖掘潜力增强服务能力等有十分重要的意义。由于医院病床的安排评价指标具有多系统和多层次的性质，因此评价医院的病床安排工作不能从单一的指标入手，要使用多指标综合评价，才能比较全面地反映实际病床工作的综合情况。只有合理配置医院床位，适度规模，科学发展，才能有效利用病床资源发挥其作用。
( U( d/ N2 l! L( \, Y4.1.1问题一中评价指标的选择
我们从问诊时间、病人住院时间及病人病情是否紧急三方面出发，考虑到病人病情对住院的需求度；又兼顾出院人数、住院人数、床位日平均占用数等方面，考虑到医院的床位有效利用率；另外根据本题医院的实际情况，排除了一些虽然存在但对本题所给医院的病床安排评价的影响不是很大或表现不是很明显的影响指标。其中，由于病人的手术费远高于病人住院期间的观察费，因此在考虑提高医院效益时认为观察期越小，床位周转率越大，病人入院率越高，手术收益越大。并考虑到模型的可靠性，即尽量满足各指标独立性的要求，排除了这些因素之间的横向与纵向联系，且符合实际情况。据此得出评价医院床位安排合理度的三个主要指标为：病人类型、等待住院时间、床位有效利用率。
4 ], _% Z3 Y7 l4.1.2 “层次分析”模型
7 L% @/ [) L5 h& Z利用“层次分析法”模型对病床利用情况进行评价，为医院管理提供决策依据。
将选取的病人类型、等待住院时间、床位有效利用率三个指标分别记作 , , ,得到成对比较矩阵
4 Z! D! T3 o5 c5 |

求行和得列向量
，
; f4 e! x+ ]6 f) S归一化得权向量
" F: l4 ^8 I) H# W 。
* |+ ~7 u3 ~: n; l计算成对比较矩阵的最大特征值
3.0554
计算一致性指标
0.0277（ 为因素个数，此题中 为3）

, N1 g0 @* G' `" @0 A% \为了确定成对比较矩阵的不一致程度的容许范围，Saaty引入了随机一致性指标 。Saaty给出的 数值如下：
: q" a4 ~+ Q5 k* _# q" T  ln        1        2        3        4        5        6        7        8        9        10        11
RI        0        0        0.58        0.90        1.12        1.24        1.32        1.41        1.45        1.49        1.51
3 {- ]& w! J+ b6 m2 _1 P& h! @! _$ ]当 >2时，称 为一致性比率，求出 0.0477<0.1（其中0.1为选取的主观信度）
4 G, c. _* T' W8 o( C由此可知，一致性检验通过，由成对比较矩阵计算出的权系数向量可以作为计算有权综合合理度的权系数。
4.2问题二模型的建立与求解

% g5 t% L4 G3 H4.2.1问题二的分析及数据处理
    医院床位安排对病人的合理度可以简要的定义为：病人对医院床位安排的满意程度。
4 K/ v) G( A/ x; p6 x# _' L$ G9 g, ^我们对附录中的数据利用Excel进行统计、汇总及分类运算处理，结合病人实际情况分析得出：第 因素在第 日病床安排中对第 人的合理度评分值，记为 。其中 =1表示病人类型； =2表示等待住院时间； =3表示病床周转率
$ O2 v. w$ P5 c    下面是根据实际情况考虑到各方面因素的影响，对 类指标合理度评分的选取做一详细的介绍。
# m6 S' u* N8 V& S通过分析疾病类型对病人生活的影响程度及病情的紧急程度，给出病床安排对4类病人(病情类型)的合理度的评分为
  
1 E; S* T4 O3 D! Z# \类似的可以求出其余两指标在第 日病床安排时对第 人的合理度评分
& g/ D/ w. s) i ×10
其中 为截止到第 日最长的病人等待入院时间
      
由问题一中根据层次分析模型得出的三指标在医院床位安排合理度中的权重表示为 结合医院床位安排对（特定时间特定病情的）病人合理度的评分，建立有权综合合理度评价模型。
4.2.2有权综合合理度评价模型
利用有权综合合理度评价指标得出第 日病床安排对第 个待住院者的合理度分值

问题二要使某日病床安排合理，只需该天病床安排的综合合理度评分（即第 日所有等待住院的病人对医院床位安排满意程度的和）取得最大值即可。
即得第 日病床的安排对全部待住院病人合理度评分的函数：
4.2.3 0-1规划模型
令
  
  y% Z: I! H% B则由上述分析得到如下的0-1规划模型：
目标函数:
8 x: G3 d2 c# H7 }$ Dmax
约束条件:
% H- E& `* P3 O. As.t.
其中 表示到第 日前待住院病人数;  表示第 日预出院病人数;
根据此模型可求得 ，即可根据第二天预出院病人数来确定第二天安排哪些病人住院较合理。
4.2.4问题二模型的应用评价
1 z0 c- \! j9 x& X; t3 y) P9 W7 L, \( r对于上述的模型，利用Lingo或Matlab软件编程求解0-1线性规划模型，并以9月1日分配住院为例，给出病床的合理分配，如下：
, i( I1 N: L- v$ o3 r门诊时间        8-19        8-19        8-19        8-19        8-20        8-20
: v- Z* j+ Z0 F/ i8 l类型        视网膜疾病        视网膜疾病        视网膜疾病        白内障        视网膜疾病        视网膜疾病

3 ?2 _0 W+ z5 N2 N9 x" E9 V* M4.3问题三模型的建立与求解
, M0 W, R8 a1 x. b+ O6 C# v4.3.1问题的分析
   首先利用Excel统计出8月1号至8月30号每种类型病人的等待时间 ，该类病人门诊时已住院人数  及待住院人数 ，下面以视网膜类病人为例，给出数据列表如下：

2 c& e9 Y- L9 R7 z
" M  M" d6 r1 g
* N1 A& k/ u0 B! S21        24        24        24        24        26        30        28        34        32        31        31        31        33        37
' o# K  P& O. ?+ E( v. w( ?( ]
  W  O7 P+ S1 m# x33        35        20        38        39        43        43        44        43        36        38        36        41        46        41

! y# s- Q8 X: v" ?/ j9 J% O# h12        12        13        12        11        13        13        14        14        13        12        12        14        13        14

+ m2 i9 ]1 `5 C4 P
% y! h8 j2 Y( I4 u" T5 `41        39        38        37        37        37        38        42        40        38        36        36        36        34        34

6 Q; w8 u! E& D5 A: e* h+ R/ ~39        35        36        38        34        41        37        31        31        31        34        34        35        32        35

: z% G% v) _" K& \4 L/ k1 U13        14        14        16        16        16        16        16        16        16        16        16        16        16        16
(表1)
利用matlab画出 与 和 的散点图及三次拟合曲线如下：
' Q0 s# g1 K9 G: u% H  Z# q2 G6 b
, o! z! x1 G5 o      
, P+ x! K4 m8 [' O; b4.3.2多元多项式回归模型
从图中可以看出该类病人的等待住院时间 随 及 的变化曲线为二次函数模型，并考虑到 与 之间的相互关系，我们建立了多元完全二次多项式回归模型：

利用matlab统计工具箱中的rstool命令求解，并且以交互式画面给出 的估计值 和预测区间，模型的回归曲线及回归系数的估计值为
+ ~0 N7 }0 X, v
9 Y! G+ I9 V" T/ p/ H
" w- X3 ~. F$ R' j- }: J
; x4 B. D# W4 a" y" E0 H

6 z; L# Y# g; I' `1 |0 a* x

( B9 ?5 F( @; s) T- r



' G6 x! V3 Z: V- z( ^4 N-0.5029        1.3368        0.0263        -0.0106        -0.0120        -0.0027
/ l4 f' ^; J% h6 Z$ `# B% d(表2)
0 n& P7 `* m& K& j6 H4 @故得出回归多项式为：
( x5 f& i) c& P$ x" U# y& O8 M, f
, I9 J6 z' [$ }- [1 d; x以此多项式我们验证9月1号数据。
# \" F( }' d. J6 q当 =32， =35时， 的估计区间为[12.8956，16.0396],即9月1号门诊的视网膜病人等待住院的时间区间为[12.8956，16.0396]
类似地，我们得出外伤病人的回归曲线方程及交互式图形如下：

0 b9 L% _+ @: n, w/ w  d
青光眼病人的回归曲线方程：

* O3 ^0 b" Z3 l; ?' i白内障病人的回归曲线方程及交互式曲线：
* Z  q4 h7 Y5 N- b: K1 {
4 f% J5 m! g. H3 ?) Z1 U具体回归曲线方程及交互式图形见附件1
根据不同类型病人门诊时的已住院 及待住院人数 代入各自的回归模型，即可得到该人等待住院的时间区间，以9月1号为例，各类病人等待时间区间为：
8 R" z% w3 p4 {6 w& g- m: d( y" @) c
视网膜        外伤        青光眼        白内障
2 F( ^4 X) G' W. e7 h[12.8956，16.0396]        [0.59786，1.3602]        [11.2415，18.5579]        [9.7619，16.1473]

(表3)
0 S: s- m% u; o6 Q$ ]  u$ ]$ Q4.4问题四模型的建立与求解
% _  S0 u% B* X0 F8 f
4.4.1问题的分析与模型的建立
    为考虑手术时间对病床合理度的影响，我们以一周为一单位时间将问题二的模型修改如下：
! h% P- ^: v, y4 s5 S线性规划模型：
目标函数：
' Q! ]) ?3 w% H' S        % 使病床安排合理度最大
其中 分别表示住院及手术满意度矩阵， 数据见附件2
约束条件：
（1）       %每天安排住院人数为该天出院人数
（2）      %每人在该周内至多既住院又接受手术
（3）当 时，若 ,则 ，   
当 时，若 ,则 ，  
: p* A  ^) p5 B) n0 W当 时，若 ,则 ，     %对类型排序， “外伤人数”， “白内障人数”， “其他疾病”
* C3 i8 k7 I1 Q$ e, x+ r（4） ,     %周一、周三其他眼科疾病不手术
2 h$ X9 Q4 w1 ^' q利用Lingo编程（详见附件2）运行结果以9月1号至9月7号为例。


4.5问题五模型的建立与求解
4.5.1问题五的分析     
0 q4 ]$ `3 O, D- V( c    分析病人在眼科医院的以下过程（如图）
; W' ]! N$ ]; ]5 e& b* S0 k3 @
/ }5 \$ U3 G) b! E可以将这个以住院为输入流以出院为输出流的过程作为一个新系统，由于每天病人住院人数及出院人数均是随机的服从泊松分布，其随机过程具有马尔科夫性质，可以采用经典的马尔科夫排队模型来描述。因此，每种类型的病人从住院到出院的过程可以近似地采用若干个相对独立的马尔科夫排队模型来描述.
4.5.2 M|M|C模型
! O7 M. H) Z, @; B# y利用Excel统计各类病人的平均住院率 ，系统的日平均服务率为 ，根据M|M|C(C≥2)模型，将病床视为服务台，病床之间是并联的，且独立工作，其每类病床的服务率相同，则系统的第 类病床的平均服务率为 ，在统计平衡状态下，服务强度为
7 f' }. [7 h5 M: p% l) B2 Q
" G) Y9 \% k: s4 k" L1 g: L此时，系统的稳态概率为
/ _  u8 g! Z( Z$ ~
5 A; J8 i) D2 l; V! c/ HM|M|C模型的主要指标如下：
（1）        平均队列长
5 i% Y9 s2 m+ Z! G( Q# \/ C% ?
2 O" o  b! W& a: Y& f! \/ B（2）平均队长

/ M9 \- P2 h* n: C+ Y（3）患者在系统中平均逗留时间
: v" _& f' S9 R. _. t
) X$ t; X0 j* S  A1 }( }又根据 <1及总床位数为79的条件给出分给每类病人的床位数C的有效范围为
+ g# c/ P. P1 D$ S* d4 t病人类型        病床数C的有效取值范围
: `# h# V% B8 v/ G2 |; Z外伤        [7,31]
8 \0 b) K( b0 [7 t( l视网膜        [22,46]
白内障单眼        [7,31]
白内障双眼        [13,37]
青光眼        [8,32]
. z' X* m' @. y: U(表4)
5 c  w% O1 ?% r7 Q1 F0 _7 M* [, b用matlab编程，求解 （各类病人在系统中的平均逗留时间），下表所示为运算结果。(程序见附件3)
表
7 U4 z$ D- E3 T8 \9 U  V: p8 g9 ]+ d& Q外伤        ws        视网膜        ws        青光眼        ws        白内障双        ws        白内障单        ws
7        14.9995        22        26.6608        8        23.7473        13        15.517        7        10.5578
% ~0 H6 C" j/ T) J  Y& f8        9.2543        23        17.3014        9        14.1607        14        11.0998        8        6.7745
9        7.9056        24        14.8688        10        12.0058        15        9.7775        9        5.8366
10        7.3933        25        13.8186        11        11.18        16        9.2031        10        5.4762
- N: [) t5 S" p  Z* H2 k11        7.1719        26        13.273        12        10.8125        17        8.9174        11        5.3202
12        7.0735        27        12.9635        13        10.6417        18        8.7674        12        5.2511
13        7.0303        28        12.78        14        10.5625        19        8.6873        13        5.221
14        7.012        29        12.6689        15        10.5266        20        8.6447        14        5.2082
0 b" B$ N* _. F' @# q15        7.0045        30        12.6011        16        10.5109        21        8.6223        15        5.2031
16        7.0016        31        12.5598        17        10.5043        22        8.6109        16        5.2011
17        7.0006        32        12.5349        18        10.5016        23        8.6051        17        5.2004
18        7.0002        33        12.52        19        10.5006        24        8.6024        18        5.2001
19        7.0001        34        12.5113        20        10.5002        25        8.601        19        5.2
# Q% E  y; A- K7 S20        7        35        12.5062        21        10.5001        26        8.6005        20        5.2
21        7        36        12.5034        22        10.5        27        8.6002        21        5.2
) |) R. a/ w( O% |22        7        37        12.5018        23        10.5        28        8.6001        22        5.2
23        7        38        12.5009        24        10.5        29        8.6        23        5.2
24        7        39        12.5005        25        10.5        30        8.6        24        5.2
25        7        40        12.5002        26        10.5        31        8.6        25        5.2
26        7        41        12.5001        27        10.5        32        8.6        26        5.2
) H# Q; v; g9 Q& |2 c% }) Z27        7        42        12.5001        28        10.5        33        8.6        27        5.2
1 a2 P: J1 T  l: e( P3 P9 \28        7        43        12.5        29        10.5        34        8.6        28        5.2
29        7        44        12.5        30        10.5        35        8.6        29        5.2
30        7        45        12.5        31        10.5        36        8.6        30        5.2
31        7        46        12.5        32        10.5        37        8.6        31        5.2

  f5 X4 q1 a3 q; L- v+ K: G0 d1 h1 l! b (表5)
4.4.3非线性规划模型
# x2 _; z! K; b在系统构造时，除原有新系统外还需考虑每类病人入院前的平均等待入院时间，利用excel统计数据，然后利用非线性规划建立模型，构造在各类病人在占用病床比例大致固定的条件下，所有病人在系统中的平均逗留时间最短的病床比例分配模型。
设 ,其中 指各类病床分配比例
$ a& D4 K: X, a1 V* {( i) J( X& d目标函数：
min            指第 类病人的平均等待时间
. E4 D3 J6 Q1 M4 K; P7 }约束条件：

4 b' ?% D* H+ l. k4 t但此时当 时， 计算起来比较麻烦，为方便数据处理，我们也可将模型转化为0-1线性规划模型。
4.4.4 0-1线性规划模型
根据表5数据及题目要求构造0-1线性规划模型。
目标函数：
/ z, k4 [( n8 E4 ?* Cmin +       其中 为上表中的数据
约束条件：
  E( s5 }, d. z+ V  L4 us.t.   .
利用Lingo编程并运行得到结果如下，当病床按13：27：12：17：10比例分配时所有病人在系统中平均逗留时间最短。程序及运行结果详见附件4
, {. t& L4 A: z0 e
五模型的评价与分析
+ V( m/ A8 i( W1 H（1）层次分析模型：
优点：此模型能定性定量的反映各因素对合理度影响的大小。
8 R6 \, P, s' h, q( o( E* g缺点：由于时间所限，考虑因素不全面；权值存在主观性。
（2）有权综合合理度评价模型：
本模型用EXCEL处理了大量数据，利用各因素得分的综合值比较准确的、直观的反映出每个因素对某天病床安排的满意度。
/ B# i  d7 }' l6 Q2 O& |/ ^（3）线性规划模型：
优点：1）本模型适用性强。
2）利用计算机处理大量数据快捷、方便。
缺点：考虑的影响因素不全面，可进行改进。
（4）排队论：
: Z1 r# s, o3 @+ j创新点：排队论系统的选取。
" T2 ?5 ^. f1 Q! ?+ J/ S' ]2 N- V

5 X5 h" i  u/ n* a% _0 I参考文献：

[1] 蒋启源，谢金星，叶俊，数学模型（第三版），北京：高等教育出版社，2005
% h( N' L' B& b" K/ k[2] 苏金明，王永利，MATLAB7.0实用指南，北京：电子工业出版社，2004
9 m- T" t6 v0 ?! W7 l3 z5 j9 p[3] 袁新生，邵大宏，郁时炼，Lingo和Excel在数学建模中的应用，北京：科学出版社，2007
[4]岳德权，吕胜利，李静铂，一个修理工的M/M/N可修排队[J].燕山大学学报，2003，27（3）：197-203
  b# X: g/ L( W+ t1 C. K7 u

/ {1 H6 x9 F$ _) Y  t
2 E, ~% g- p- n/ h$ g% C5 N* |- ^4 _
, D4 D" S* u' E! _3 y2 [

" z0 k' \8 W, d* ^
5 A) a; w$ Y; b5 l
3 P( Y6 T; G, k2 i9 p* G



8 r  i5 [$ y5 n0 \% Q! L' }附件1：
) _& T' U0 {- \
3 \, |% J# W& k. ^x1q=[8 9 9 9 8 8 7 9 15 18 18 18 19 21 21 17 15 13 9 7 7 6 6 5 5 5 6 8 8 9];
- U$ ]  s( V+ X0 [% nx2q=[20 19 37 20 20 19 19 19 16 11 9 8 9 11 9 9 8 9 9 9 9 10 9 8 10 10 9 8 10 10];
yq=[12 12 0 12 11 12 15 15 14 13 12 0 12 13 0 13 15 15 15 15 15 15 13 13 13 13 13 13 13 13 ] ;
xq=[x1q' x2q']
. |7 P7 j: O/ M5 y% ^* irstool(xq,yq,'quadratic')

2 y  M: q$ Q+ A. C
1 y8 w' X0 @* k! ^* N: @# j  Zbeta6,rmse6
  O5 w7 n8 J/ T# W/ f
! B5 ]4 X( j- l; ]! L2 z9 Sbeta6 =
9 X  u: b" ^2 h. H7 g' ~; H
$ I( F0 X) h- r9 q5 O* H* @   13.6974
/ X: P- i6 _4 D" s. B, m    0.2513
   -0.3200
    0.1068
- L: o; g5 z/ _& H3 W   -0.0648
   -0.0253
. _% [' L" x3 n. {) @7 H
" z  u5 u! {  A& I& r, d& B( I% Z
rmse6 =

    2.9225
# C8 E) A8 ~9 g/ D; r' ~7 G: F0 d

x1b=[24 25 29 33 36 37 35 34 23 23 23 24 24 21 17 18 21 22 25 28 26 25 22 26 29 31 30 27 29 26];
5 x5 {2 J0 _  T8 E2 \! px2b=[46 46 45 43 41 47 51 50 45 44 41 47 52 56 52 53 45 52 52 46 48 51 49 45 57 48 48 48 48 42];
yb=[12 12 12 12 11 13 14 14 14 13 12 12 13 14 13 13 14 14 16 16 16 16 16 16 16 16 16 16 16 16];
xb=[x1b' x2b']
rstool(xb,yb,'quadratic')
beta6,rmse6
, o* [4 k/ t$ f( Q
3 q( k  ?2 v& {' C+ G2 kbeta6 =

7 @+ ]) q9 M9 w* P1 H4 y  -38.1820
    1.1368
# _& E# u6 T# T6 J' z! }; |% w# a; s    1.3568
    0.0063
   -0.0257
- u- r7 M- V" j( e% C& h, g   -0.0139


: M; `' V1 E: Armse6 =
7 e3 |) r  H9 a8 u& W
, Z4 n: v- D! y; E; ?  Q0 D    1.4672
6 W3 O- S/ l% B! j




' X* X2 J- f9 n: {' B1 v
; ]) m, ~2 B9 O9 v/ u

% u$ L. W! _) X% A! t8 F) W8 P附件2：
model:
9 I3 R& L+ ?  y7 y/ B$ dsets:
ren/1..66/:n;
/ u& ~) Q7 W( h( I9 ttian/1..7/:m;
link(ren,tian):a,b,z,s;
7 m6 e( @, Q1 q  y) V" q0 f% A0 ]endsets
4 Y  z, _; Z0 c; w1 V0 u3 \: ^data:
* u! m/ P* ~- ~; c" r, T( h- H3 Wm=6 2 5 9 13 17 10;
n=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 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66;
$ {0 g, A+ l/ aa=6.97441 6.68291 6.39141 6.09991 0 0 0
- B- G# H7 J% }& n  z2 p  0 6.97441 6.68291 6.39141 6.09991 0 0
8 ]- X, B0 h1 U! k0 w) u5 O3 j3 m  0 0 6.97441 6.68291 6.39141 6.09991 0
  0 0 0 6.97441 6.68291 6.39141 6.09991
) U3 F" w  g1 y# t' a  0 0 0 0 6.97441 6.68291 6.39141
/ k; z# x3 \2 G$ E4 D  0 0 0 0 6.97441 6.68291 6.391410
' |1 f; @4 b/ w2 I  0 0 0 6.97441 6.68291 6.39141
  0 0 0 0 0 6.97441 6.68291
! @# J/ J# ]- e' d  ……
  5.51261 5.22111 4.92961 4.63811 4.34661 4.05511 3.76361
  5.22111 4.92961 4.63811 4.34661 4.05511 3.76361 3.47211;
* Z) E  d  w' G4 T) N
b=6.97441 6.68291 6.39141 6.09991 0 0 0
. P0 {- a, ?5 I" N8 X, m  0 6.97441 6.68291 6.39141 6.09991 0 0
  0 0 6.97441 6.68291 6.39141 6.09991 0
  0 0 0 6.97441 6.68291 6.39141 6.09991
2 h* W+ a3 v0 Z4 R( u4 @1 P  0 0 0 0 6.97441 6.68291 6.39141
  0 0 0 0 6.97441 6.68291 6.391410
……
+ l* n* J# N# z# M  
8 x) U& R- x  w# w  5.22111 4.92961 4.63811 4.34661 4.05511 3.76361 3.47211;
7 V# j; F9 v: K4 Senddata

max=@sum(link:a*z+b*s);
: a2 @8 `& o/ l& |
@for(ren(i)sum(tian(j):z(i,j))<=m(j));
4 {* A. k/ Z- G8 X) J@for(tian(j)sum(ren(i):z(i,j)+s(i,j))>=0);
6 f7 @. _4 K+ H' ]' m% `# s@for(tian(j)sum(ren(i):z(i,j)+s(i,j))<=2);
@if(n(i)#le#8,@if(link(i,j):z<=s,s(i,j)=0,n(i)#le#21,@if(link(i,j):z<=s+1,s(i,j)=0,n(i)#le#37,@if(link(i,j):z<=s+1,s(i,j)=0,n(i)#le#60,@if(link(i,j):z<=s+2,s(i,j)=0,n(i)#le#66,,@if(link(i,j):z<=s+2,s(i,j)=0))))));
……
End
附件3：
% x0 ?- w3 s1 b  C: J: M- z( h以白内障单眼c=7程序数据为例
( [" {; @" q: ~; h3 _0 J5 b
程序：
: O0 w9 Z* e) p/ G; Z" M" ~4 b- na=1.2;
& g# G. B; J( }$ B3 p. Bb=1/5.2;
k=[0:c-1];
x=sum(1./gamma(k).*((a/b).^k));
4 V& s( r" c; sy=1./gamma(c)*(a/b).^c*b*c/((c*b)-a);
z=sum(x+y);
p0=inv(z)
  [( h' X6 ]( }( w9 J; V/ blq=(a/b).^c/(gamma(c)*(1-a/(c*b)).^2)*a/(c*b)*p0
& A. `9 W; l! m* ?ls=lq+a/b
ws=ls/a
* Y3 `1 ]1 m5 r# e+ s1 Cwq=lq/a
- n2 G2 `8 H& k6 F" _
运行结果：
4 L' C# V5 p& h' G# P
c =

$ G% o6 B2 P8 p+ k/ X# k0 o" C; W2 q     7
2 }% r; Q5 H1 ?
6 P+ t6 R5 J" Q$ w6 U$ u, t
2 ^) W9 k7 h3 P/ I) lp0 =
0 r2 ]" g, K% `
" k% R( h9 _% N9 _& h5 B$ u1 H  1.6617e-004
  y' J# I( k1 N  \, \/ f* |
  H5 ~5 y2 a: ]# J% A0 u
) I. L/ E9 Q4 w3 a8 p3 Y% x1 d  Ylq =

* E9 i; \2 x: f; I6.4294

" Q# S0 M; ^" y" L$ Zls =
0 ~0 {" x* I# j. z; ?
   12.6694
6 F# v* I1 m( d5 @1 j6 Y, q7 j
- Z  A% s6 G, u6 ^) f# L: z. Xws =

   10.5578
. U6 o8 w3 p8 _3 s2 E9 r. C
wq =

6 L5 b! J( k  M% R3 U    5.3578
* Z; M  L0 v- D, Z, }6 n+ b: T
# l6 a" E/ r  T* g, K附件4：
model:
2 k: B  W$ W9 r8 Y  cMin=14.9995*x11+9.2543*x21+7.9056*x31+7.3933*x41+7.1719*x51+7.0735*x61+7.0303*x71+7.012*x81+7.0045*x91+7.0016*x101+7.0006*x111+7.0002*x121+7.0001*x131+7*x141+7*x151+7*x161+7*x171+7*x181+7*x191+7*x201+7*x211+7*x221+7*x231+7*x241+7*x251+
26.6608*x12+17.3014*x22+14.8688*x32+13.8186*x42+13.273*x52+12.9635*x62+12.78*x72+12.6689*x82+12.6011*x92+12.5598*x102+12.5349*x112+12.52*x122+12.5113*x132+12.5062*x142+12.5034*x152+12.5018*x162+12.5009*x172+12.5005*x182+12.5002*x192+12.5001*x202+12.5001*x212+12.5*x222+12.5*x232+12.5*x242+12.5*x252+
/ T% q; X7 E) ^0 x23.7473*x13+14.1607*x23+12.0058*x33+11.18*x43+10.8125*x53+10.6417*x63+10.5625*x73+10.5266*x83+10.5109*x93+10.5043*x103+10.5016*x113+10.5006*x123+10.5002*x133+10.5001*x143+10.5*x153+10.5*x163+10.5*x173+10.5*x183+10.5*x193+10.5*x203+10.5*x213+10.5*x223+10.5*x233+10.5*x243+10.5*x253+
/ \' I9 z5 L0 `9 g9 U$ Q15.517*x14+11.0998*x24+9.7775*x34+9.2031*x44+8.9174*x54+8.7674*x64+8.6873*x74+8.6447*x84+8.6223*x94+8.6109*x104+8.6051*x114+8.6024*x124+8.601*x134+8.6005*x144+8.6002*x154+8.6001*x164+8.6*x174+8.6*x184+8.6*x194+8.6*x204+8.6*x214+8.6*x224+8.6*x234+8.6*x244+8.6*x254+
10.5578*x15+6.7745*x25+5.8366*x35+5.4762*x45+5.3202*x55+5.2511*x65+5.221*x75+5.2082*x85+5.2031*x95+5.2011*x105+5.2004*x115+5.2001*x125+5.2*x135+5.2*x145+5.2*x155+5.2*x165+5.2*x175+5.2*x185+5.2*x195+5.2*x205+5.2*x215+5.2*x225+5.2*x235+5.2*x245+5.2*x255+
$ [3 D, L5 C9 D1 ]
1/79*(7*x11+8*x21+9*x31+10*x41+11*x51+12*x61+13*x71+14*x81+15*x91+16*x101+17*x111+18*x121+19*x131+20*x141+21*x151+22*x161+23*x171+24*x181+25*x191+26*x201+27*x211+28*x221+29*x231+30*x241+31*x251)+
7 N* f8 c7 D; D! p8 U# F12.7/79*(22*x12+23*x22+24*x32+25*x42+26*x52+27*x62+28*x72+29*x82+30*x92+31*x102+32*x112+33*x122+34*x132+35*x142+36*x152+37*x162+38*x172+39*x182+40*x192+41*x202+42*x212+43*x222+44*x232+45*x242+46*x252)+
! l( d4 U6 t4 Z12.3/79*(8*x13+9*x23+10*x33+11*x43+12*x53+13*x63+14*x73+15*x83+16*x93+17*x103+18*x113+19*x123+20*x133+21*x143+22*x153+23*x163+24*x173+25*x183+26*x193+27*x203+28*x213+29*x223+30*x233+31*x243+32*x253)+
12.7/79*(13*x14+14*x24+15*x34+16*x44+17*x54+18*x64+19*x74+20*x84+21*x94+22*x104+23*x114+24*x124+25*x134+26*x144+27*x154+28*x164+29*x174+30*x184+31*x194+32*x204+33*x214+34*x224+35*x234+36*x244+37*x254)+
) _0 d1 X/ e1 M4 q* O. `+ O+ H12.7/79*(7*x15+8*x25+9*x35+10*x45+11*x55+12*x65+13*x75+14*x85+15*x95+16*x105+17*x115+18*x125+19*x135+20*x145+21*x155+22*x165+23*x175+24*x185+25*x195+26*x205+27*x215+28*x225+29*x235+30*x245+31*x255);

x11+x21+x31+x41+x51+x61+x71+x81+x91+x101+x111+x121+x131+x141+x151+x161+x171+x181+x191+x201+x211+x221+x231+x241+x251=1;

  S. \- b# B- ]1 n- @X12+x22+x32+x42+x52+x62+x72+x82+x92+x102+x112+x122+x132+x142+x152+x162+x172+x182+x192+x202+x212+x222+x232+x242+x252=1;
/ t* z0 C" Q3 Q- r) Y0 t+ i) v
6 l" E6 o5 K: k  p/ P! e% nX13+x23+x33+x43+x53+x63+x73+x83+x93+x103+x113+x123+x133+x143+x153+x163+x173+x183+x193+x203+x213+x223+x233+x243+x253=1;

X14+x24+x34+x44+x54+x64+x74+x84+x94+x104+x114+x124+x134+x144+x154+x164+x174+x184+x194+x204+x214+x224+x234+x244+x254=1;

X15+x25+x35+x45+x55+x65+x75+x85+x95+x105+x115+x125+x135+x145+x155+x165+x175+x185+x195+x205+x215+x225+x235+x245+x255=1;
! B/ V# g( z8 u# W; j) G
: {3 s3 n0 n# x( |" {7*x11+8*x21+9*x31+10*x41+11*x51+12*x61+13*x71+14*x81+15*x91+16*x101+17*x111+18*x121+19*x131+20*x141+21*x151+22*x161+23*x171+24*x181+25*x191+26*x201+27*x211+28*x221+29*x231+30*x241+31*x251+
22*x12+23*x22+24*x32+25*x42+26*x52+27*x62+28*x72+29*x82+30*x92+31*x102+32*x112+33*x122+34*x132+35*x142+36*x152+37*x162+38*x172+39*x182+40*x192+41*x202+42*x212+43*x222+44*x232+45*x242+46*x252+
" U5 {2 v$ \: g5 `$ y  n8*x13+9*x23+10*x33+11*x43+12*x53+13*x63+14*x73+15*x83+16*x93+17*x103+18*x113+19*x123+20*x133+21*x143+22*x153+23*x163+24*x173+25*x183+26*x193+27*x203+28*x213+29*x223+30*x233+31*x243+32*x253+
13*x14+14*x24+15*x34+16*x44+17*x54+18*x64+19*x74+20*x84+21*x94+22*x104+23*x114+24*x124+25*x134+26*x144+27*x154+28*x164+29*x174+30*x184+31*x194+32*x204+33*x214+34*x224+35*x234+36*x244+37*x254+
7*x15+8*x25+9*x35+10*x45+11*x55+12*x65+13*x75+14*x85+15*x95+16*x105+17*x115+18*x125+19*x135+20*x145+21*x155+22*x165+23*x175+24*x185+25*x195+26*x205+27*x215+28*x225+29*x235+30*x245+31*x255=79;

0 E; J8 q9 w/ z! f# E3 W% XX11>=0;x21>=0;x31>=0;x41>=0;x51>=0;x61>=0;x71>=0;x81>=0;x91>=0;x101>=0;x111>=0;x121>=0;x131>=0;x141>=0;x151>=0;x161>=0;x171>=0;x181>=0;x191>=0;x201>=0;x211>=0;x221>=0;x231>=0;x241>=0;x251>=0;
x12>=0;x22>=0;x32>=0;x42>=0;x52>=0;x62>=0;x72>=0;x82>=0;x92>=0;x102>=0;x112>=0;x122>=0;x132>=0;x142>=0;x152>=0;x162>=0;x172>=0;x182>=0;x192>=0;x202>=0;x212>=0;x222>=0;x232>=0;x242>=0;x252>=0;
/ ]9 h$ d0 E, y& K% I% @3 ex13>=0;x23>=0;x33>=0;x43>=0;x53>=0;x63>=0;x73>=0;x83>=0;x93>=0;x103>=0;x113>=0;x123>=0;x133>=0;x143>=0;x153>=0;x163>=0;x173>=0;x183>=0;x193>=0;x203>=0;x213>=0;x223>=0;x233>=0;x243>=0;x253>=0;
x14>=0;x24>=0;x34>=0;x44>=0;x54>=0;x64>=0;x74>=0;x84>=0;x94>=0;x104>=0;x114>=0;x124>=0;x134>=0;x144>=0;x154>=0;x164>=0;x174>=0;x184>=0;x194>=0;x204>=0;x214>=0;x224>=0;x234>=0;x244>=0;x254>=0;
x15>=0;x25>=0;x35>=0;x45>=0;x55>=0;x65>=0;x75>=0;x85>=0;x95>=0;x105>=0;x115>=0;x125>=0;x135>=0;x145>=0;x155>=0;x165>=0;x175>=0;x185>=0;x195>=0;x205>=0;x215>=0;x225>=0;x235>=0;x245>=0;x255>=0;

end
) `8 c8 t2 i* U% v7 a   

1 N6 I  p) q) ^5 c5 Z3 j4 Q* Y结果：
) ~: g, g: `$ W1 L  Global optimal solution found.
  Objective value:                              55.91382
( _5 u6 n- R. G  Infeasibilities:                              0.000000
- m( A( s; n" T9 o- V  Total solver iterations:                            65
) b% D/ z+ r+ n" ?! U
: V2 b$ Q) X+ Y0 k9 j
                       Variable           Value        Reduced Cost
4 y4 v) D, I* P8 _                            X11        0.000000            7.756808
                            X21        0.000000            2.047006
( N' [: R- y$ P7 F2 U1 c- k                            X31        0.000000           0.7337051
                            X41        0.000000           0.2568038
( O% b# d: z" n7 l2 X                            X51        0.000000           0.7080253E-01
                            X61        0.000000           0.7801266E-02
                            X71        1.000000            0.000000
; j; _6 h" f3 g( f                            X81        0.000000           0.1709873E-01
                            X91        0.000000           0.4499747E-01
0 p# R& L3 ^; B% o7 Y# j                           X101        0.000000           0.7749620E-01
                           X111        0.000000           0.1118949
% i" O% n% U# e: t- A2 ^                           X121        0.000000           0.1468937
                           X131        0.000000           0.1821924
                           X141        0.000000           0.2174911
                           X151        0.000000           0.2528899
                           X161        0.000000           0.2882886
' D4 K/ M* |0 l7 i& b                           X171        0.000000           0.3236873
                           X181        0.000000           0.3590861
                           X191        0.000000           0.3944848
                           X201        0.000000           0.4298835
5 u+ q8 j4 O. Y: s                           X211        0.000000           0.4652823
                           X221        0.000000           0.5006810
, H8 h  f4 m/ y* i- t) a" Y                           X231        0.000000           0.5360797
                           X241        0.000000           0.5714785
2 d1 ]( T( J) L3 G, b! C                           X251        0.000000           0.6068772
$ r% }, w& q+ F" a                            X12        0.000000            12.77980
                            X22        0.000000            3.603900
                            X32        0.000000            1.354800
+ O0 G7 ~' ?' a( J& L                            X42        0.000000           0.4881000
: i- D2 q* V- K5 T                            X52        0.000000           0.1260000
                            X62        1.000000            0.000000
                            X72        0.000000            0.000000
                            X82        0.000000           0.7240000E-01
- |6 @# x/ T& w' D, A! @! O9 h+ o                            X92        0.000000           0.1881000
$ c0 C# Q' Y5 }5 G; P                           X102        0.000000           0.3303000
: R0 V: ]% k, R: y) j                           X112        0.000000           0.4889000
                           X122        0.000000           0.6575000
                           X132        0.000000           0.8323000
                           X142        0.000000            1.010700
5 Q$ w5 Z% m) v2 g# I                           X152        0.000000            1.191400
) v" W9 L6 }% I) c; X                           X162        0.000000            1.373300
                           X172        0.000000            1.555900
  w2 A+ a9 g3 h5 Z9 |: ~; l                           X182        0.000000            1.739000
$ e( q9 }# |0 d                           X192        0.000000            1.922200
                           X202        0.000000            2.105600
9 `2 X1 [4 W# G5 z( D$ m4 N; Z  l                           X212        0.000000            2.289100
7 W( l8 E. W( X/ e$ f0 a! W                           X222        0.000000            2.472500
6 b$ _% M0 M! p: U  Y8 U% D# i$ f                           X232        0.000000            2.656000
                           X242        0.000000            2.839500
0 X3 y, t# F# n5 z$ v                           X252        0.000000            3.023000
3 N3 K2 b+ ~0 m) P" S3 |& |                            X13        0.000000            12.22105
6 E6 U) M: I. Q" J& `3 ?7 I  Y                            X23        0.000000            2.812890
                            X33        0.000000           0.8364266
% Z  ]" ?5 v7 D; n1 t                            X43        0.000000           0.1890633
0 A3 z& V# ^& m. c0 q3 [                            X53        1.000000            0.000000
2 X, \- P+ v3 q" S, M5 l0 G! W                            X63        0.000000           0.7636709E-02
                            X73        0.000000           0.1068734
$ n' y& k' a! l) M+ ^  `/ X                            X83        0.000000           0.2494101
* t: ?6 G$ k) ~5 i                            X93        0.000000           0.4121468
                           X103        0.000000           0.5839835
3 F- K9 R! X8 L                           X113        0.000000           0.7597203
                           X123        0.000000           0.9371570
                           X133        0.000000            1.115194
                           X143        0.000000            1.293530
3 f2 Y9 q/ n: x3 a                           X153        0.000000            1.471867
' s! d/ \" \+ A; l                           X163        0.000000            1.650304
/ f1 w! O, W3 J3 k                           X173        0.000000            1.828741
% Y1 b' o# K; |( h                           X183        0.000000            2.007177
! Q. W' @- J& p# b5 A3 w8 c                           X193        0.000000            2.185614
% I7 ?) H( {& l- S, ]( ?5 I                           X203        0.000000            2.364051
4 f& n! |; w, Q2 V                           X213        0.000000            2.542487
                           X223        0.000000            2.720924
                           X233        0.000000            2.899361
                           X243        0.000000            3.077797
                           X253        0.000000            3.256234
                            X14        0.000000            5.865600
' S% O' Z+ T& A                            X24        0.000000            1.631900
                            X34        0.000000           0.4931000
* |0 g$ L1 L- |. w% Y8 d& s                            X44        0.000000           0.1022000
$ }% `; C" j; k# w* g+ E                            X54        1.000000            0.000000
+ _6 [/ L0 c7 E) T- N                            X64        0.000000           0.3350000E-01
                            X74        0.000000           0.1369000
; h! J1 P& I+ M* z* b* b                            X84        0.000000           0.2778000
                            X94        0.000000           0.4389000
  |8 T' B' U; U# {" E; E                           X104        0.000000           0.6110000
                           X114        0.000000           0.7887000
                           X124        0.000000           0.9695000
                           X134        0.000000            1.151600
# Q9 a+ r/ Q. K& r                           X144        0.000000            1.334600
                           X154        0.000000            1.517800
' g8 `" c' F" z9 J9 {; e8 B                           X164        0.000000            1.701200
                           X174        0.000000            1.884600
$ A9 k7 S( c, d3 D, ]9 x                           X184        0.000000            2.068100
                           X194        0.000000            2.251600
                           X204        0.000000            2.435100
+ N5 Q/ r! N# J0 M& Y% v                           X214        0.000000            2.618600
                           X224        0.000000            2.802100
, c7 p8 `) h5 i# k: W% F                           X234        0.000000            2.985600
                           X244        0.000000            3.169100
$ X( P! J  i8 `2 p/ B( r6 \                           X254        0.000000            3.352600
                            X15        0.000000            4.531100
# \! f, x( t  `8 \! c5 }                            X25        0.000000           0.9313000
) h; Y( z6 T8 C                            X35        0.000000           0.1769000
                            X45        1.000000            0.000000
& l( W; F: I# N8 r                            X55        0.000000           0.2750000E-01
                            X65        0.000000           0.1419000
2 v; w5 {  O( Y: P3 l, _                            X75        0.000000           0.2953000
; m  C- O5 q4 {4 e& R                            X85        0.000000           0.4660000
                            X95        0.000000           0.6444000
3 R8 N' i, y% p                           X105        0.000000           0.8259000
                           X115        0.000000            1.008700
                           X125        0.000000            1.191900
- ~) e. |: I" M0 P, l                           X135        0.000000            1.375300
  @/ X- D4 G' C' U                           X145        0.000000            1.558800
                           X155        0.000000            1.742300
                           X165        0.000000            1.925800
                           X175        0.000000            2.109300
                           X185        0.000000            2.292800
                           X195        0.000000            2.476300
                           X205        0.000000            2.659800
                           X215        0.000000            2.843300
                           X225        0.000000            3.026800
                           X235        0.000000            3.210300
                           X245        0.000000            3.393800
8 ]; t8 d; @2 \' U/ l2 [. J                           X255        0.000000            3.577300

                            Row    Slack or Surplus      Dual Price
                              1        55.91382           -1.000000
                              2        0.000000           -7.490484
: }$ d; @: U  R                              3        0.000000           -17.91800
                              4        0.000000           -12.95374
                              5        0.000000           -12.03690
9 `: q3 D- Y& B0 G7 |8 M                              6        0.000000           -7.311200
  u, M0 e9 e8 X0 o! |: |9 }# c                              7        0.000000           0.2274051E-01
% x9 A* q, R) W                              8        0.000000            0.000000
                              9        0.000000            0.000000
6 G: f6 Z8 D8 b$ W" X- B( N                             10        0.000000            0.000000
                             11        0.000000            0.000000
                             12        0.000000            0.000000
                             13        0.000000            0.000000
                             14        1.000000            0.000000
                             15        0.000000            0.000000
                             16        0.000000            0.000000
1 b+ b" c! s+ K! e# k                             17        0.000000            0.000000
                             18        0.000000            0.000000
' E# C9 f, H% n6 H3 X; ]7 }7 C9 {: X' e                             19        0.000000            0.000000
                             20        0.000000            0.000000
                             21        0.000000            0.000000
* I/ C  _1 o. V/ W* p: I% C7 S$ n                             22        0.000000            0.000000
  s: R) }( G2 E                             23        0.000000            0.000000
, r% K0 t# u" ^* T                             24        0.000000            0.000000
* z' C) j% {# n, K( W& g                             25        0.000000            0.000000
- p* O, ?, j% x5 i' P! u                             26        0.000000            0.000000
                             27        0.000000            0.000000
                             28        0.000000            0.000000
9 t! X5 _9 x  w                             29        0.000000            0.000000
                             30        0.000000            0.000000
7 ~0 \, M% T/ Q+ X+ p. x1 k                             31        0.000000            0.000000
6 l. N  V, \6 b* s) X2 a! A% v                             32        0.000000            0.000000
                             33        0.000000            0.000000
                             34        0.000000            0.000000
: o) T& J: Y* v( |5 x- E; ?8 f                             35        0.000000            0.000000
: Q3 J3 ]: ^- N" p* C/ W- ^                             36        0.000000            0.000000
                             37        0.000000            0.000000
9 N- b2 D* |2 E3 G$ F1 @( o( r                             38        1.000000            0.000000
                             39        0.000000            0.000000
                             40        0.000000            0.000000
                             41        0.000000            0.000000
* l/ R2 R: o9 a3 k' ~                             42        0.000000            0.000000
                             43        0.000000            0.000000
, x7 @- T" @/ g; |( U                             44        0.000000            0.000000
                             45        0.000000            0.000000
2 `* `  ~5 u/ A" N' _" C                             46        0.000000            0.000000
                             47        0.000000            0.000000
' u, Z8 U1 n. X/ W4 W7 y7 E- e                             48        0.000000            0.000000
                             49        0.000000            0.000000
                             50        0.000000            0.000000
% H% F9 o' a  w5 e                             51        0.000000            0.000000
- h% P) f6 S( G) `0 `! J+ z  W                             52        0.000000            0.000000
                             53        0.000000            0.000000
& O) P8 r- R$ L+ P+ d' K                             54        0.000000            0.000000
# m; _* E: i9 ^* ^" D                             55        0.000000            0.000000
( `$ r- }/ P& t% I/ g                             56        0.000000            0.000000
' M( Q( G& L3 V' w! X                             57        0.000000            0.000000
                             58        0.000000            0.000000
& C' ^# S( K: R6 N( {                             59        0.000000            0.000000
                             60        0.000000            0.000000
                             61        0.000000            0.000000
                             62        1.000000            0.000000
6 e+ m2 V! b6 f" D3 ^2 r                             63        0.000000            0.000000
                             64        0.000000            0.000000
                             65        0.000000            0.000000
                             66        0.000000            0.000000
# B" }+ H) @8 F$ [9 W- b                             67        0.000000            0.000000
                             68        0.000000            0.000000
3 H+ y/ \: ^3 }1 X/ W0 g6 s6 d  K; |                             69        0.000000            0.000000
# W; d& S  ?4 F1 ~; p                             70        0.000000            0.000000
9 _5 k7 N$ O: ]                             71        0.000000            0.000000
3 B$ e9 W/ j5 N4 W( Z1 [! C. i0 B7 `                             72        0.000000            0.000000
                             73        0.000000            0.000000
7 c. E6 @* h9 T6 B3 w" R5 x                             74        0.000000            0.000000
                             75        0.000000            0.000000
0 Y, O, l" U* P' _( N7 R( s0 T                             76        0.000000            0.000000
3 C" N/ J( h9 w5 B                             77        0.000000            0.000000
                             78        0.000000            0.000000
. c3 ]4 [/ U, n% g4 b8 I- x+ G1 F                             79        0.000000            0.000000
                             80        0.000000            0.000000
3 m5 E( t' ~( {5 g- ^  y7 Y6 u: q5 T                             81        0.000000            0.000000
: T: G- Z0 ?% \% q% Q. Z                             82        0.000000            0.000000
                             83        0.000000            0.000000
( _7 u: _0 Y& E# J                             84        0.000000            0.000000
                             85        0.000000            0.000000
                             86        0.000000            0.000000
* P+ G% d0 f) E! Q/ G                             87        1.000000            0.000000
                             88        0.000000            0.000000
                             89        0.000000            0.000000
                             90        0.000000            0.000000
                             91        0.000000            0.000000
                             92        0.000000            0.000000
9 Z. W# J$ J, r0 g0 H3 O/ b# @                             93        0.000000            0.000000
                             94        0.000000            0.000000
                             95        0.000000            0.000000
1 V& X4 V6 ]0 t( J                             96        0.000000            0.000000
                             97        0.000000            0.000000
! q9 w4 ]/ h" H& ^5 u) c+ U8 l                             98        0.000000            0.000000
2 x3 z  R$ i8 y, V, \                             99        0.000000            0.000000
4 o6 P7 Q4 ~; ~                            100        0.000000            0.000000
                            101        0.000000            0.000000
                            102        0.000000            0.000000
0 c' A9 h. d* {                            103        0.000000            0.000000
) o8 V' ?3 k% X                            104        0.000000            0.000000
                            105        0.000000            0.000000
0 G/ g$ F% l7 X% J% a9 ], q  I                            106        0.000000            0.000000
  |& @5 c3 D4 I+ f1 o" \                            107        0.000000            0.000000
  S4 |3 F; f% w' W                            108        0.000000            0.000000
                            109        0.000000            0.000000
% H3 l) L, q! L1 q2 o                            110        0.000000            0.000000
                            111        1.000000            0.000000
                            112        0.000000            0.000000
0 E! g/ A* |7 B0 ~% \                            113        0.000000            0.000000
                            114        0.000000            0.000000
7 X3 \  [! K8 }! h, Q* j                            115        0.000000            0.000000
' e% }$ C; t+ g: F; B" r, J/ \1 Y. E3 a                            116        0.000000            0.000000
4 t& o/ c7 W& @, t, B                            117        0.000000            0.000000
4 C' g' z) [1 h+ d3 o! H                            118        0.000000            0.000000
  O" d# R' i0 A                            119        0.000000            0.000000
                            120        0.000000            0.000000
                            121        0.000000            0.000000
                            122        0.000000            0.000000
6 r) W) R; K+ ]: o- \. s$ ?1 v                            123        0.000000            0.000000
                            124        0.000000            0.000000
; L: V% N; f* Y( d$ b+ t1 y+ d                            125        0.000000            0.000000
" i1 W$ i  C6 I" t                            126        0.000000            0.000000
                            127        0.000000            0.000000
                            128        0.000000            0.000000
                            129        0.000000            0.000000
                            130        0.000000            0.000000
% h$ a; {" k6 R& S                            131        0.000000            0.000000
$ L! d% V& F. o' V                            132        0.000000            0.000000