长江水质的评价与预测.doc

长江水质的评价和预测摘要在长江水质日趋下滑的情况下，经过对长江水质的评价和预测，对长江水质污染问题提出合理的建议和意见。针对不同的问题，建立不同的模型，提出不同的解决方案。对于问题1，为不同的水质设置不同的权重值，经过分析对水质做出定量的综合评价；对各地区水质的污染情况，在对每一个地区在近两年多的28个月中水质情况统计的基础上，计算出该地区污染的种类及各种类污染出现的频率，通过计算分析每个地区的污染情况。对于问题2，用Matlab程序解出时间矩阵T，用污染物的含量减去降解后的剩余量得到自排量。比较每一个观测点的自排量，可以得出长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区。对于问题3，运用灰色理论的GM（1，1）模型，解得关联度e =0.9993，灰色关联度精度等级为1级（好），最终得出如不治理，长江水质将会越来越差。对于问题，取19952004年文年水的IV和V类水质分别所占百分比和长江的年总流量作为预测原始值，预测出20052014年的IV和V类水质分别所占百分比和长江的年总流量。设长江污水年排放量与长江的年总流量的比作为因变量Y，利用spss 软件的分析功能，可以估算出未来10年的年处理污水量。问题，在问题1对长江近两年多的水质情况定量的综合评价，并分析各地区水质的污染状况；问题2研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区；问题3对长江未来水质污染的发展趋势做出预测分析的基础上，对解决长江水质污染问题提出一些切实可行的建议和意见。关键词：水质问题 权重 灰色模型预测法 变量剔除 线性回归一、问题重述长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。2004年10月，由全国政协与中国发展研究院联合组成“保护长江万里行”考察团，从长江上游宜宾到下游上海，对沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心。为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃”，并发出了“拿什么拯救癌变长江”的呼唤。通常认为一个观测站的水质污染主要来自于本地区的排污和上游的污水。一般说来，江河自身对污染物都有一定的自然净化能力，即污染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。反映江河自然净化能力的指标称为降解系数。事实上，长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.10.5之间，比如可以考虑取0.2(单位：1/天)。现需讨论下列问题：（1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?（3）假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。（4）根据预测分析，如果未来10年内每年都要求长江干流的类和类水的比例控制在20%以内，且没有劣类水,那么每年需要处理多少污水？（5）对解决长江水质污染问题有什么切实可行的建议和意见。二、问题分析 (一)问题一的分析问题一可以分解为两个方面进行讨论。对于近两年多水质情况的定量评价方面，引进评价指数和权。将三类水质在不同流域的评价指数设为ij，各类水质的百分比设为i；把I、II、III、IV、V和劣V水质的权分别记为6、5、4、3、2、1。计算出i和j，通过分析比较，定量做出水质的综合评价。对于各地区水质的污染状况，在上述对每一个地区在近两年多的28个月中水质情况统计的基础上，再找出该地区污染的种类及各种类污染出现的频率，建立表格，通过计算分析得出结论。(二)问题2的分析设7的观测点分别为， Si为i到i+1观测点的距离，模型里假设各段之间水流为匀速，以两观测点的平均速度作为Si的速度，则 ，。用Matlab程序解出时间矩阵T，高锰酸盐量（或氨氮含量）Q=c*q，设降解后的剩余量为，则观测点的自排量M=Q-P。比较每一个观测点的自排量，可以得出长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区。(三)问题3的分析运用灰色预测的原理，建立GM (1, 1) 模型，编写程序，得出未来年的水质情况。(四)问题4的分析取19952004年文年水的IV和V类水质分别所占百分比和长江的年总流量作为预测原始值，预测出20052014年的IV和V类水质分别所占百分比和长江的年总流量。设长江污水年排放量与长江的年总流量的比作为因变量Y，利用spss 软件的线性回归分析功能，对19952004 年的每年长江干流水文年各类水质所占河长的百分比做为历史数据进行回归分析，得到各个因子的常量系数，即可得到Y的表达式，由此得出，且未来的年处理污水量=期望长江年总流量，经过上述方法可以估算出未来10年的年处理污水量。(五)问题5的分析在问题1对长江近两年多的水质情况定量的综合评价，并分析各地区水质的污染状况；问题2研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区；问题3对长江未来水质污染的发展趋势做出预测分析的基础上，对解决长江水质污染问题提出一些切实可行的建议和意见。三、模型假设1. 每一个观测站的水质污染物主要来自本地区的排污河上游的污水，河水在通过观测点后立即混入下段区域的污水；2. 降解作用在整个流域中是平稳进行的，主要污染物高锰酸盐指数和氨氮的降解系数取0.2（单位：1/天）；3.不考虑洪水干旱等特殊气候对水质的影响；42005年2014年的评价河长与2004年的评价河长相等；5.任意两观测点之间的江水流速是所观测的两点流速的平均值。四、定义与符号说明：三类水质在不同流域的评价指数（）；：第i类水质评价指数的平均值；：第j类流域评价指数的平均值；：各类水质的百分比（）；：第i个观测点（）；：第i个观测点到i+1个观测点的距离（）；：第i个观测点的水流速度（）；：从i观测点到i+1观测点的时间（）；：第i段流域的长度（）；：第i个观测点自排量（）；：降解系数=0.2/天；：浓度；：水流量；：年废水排放总量。五、模型的建立与求解（一）问题一模型的建立和求解：对于问题1，根据附件4中03年、04年的报告表，运用权重的表达式： 计算出百分比如下表1： 年份时段全流域干流支流平均值2003年枯水期3.7223.7883.7653.758丰水期3.8354.0713.7873.898水文年4.0204.2223.9784.0752004年枯水期3.6733.5483.6993.640丰水期3.7033.6783.6833.688水文年3.6883.6973.6843.690表1 近两年来长江流域水质报告图表2 2003年长江流域水质折线图表3 2004年长江流域水质折线图对于问题2，根据附件3近两年的28个月中的水质情况进行统计，再找出该地区污染的种类及该种类污染出现的频率。统计结果如下表4：观测站pHDOCODMnNH3-N四川攀枝花龙洞001/281/28重庆朱沱0000湖北宜昌南津关0000湖南岳阳城陵矶0000江西九江河西水厂0000安徽安庆皖河口0000江苏南京林山0000四川乐山岷江大桥06/284/289/28四川宜宾凉姜沟001/283/28四川泸州沱江二桥05/281/285/28湖北丹江口胡家岭0000湖南长沙新港01/28010/28湖南岳阳岳阳楼003/280湖北武汉宗关01/2800江西南昌滁槎1/288/282/2821/28江西九江蛤蟆石004/280江苏扬州三江营01/2800表4 各地污染频率统计表从上表可以看出江西南昌滁槎出现了PH值超标情况；四川乐山岷江大桥，四川泸州沱江二桥，江苏扬州三江营，湖北武汉宗关，江西南昌滁槎，湖南长沙新港这6个地区出现了含氧量超标情况；四川攀枝花龙洞，四川乐山岷江大桥，四川宜宾凉姜沟，四川泸州沱江二桥，湖南岳阳岳阳楼，江西南昌滁槎，江西九江蛤蟆石这7个地区出现了高锰酸盐超标情况；四川攀枝花龙洞，四川乐山岷江大桥，江西南昌滁槎，四川泸州沱江二桥，湖南长沙新港，四川宜宾凉姜沟这6给地区出现了氨氮超标情况。经过上述分析统计，结合图表可以看出2003年到2004年的水质情况下降，而且在水文年的水质较好，枯水期的水质最差。四川攀枝花龙洞CODMn，NH3-NC未达标率分别为：1/28，1/28，即该地区在抽样的28个月内有1个月出现高锰酸盐超标情况，1个月出现氨氮超标情况。四川乐山岷江大桥DO，CODMn，NH3-N的未达标率分别为：6/28，4/28，9/28，即该地区在抽样的28个月内有6个月出现含氧量超标情况，4个月出现高锰酸盐超标情况，9个月出现氨氮超标情况。四川宜宾凉姜沟CODMn，NH3-NC未达标率分别为：1/28，3/28，即该地区即该地区在抽样的28个月内有1个月出现高锰酸盐超标情况，3个月出现氨氮超标情况。四川泸州沱江二桥DO，CODMn，NH3-N的未达标率分别为：5/28，1/28，5/28，即该地区在抽样的28个月内有5个月出现含氧量超标情况，1个月出现高锰酸盐超标情况，5个月出现氨氮超标情况。湖南长沙新港DO，NH3-N的未达标率分别为：1/28，10/28，即该地区在抽样的28个月内有1个月出现含氧量超标况，10个月出现氨氮超标情况。湖南岳阳岳阳楼CODMn，未达标率为：3/28，即该地区即该地区在抽样的28个月内有3个月出现高锰酸盐超标情况。湖北武汉宗关湖南岳阳岳阳楼DO未达标率为：1/28，即该地区即该地区在抽样的28个月内有1个月出现含氧量超标情况。江西南昌滁槎PH，DO，CODMn，NH3-NC未达标率分别为：1/28，8/28，2/28，21/28即该地区在抽样的28个月内有1个月出现PH超标情况，8个月出现含氧量超标情况，2个月出现高锰酸盐超标情况，21个月出现氨氮超标情况。江西九江蛤蟆石CODMn，未达标率为：4/28，即该地区即该地区在抽样的28个月内有4个月出现高锰酸盐超标情况。湖北武汉宗关湖南岳阳岳阳楼DO未达标率为：1/28，即该地区即该地区在抽样的28个月内有1个月出现含氧量超标情况。（二）问题二模型的建立和求解：模型里假设各段之间水流为匀速，以两观测点的平均速度作为Si的速度，则，记对应的流域段长度为，则。分别运用Matlab程序解出：；速度矩阵江水流过所需时间矩阵单位时间内高锰酸盐量（单位：g）：，运用Matlab程序解得：单位时间内氨氮量（单位：g）：，运用Matlab程序解得：上游高锰酸盐经过降解到下游时剩余高锰酸盐量，上游氨氮经过降解到下游相邻观测点时剩余氨氮量，运用Matlab程序解得：Q1 = 19.0, 3.09, 13.0, 4.68, .443e4, 45.0 28.6, .627, 10.3, 6.62, .571e4, 5.88 27.3, 6.50, 213., 103., .977e4, 196. 44.1, 21.3, 456., 203., .985e4, 290. 53.6, 110., 385., 168., .946e4, 392. .108e4, .243e4, .305e4, .200e4, .456e5, .934e4 3.73, 12.6, 633., 230., .740e4, 538. 1.92, .234, 1.10, .542, .108e4, 4.20 .694, .487e-1, .968, .331, 669., .917 .458e-1, .268e-2, .135e-1, .251e-1, 334., .316 .415e-2, .624e-5, .369e-2, .154e-1, 407., .366 .223e-2, .173e-5, .445e-2, .478e-1, 568., 2.49 .156e-1, .922e-4, .962e-3, .384e-2, 467., .948Q2 = 1.24, .185, 1.62, .369, 412., 3.53 .466, .770e-1, 1.03, .589, 405., 4.63 .436, .390, 19.1, 10.3, 573., 20.8 .734, 1.16, 22.8, 17.4, 606., 37.5 9.24, 9.34, 23.6, 13.9, .113e4, 33.0 15.9, 88.2, 260., 154., .599e4, 817. .373, 1.81, 39.3, 23.7, 740., 43.5 .410e-1, .316e-1, .139, .751e-1, 108., .229 .347e-1, .191e-1, .645e-1, .293e-1, 34.9, .755e-1 .267e-2, .111e-2, .156e-2, .207e-2, 29.2, .243e-1 .691e-3, .191e-5, .332e-3, .206e-2, 17.1, .244e-1 .528e-3, .493e-6, .402e-3, .490e-2, 34.7, .188 .142e-2, .291e-4, .521e-4, .556e-3, 32.7, .619e-1任意地区自身排放高锰酸盐和氮氨排放量分别为即为某地区高锰酸盐和氮氨总量减去上游相邻城市经过降解后剩下即可，分别如下：，利用matlab解得：如下所示：时间A1A2A3A4A5A6A720040484874828146197844678446778170446552004051599628791593997379092373509062094200406100254257377134788877529741590660042004071118454076726191.01E+056999744380575102004082169221146742901.03E+05640325400059278200409383082.08E+051.79E+052.07E+051.36E+05731201.77E+05200410260825916706577741746890577007242220041142001633820140311993211936320352762004121521.68514.31776037449303602330129799200501854.45226868341769261603164636000200502550.86485.4902025536319303658334730200503685.390061087828864257406450238698200504706.273001296020996302004599333149平均值8986.02337077.950708.2370086.4657375.7745568.8557428.85如下所示：时间A1A2A3A4A5A6A7200404553.52896.85759.86654.47586.66078592.47200405260.43536.55939.965596555.480953100.4200406160.42555.66901.67890.94414.748841965.2200407186.42951.33630.88653.24302.664121567.520040837401792.84550.78523.47690.142141722200409565.27600.115427158801804643985663200410260.83725.64391.27988.74688.345301492.5200411902205.925444319.93211.91932979.7720041276.083340.511843316.915841939.11489.920050149.842170.81005.43449.822892430.8316820050291.81981.6811.83431.413392448.94681200503161.982559.6984.22956.815734887.3429.8120050464.22299.57023040.821142997.34419.9平均值481.58463047.4314140.9546358.8625030.3544249.7232405.496通过比较平均值的大小，得出高锰酸盐超标的污染A4污染最严重，A7次之，A5为第三，A1最轻。氨氮超标的污染A4最严重，A5次之，A3为第三，A1最轻。（三）问题三模型的建立和求解：GM (1, 1) 模型：灰色系统理论的实质是将无规律的原始数据进行累加生成，得到规律性较强的生成数列后再重新建模。由生成模型得到的数据再通过累加生成的逆运算累减生成得到还原模型，再还原模型作为预测模型。灰色模型是预测工作的基础模型。对于采集的原始数据序列:记为原始序列，为由经过一次累加生成的序列，其中，表示的均值生成序列， 表5通过运用灰色原理，GM模型解得关联度e =0.9993，灰色关联度精度等级为1级（好），画出图像。最终得出未来10年的值分别为：4.3314 ，4.2580，4.1857，4.1147，4.0449，3.9763，3.9089，3.8425，3.7774，3.7133。从整体分析，如果不加治理，长江水质在未来10年会持续下降。（四）问题四模型的建立和求解：取19952004年文年水的IV和V类水质分别所占百分比和长江的年总流量作为预测原始值，预测出20052014年的IV和V类水质分别所占百分比和长江的年总流量。长江污水年排放量与长江的年总流量的比作为因变量Y，利用spss 软件的线性回归分析功能对19952004 年的每年长江干流水文年各类水质所占河长的百分比做为历史数据进行回归分析，得到各个因子的常量系数，分别为系统自动剔除第三变量，从而得到：其中由此得出，且未来的年处理污水量=期望长江年总流量年份期望长江总流量IV类河长所占百分比% ()V类河长所占百分比% ()需要处理的污水量（亿吨）2005102657.651.570.004344.13952006101908.351.950.005354.0072007101169.112.400.006565.7542008100429.932.970.008181.34022009996910.833.670.0199.692010989611.824.530.012118.7522011982412.895.600.015147.362012975314.066.920.019185.3072013968215.348.540.023222.6862014961116.7310.550.029278.719表6 每年需要处理的污水量（五）问题五对解决长江水质污染问题的一些切实可行的建议和意见：（1）对保护长江立法，让人们的意识受到法律的约束。一方面可以制定促进循环经济发展的政策和法律法规；另一方面，还可以加大对违法排污行为的处罚力度。（2）对于生活污水也需要进行一定的控制。唤醒民众、领导干部的环保意识。由于长江两岸农田过量喷洒的农药被雨水冲刷到河里，出现了氨氮含量过高、河水富营养化、蓝藻疯长等情况，严重破坏了水资源。迫切需要提高人民大众的环保意识，在保证收成的情况下，尽量少使用农药毋将农药水随便排放。同时需要政府拨款在长江沿岸多建水净化厂，将需要排放到长江里的水进行净化处理后再排放。（3）在长江上游地区植树造林，据统计长江上游地区森林覆盖率曾达到60%至80%，现如今却只有5%至7%，因此每年流入长江泥沙近24亿吨，长江源头出现沙漠化及断流，污染还殃及近海。2003年出现赤潮近40次，赤潮导致“氧亏”，进而影响渔业资源及品质，水生生态系统濒于崩溃。因此应该鼓励在长江上游地区植树造林。（4）对于长江沿线的钢铁厂、化工厂、造纸厂、造船厂、拆船厂这样的重污染企业禁止直接将污水排入长江。严格要求工厂将污水经过处理设备处理后，再对各种污染物含量进行检测，检测所要排放的水属于可饮用水，才可排放。若发现工厂将检验未达标水排入长江则处以罚款或者是责令他们拆迁工厂。六、模型评价与推广模型中运用了灰色模型预测法，在预测实践中以中短期预测的效果为佳，弱化了预测对象的随机性变化，突出了趋势性。忽略了次要因素对水质的影响，简洁、易懂，具有一定的合理性和实用性。模型还考虑到了变量剔除问题，正确的选择了变量，减少异常点对模型结果的影响。但是模型对于水质的净化只考虑了自然降解，而对于实际情况比如是否存在污水处理厂等认为净化设施的情况没有考虑。对于降解系数选取，认为全流域相等，与实际情况存在一定的误差。对于污水在监测站之后马上汇入江水的假设，与实际沿江不规律分布的污水排放不符，也会产生误差。实际情况下，污染物排放量与诸多因素有关，仅仅估计其与时间的关系会造成偏差。用给出的前十年的水质数据来预测后十年的水质变化情况，数据不够充足，容易造成误差。七、参考文献 1 邓聚龙1 灰色预测与决策M 1 武汉: 华中理工大学出版社，1986.1 2 姜启源，谢金星，叶俊，数学模型，北京：高等教育出版社，2004.4八、附件第二题程序1clear;clc;a1=2.9 2.8 3.0 3.2 2.95 4.95 2.7 2.3 2.3 1.8 1.5 1.4 1.65;a2=1.5 1.35 1.65 1.85 2.75 3.25 1.9 1.3 1.2 1.0 0.7 0.65 0.8;a3=0.9 0.85 1.25 1.45 1.4 1.8 1.55 0.75 0.75 0.55 0.5 0.5 0.45;a4=0.95 1.0 1.4 1.5 1.45 2.0 1.6 0.85 0.8 0.65 0.65 0.7 0.6;a5=1.05 1.1 1.5 1.55 1.6 2.75 1.65 0.9 0.8 0.7 0.7 0.8 0.75;a6=1.15 1.15 1.55 1.65 1.7 3.4 1.8 0.95 0.85 0.75 0.75 0.85 0.8;b1=950000;b2=778000;b3=395000;b4=500000;b5=164000;b6=464000;c1=b1./(86400.*a1);c2=b2./(86400.*a2);c3=b3./(86400.*a3);c4=b4./(86400.*a4);c5=b5./(86400.*a5);c6=b6./(86400.*a6);c1=vpa(c1,3)c2=vpa(c2,3)c3=vpa(c3,3)c4=vpa(c4,3)c5=vpa(c5,3)c6=vpa(c6,3)运行结果：第二题程序2clear;clc;a1=8487 48300 46200 84480 81490 82600;. 15996 28820 59400 73800 92380 10800;. 10025 42600 77140 79100 75400 51360;. 11184 54120 72640 101220 70200 54230;. 21692 21200 74400 103600 64200 63460;. 38308 209440 181900 209820 138320 118720;. 2608 25920 70670 78050 47120 65100;. 4200 16340 20140 31200 32120 37400;. 1521.6 8515 17760 37450 30360 23970;. 854.4 5226 8683 41769 26160 31980;. 550.8 6485.4 9020 25536 31930 36990;. 685.3 9006 10878 28864 25740 65070;. 706.2 7300 12960 20996 30200 46460;a2=553.5 2898 5760 6656 7587 6490;. 260.4 3537 5940 6560 6556 8500;. 160.4 2556 6902 7910 4425 5457;. 186.4 2952 3632 8676 4320 7018;. 3740 1802 4560 8547 7704 5344;. 565.2 7616 15515 16140 18200 10388;. 260.8 3726 4393 8028 4712 5270;. 90 2205.9 2544 4320 3212 2040;. 76.08 3340.5 1184 3317 1584 1974;. 49.84 2170.8 1005.4 3449.8 2289 2460;. 91.8 1981.65 811.8 3431.4 1339 2466;. 161.98 2559.6 984.2 2956.8 1573 4922;. 64.2 2299.5 702 3040.8 2114 3030;t=3.79 6.00 5.08 6.09 1.81 4.67;. 3.93 6.67 5.38 5.79 1.73 4.67;. 3.67 5.46 3.66 4.13 1.27 3.46;. 3.44 4.87 3.15 3.86 1.22 3.25;. 3.73 3.27 3.27 3.99 1.19 3.16;. 2.22 2.77 2.54 2.89 0.69 1.58;. 4.07 4.74 2.93 3.62 1.15 2.98;. 4.78 6.93 6.10 6.81 2.11 5.65;. 4.78 7.50 6.10 7.23 2.37 6.32;. 6.11 9.00 8.31 8.90 2.71 7.16;. 7.33 12.9 9.14 8.90 2.71 7.16;. 7.85 13.9 9.14 8.27 2.37 6.32;. 6.66 11.3 10.2 9.64 2.59 6.71;x=0.2.t;y1=a1.*x;y2=a2.*x;Q1=vpa(y1,3)Q2=vpa(y2,3)运行结果：Q1 = 19.0, 3.09, 13.0, 4.68, .443e4, 45.0 28.6, .627, 10.3, 6.62, .571e4, 5.88 27.3, 6.50, 213., 103., .977e4, 196. 44.1, 21.3, 456., 203., .985e4, 290. 53.6, 110., 385., 168., .946e4, 392. .108e4, .243e4, .305e4, .200e4, .456e5, .934e4 3.73, 12.6, 633., 230., .740e4, 538. 1.92, .234, 1.10, .542, .108e4, 4.20 .694, .487e-1, .968, .331, 669., .917 .458e-1, .268e-2, .135e-1, .251e-1, 334., .316 .415e-2, .624e-5, .369e-2, .154e-1, 407., .366 .223e-2, .173e-5, .445e-2, .478e-1, 568., 2.49 .156e-1, .922e-4, .962e-3, .384e-2, 467., .948 Q2 = 1.24, .185, 1.62, .369, 412., 3.53 .466, .770e-1, 1.03, .589, 405., 4.63 .436, .390, 19.1, 10.3, 573., 20.8 .734, 1.16, 22.8, 17.4, 606., 37.5 9.24, 9.34, 23.6, 13.9, .113e4, 33.0 15.9, 88.2, 260., 154., .599e4, 817. .373, 1.81, 39.3, 23.7, 740., 43.5 .410e-1, .316e-1, .139, .751e-1, 108., .229 .347e-1, .191e-1, .645e-1, .293e-1, 34.9, .755e-1 .267e-2, .111e-2, .156e-2, .207e-2, 29.2, .243e-1 .691e-3, .191e-5, .332e-3, .206e-2, 17.1, .244e-1 .528e-3, .493e-6, .402e-3, .490e-2, 34.7, .188 .142e-2, .291e-4, .521e-4, .556e-3, 32.7, .619e-1第二题程序3clear;clc;a1=48300 46200 84480 81490 82600 44700;. 28820 59400 73800 92380 10800 62100;. 42600 77140 79100 75400 51360 66200;. 54120 72640 101220 70200 54230 57800;. 21200 74400 103600 64200 63460 59670;. 209440 181900 209820 138320 118720 186300;. 25920 70670 78050 47120 65100 72960;. 16340 20140 31200 32120 37400 35280;. 8515 17760 37450 30360 23970 29800;. 5226 8683 41769 26160 31980 36000;. 6485.4 9020 25536 31930 36990 34730;. 9006 10878 28864 25740 65070 38700;. 7300 12960 20996 30200 46460 33150;a2=2898 5760 6656 7587 6490 596;. 3537 5940 6560 6556 8500 3105;. 2556 6902 7910 4425 5457 1986;. 2952 3632 8676 4320 7018 1605;. 1802 4560 8547 7704 5344 1755;. 7616 15515 16140 18200 10388 6480;. 3726 4393 8028 4712 5270 1536;. 2205.9 2544 4320 3212 2040 980;. 3340.5 1184 3317 1584 1974 1490;. 2170.8 1005.4 3449.8 2289 2460 3168;. 1981.65 811.8 3431.4 1339 2466 4681;. 2559.6 984.2 2956.8 1573 4922 430;. 2299.5 702 3040.8 2114 3030 4420;x1=19.0, 3.09, 13.0, 4.68, .443e4, 45.0;. 28.6, .627, 10.3, 6.62, .571e4, 5.88;. 27.3, 6.50, 213., 103., .977e4, 196.;. 44.1, 21.3, 456., 203., .985e4, 290.;. 53.6, 110., 385., 168., .946e4, 392.;. .108e4, .243e4, .305e4, .200e4, .456e5, .934e4;. 3.73, 12.6, 633., 230., .740e4, 538.;. 1.92, .234, 1.10, .542, .108e4, 4.20;. .694, .487e-1, .968, .331, 669., .917;. .458e-1, .268e-2, .135e-1, .251e-1, 334., .316;. .415e-2, .624e-5, .369e-2, .154e-1, 407., .366;. .223e-2, .173e-5, .445e-2, .478e-1, 568., 2.49;. .156e-1, .922e-4, .962e-3, .384e-2, 467., .948;x2 = 1.24, .185, 1.62, .369, 412., 3.53;. .466, .770e-1, 1.03, .589, 405., 4.63;. .436, .390, 19.1, 10.3, 573., 20.8;. .734, 1.16, 22.8, 17.4, 606., 37.5;. 9.24, 9.34, 23.6, 13.9, .113e4, 33.0;. 15.9, 88.2, 260., 154., .599e4, 817.;. .373, 1.81, 39.3, 23.7, 740., 43.5;. .410e-1, .316e-1, .139, .751e-1, 108., .229;. .347e-1, .191e-1, .645e-1, .293e-1, 34.9, .755e-1;. .267e-2, .111e-2, .156e-2, .207e-2, 29.2, .243e-1;. .691e-3, .191e-5, .332e-3, .206e-2, 17.1, .244e-1;. .528e-3, .493e-6, .402e-3, .490e-2, 34.7, .188;. .142e-2, .291e-4, .521e-4, .556e-3, 3