数值分析-求解线性方程组

一、 习题分析

分析题目可知，要解出A和b涉及到非线性方程组的迭代解法，这里我们考虑用牛顿法和弦割法。牛顿法求解非线性方程组原理如下：

分析题目可知，解出A和b之后，需要解线性方程组Ay=b，这里我们考虑用列主元的高斯消去法和选主元的Doolittle分解法，列主元高斯消去法求解线性方程组原理如下：

选主元的Doolittle分解法求解线性方程组原理如下：

二、 编程思路

1、步骤一：求a1、a2……a999、a1000

非线性方程组为：

其中

对应该方程组在python中可定义为：

牛顿法求解非线性方程组编程思路：

牛顿法在python中可编写为：

其中迭代精度设为：0.000000000001、迭代次数为200次。args为未知量，A为系数矩阵，x_0为未知量初值。

2、步骤二：求解b1、b2……b999、b1000

也用Newton法求解

3、步骤三：知道a1、a2……a999、a1000和b1、b2……b999、b1000后，解线性方程组：Ay=b可得y1、y2……y999、y1000。

列主元高斯消去编程思路：

列主元高斯消去法的python编程：

选主元的Doolittle分解法解线性方程组编程思路：

选主元的Doolittle分解法解线性方程组Python编程：

三、 结果分析

程序运行截图：

列主元高斯消去法解线性方程组部分结果如下，可见y991_{y1000可以满足精度需求，是正确的解，而在回代求解的过程中出现误差积累，y1}y10严重偏离真实值，与真实值大约相差74个数量级。可见，在未知量比较大的情况下，高斯法误差会比较大。

y1=-8.076186e+75 y2=-3.073498e+76 y3=-1.012839e+76 y4=1.652025e+76 y5=6.962022e+75

y6=-1.240855e+76 y7=-1.080654e+76 y8=1.114697e+76 y9=7.532522e+75 y10=-1.612140e+75

………………………………………………………………………………………………………………….

y991=1.866689e-02 y992=9.389536e-02 y993=6.602917e-02 y994=-3.129197e-03 y995=1.030094e-02

y996=2.140479e-02 y997=9.156065e-02 y998=6.851463e-02 y999=-8.960650e-03 y1000=2.085837e-02

为此，我们采用选主元的Doolittle分解法来解该线性方程组

程序运行结果：

a1=5.160498e+00 a2=1.567652e+01 a3=5.302487e+00 a4=1.500328e+01 a5=2.999834e+01

a6=5.165717e+00 a7=1.560918e+01 a8=5.343773e+00 a9=1.500691e+01 a10=2.993440e+01

a11=5.171030e+00 a12=1.554185e+01 a13=5.385009e+00 a14=1.501046e+01 a15=2.987050e+01

a16=5.176437e+00 a17=1.547451e+01 a18=5.426195e+00 a19=1.501395e+01 a20=2.980663e+01

a21=5.181937e+00 a22=1.540717e+01 a23=5.467332e+00 a24=1.501736e+01 a25=2.974280e+01

a26=5.187529e+00 a27=1.533983e+01 a28=5.508417e+00 a29=1.502070e+01 a30=2.967900e+01

a31=5.193215e+00 a32=1.527248e+01 a33=5.549452e+00 a34=1.502397e+01 a35=2.961524e+01

a36=5.198993e+00 a37=1.520514e+01 a38=5.590436e+00 a39=1.502717e+01 a40=2.955152e+01

a41=5.204864e+00 a42=1.513780e+01 a43=5.631369e+00 a44=1.503030e+01 a45=2.948783e+01

a46=5.210827e+00 a47=1.507046e+01 a48=5.672251e+00 a49=1.503335e+01 a50=2.942418e+01

a51=5.216882e+00 a52=1.500312e+01 a53=5.713080e+00 a54=1.503634e+01 a55=2.936058e+01

a56=5.223030e+00 a57=1.493579e+01 a58=5.753858e+00 a59=1.503926e+01 a60=2.929701e+01

a61=5.229269e+00 a62=1.486846e+01 a63=5.794583e+00 a64=1.504211e+01 a65=2.923348e+01

a66=5.235600e+00 a67=1.480113e+01 a68=5.835256e+00 a69=1.504489e+01 a70=2.916999e+01

a71=5.242023e+00 a72=1.473381e+01 a73=5.875875e+00 a74=1.504760e+01 a75=2.910654e+01

a76=5.248538e+00 a77=1.466649e+01 a78=5.916442e+00 a79=1.505024e+01 a80=2.904314e+01

a81=5.255143e+00 a82=1.459917e+01 a83=5.956956e+00 a84=1.505281e+01 a85=2.897978e+01

a86=5.261840e+00 a87=1.453186e+01 a88=5.997415e+00 a89=1.505532e+01 a90=2.891646e+01

a91=5.268628e+00 a92=1.446456e+01 a93=6.037821e+00 a94=1.505775e+01 a95=2.885318e+01

a96=5.275507e+00 a97=1.439727e+01 a98=6.078173e+00 a99=1.506012e+01 a100=2.878994e+01

a101=5.282476e+00 a102=1.432998e+01 a103=6.118470e+00 a104=1.506242e+01 a105=2.872675e+01

a106=5.289536e+00 a107=1.426270e+01 a108=6.158712e+00 a109=1.506465e+01 a110=2.866361e+01

a111=5.296686e+00 a112=1.419543e+01 a113=6.198900e+00 a114=1.506682e+01 a115=2.860051e+01

a116=5.303926e+00 a117=1.412817e+01 a118=6.239032e+00 a119=1.506891e+01 a120=2.853746e+01

a121=5.311256e+00 a122=1.406092e+01 a123=6.279109e+00 a124=1.507094e+01 a125=2.847445e+01

a126=5.318676e+00 a127=1.399368e+01 a128=6.319130e+00 a129=1.507290e+01 a130=2.841149e+01

a131=5.326186e+00 a132=1.392645e+01 a133=6.359096e+00 a134=1.507480e+01 a135=2.834858e+01

a136=5.333785e+00 a137=1.385923e+01 a138=6.399004e+00 a139=1.507663e+01 a140=2.828571e+01

a141=5.341474e+00 a142=1.379202e+01 a143=6.438857e+00 a144=1.507839e+01 a145=2.822290e+01

a146=5.349251e+00 a147=1.372483e+01 a148=6.478652e+00 a149=1.508009e+01 a150=2.816013e+01

a151=5.357118e+00 a152=1.365764e+01 a153=6.518391e+00 a154=1.508172e+01 a155=2.809741e+01

a156=5.365073e+00 a157=1.359048e+01 a158=6.558072e+00 a159=1.508329e+01 a160=2.803475e+01

a161=5.373117e+00 a162=1.352332e+01 a163=6.597696e+00 a164=1.508479e+01 a165=2.797213e+01

a166=5.381249e+00 a167=1.345619e+01 a168=6.637262e+00 a169=1.508622e+01 a170=2.790957e+01

a171=5.389469e+00 a172=1.338906e+01 a173=6.676771e+00 a174=1.508759e+01 a175=2.784705e+01

a176=5.397777e+00 a177=1.332196e+01 a178=6.716220e+00 a179=1.508890e+01 a180=2.778459e+01

a181=5.406174e+00 a182=1.325487e+01 a183=6.755612e+00 a184=1.509014e+01 a185=2.772219e+01

a186=5.414657e+00 a187=1.318779e+01 a188=6.794944e+00 a189=1.509132e+01 a190=2.765983e+01

a191=5.423229e+00 a192=1.312074e+01 a193=6.834218e+00 a194=1.509243e+01 a195=2.759753e+01

a196=5.431887e+00 a197=1.305370e+01 a198=6.873433e+00 a199=1.509349e+01 a200=2.753529e+01

a201=5.440633e+00 a202=1.298668e+01 a203=6.912588e+00 a204=1.509447e+01 a205=2.747310e+01

a206=5.449465e+00 a207=1.291968e+01 a208=6.951683e+00 a209=1.509540e+01 a210=2.741096e+01

a211=5.458385e+00 a212=1.285270e+01 a213=6.990719e+00 a214=1.509626e+01 a215=2.734888e+01

a216=5.467390e+00 a217=1.278574e+01 a218=7.029694e+00 a219=1.509706e+01 a220=2.728686e+01

a221=5.476482e+00 a222=1.271880e+01 a223=7.068609e+00 a224=1.509779e+01 a225=2.722490e+01

a226=5.485661e+00 a227=1.265188e+01 a228=7.107464e+00 a229=1.509847e+01 a230=2.716299e+01

a231=5.494925e+00 a232=1.258498e+01 a233=7.146257e+00 a234=1.509908e+01 a235=2.710114e+01

a236=5.504274e+00 a237=1.251811e+01 a238=7.184990e+00 a239=1.509963e+01 a240=2.703935e+01

a241=5.513710e+00 a242=1.245126e+01 a243=7.223662e+00 a244=1.510012e+01 a245=2.697762e+01

a246=5.523230e+00 a247=1.238443e+01 a248=7.262272e+00 a249=1.510055e+01 a250=2.691595e+01

a251=5.532836e+00 a252=1.231763e+01 a253=7.300820e+00 a254=1.510092e+01 a255=2.685434e+01

a256=5.542527e+00 a257=1.225085e+01 a258=7.339306e+00 a259=1.510123e+01 a260=2.679280e+01

a261=5.552302e+00 a262=1.218409e+01 a263=7.377731e+00 a264=1.510147e+01 a265=2.673131e+01

a266=5.562162e+00 a267=1.211736e+01 a268=7.416093e+00 a269=1.510166e+01 a270=2.666988e+01

a271=5.572107e+00 a272=1.205066e+01 a273=7.454392e+00 a274=1.510179e+01 a275=2.660852e+01

a276=5.582135e+00 a277=1.198398e+01 a278=7.492629e+00 a279=1.510185e+01 a280=2.654722e+01

a281=5.592248e+00 a282=1.191733e+01 a283=7.530803e+00 a284=1.510186e+01 a285=2.648598e+01

a286=5.602444e+00 a287=1.185071e+01 a288=7.568914e+00 a289=1.510181e+01 a290=2.642481e+01

a291=5.612723e+00 a292=1.178411e+01 a293=7.606962e+00 a294=1.510170e+01 a295=2.636370e+01

a296=5.623086e+00 a297=1.171755e+01 a298=7.644946e+00 a299=1.510154e+01 a300=2.630266e+01

a301=5.633532e+00 a302=1.165101e+01 a303=7.682867e+00 a304=1.510131e+01 a305=2.624168e+01

a306=5.644061e+00 a307=1.158450e+01 a308=7.720723e+00 a309=1.510103e+01 a310=2.618076e+01

a311=5.654673e+00 a312=1.151802e+01 a313=7.758516e+00 a314=1.510068e+01 a315=2.611992e+01

a316=5.665367e+00 a317=1.145157e+01 a318=7.796244e+00 a319=1.510029e+01 a320=2.605914e+01

a321=5.676144e+00 a322=1.138515e+01 a323=7.833908e+00 a324=1.509983e+01 a325=2.599843e+01

a326=5.687003e+00 a327=1.131877e+01 a328=7.871508e+00 a329=1.509932e+01 a330=2.593778e+01

a331=5.697943e+00 a332=1.125241e+01 a333=7.909042e+00 a334=1.509875e+01 a335=2.587720e+01

a336=5.708965e+00 a337=1.118609e+01 a338=7.946512e+00 a339=1.509812e+01 a340=2.581670e+01

a341=5.720069e+00 a342=1.111980e+01 a343=7.983917e+00 a344=1.509744e+01 a345=2.575626e+01

a346=5.731254e+00 a347=1.105354e+01 a348=8.021256e+00 a349=1.509670e+01 a350=2.569589e+01

a351=5.742519e+00 a352=1.098732e+01 a353=8.058530e+00 a354=1.509591e+01 a355=2.563559e+01

a356=5.753866e+00 a357=1.092113e+01 a358=8.095738e+00 a359=1.509506e+01 a360=2.557536e+01

a361=5.765293e+00 a362=1.085497e+01 a363=8.132881e+00 a364=1.509416e+01 a365=2.551520e+01

a366=5.776801e+00 a367=1.078885e+01 a368=8.169957e+00 a369=1.509320e+01 a370=2.545512e+01

a371=5.788389e+00 a372=1.072277e+01 a373=8.206968e+00 a374=1.509219e+01 a375=2.539510e+01

a376=5.800057e+00 a377=1.065672e+01 a378=8.243912e+00 a379=1.509112e+01 a380=2.533516e+01

a381=5.811805e+00 a382=1.059071e+01 a383=8.280790e+00 a384=1.509000e+01 a385=2.527529e+01

a386=5.823632e+00 a387=1.052473e+01 a388=8.317601e+00 a389=1.508882e+01 a390=2.521550e+01

a391=5.835538e+00 a392=1.045879e+01 a393=8.354346e+00 a394=1.508760e+01 a395=2.515577e+01

a396=5.847524e+00 a397=1.039289e+01 a398=8.391024e+00 a399=1.508631e+01 a400=2.509613e+01

a401=5.859589e+00 a402=1.032703e+01 a403=8.427635e+00 a404=1.508498e+01 a405=2.503655e+01

a406=5.871732e+00 a407=1.026120e+01 a408=8.464179e+00 a409=1.508359e+01 a410=2.497705e+01

a411=5.883954e+00 a412=1.019541e+01 a413=8.500655e+00 a414=1.508216e+01 a415=2.491763e+01

a416=5.896254e+00 a417=1.012967e+01 a418=8.537064e+00 a419=1.508067e+01 a420=2.485828e+01

a421=5.908632e+00 a422=1.006396e+01 a423=8.573406e+00 a424=1.507912e+01 a425=2.479901e+01

a426=5.921088e+00 a427=9.998289e+00 a428=8.609680e+00 a429=1.507753e+01 a430=2.473981e+01

a431=5.933622e+00 a432=9.932661e+00 a433=8.645887e+00 a434=1.507588e+01 a435=2.468069e+01

a436=5.946233e+00 a437=9.867074e+00 a438=8.682025e+00 a439=1.507419e+01 a440=2.462165e+01

a441=5.958921e+00 a442=9.801528e+00 a443=8.718096e+00 a444=1.507244e+01 a445=2.456268e+01

a446=5.971686e+00 a447=9.736023e+00 a448=8.754099e+00 a449=1.507064e+01 a450=2.450380e+01

a451=5.984528e+00 a452=9.670560e+00 a453=8.790033e+00 a454=1.506879e+01 a455=2.444499e+01

a456=5.997447e+00 a457=9.605140e+00 a458=8.825899e+00 a459=1.506690e+01 a460=2.438626e+01

a461=6.010442e+00 a462=9.539762e+00 a463=8.861697e+00 a464=1.506495e+01 a465=2.432761e+01

a466=6.023513e+00 a467=9.474427e+00 a468=8.897426e+00 a469=1.506295e+01 a470=2.426904e+01

a471=6.036660e+00 a472=9.409135e+00 a473=8.933086e+00 a474=1.506091e+01 a475=2.421055e+01

a476=6.049883e+00 a477=9.343887e+00 a478=8.968678e+00 a479=1.505881e+01 a480=2.415214e+01

a481=6.063181e+00 a482=9.278683e+00 a483=9.004201e+00 a484=1.505667e+01 a485=2.409381e+01

a486=6.076554e+00 a487=9.213524e+00 a488=9.039656e+00 a489=1.505448e+01 a490=2.403556e+01

a491=6.090003e+00 a492=9.148409e+00 a493=9.075041e+00 a494=1.505224e+01 a495=2.397739e+01

a496=6.103527e+00 a497=9.083339e+00 a498=9.110357e+00 a499=1.504995e+01 a500=2.391930e+01

a501=6.117125e+00 a502=9.018315e+00 a503=9.145605e+00 a504=1.504762e+01 a505=2.386130e+01

a506=6.130797e+00 a507=8.953336e+00 a508=9.180783e+00 a509=1.504524e+01 a510=2.380337e+01

a511=6.144544e+00 a512=8.888403e+00 a513=9.215892e+00 a514=1.504281e+01 a515=2.374553e+01

a516=6.158365e+00 a517=8.823517e+00 a518=9.250931e+00 a519=1.504033e+01 a520=2.368777e+01

a521=6.172260e+00 a522=8.758677e+00 a523=9.285902e+00 a524=1.503781e+01 a525=2.363010e+01

a526=6.186228e+00 a527=8.693884e+00 a528=9.320802e+00 a529=1.503524e+01 a530=2.357251e+01

a531=6.200270e+00 a532=8.629138e+00 a533=9.355634e+00 a534=1.503263e+01 a535=2.351500e+01

a536=6.214385e+00 a537=8.564440e+00 a538=9.390396e+00 a539=1.502997e+01 a540=2.345758e+01

a541=6.228573e+00 a542=8.499789e+00 a543=9.425088e+00 a544=1.502726e+01 a545=2.340024e+01

a546=6.242833e+00 a547=8.435187e+00 a548=9.459711e+00 a549=1.502451e+01 a550=2.334298e+01

a551=6.257167e+00 a552=8.370633e+00 a553=9.494265e+00 a554=1.502171e+01 a555=2.328581e+01

a556=6.271572e+00 a557=8.306127e+00 a558=9.528749e+00 a559=1.501887e+01 a560=2.322872e+01

a561=6.286050e+00 a562=8.241670e+00 a563=9.563163e+00 a564=1.501599e+01 a565=2.317172e+01

a566=6.300600e+00 a567=8.177262e+00 a568=9.597507e+00 a569=1.501306e+01 a570=2.311481e+01

a571=6.315222e+00 a572=8.112903e+00 a573=9.631782e+00 a574=1.501009e+01 a575=2.305798e+01

a576=6.329916e+00 a577=8.048593e+00 a578=9.665988e+00 a579=1.500707e+01 a580=2.300123e+01

a581=6.344680e+00 a582=7.984333e+00 a583=9.700124e+00 a584=1.500401e+01 a585=2.294457e+01

a586=6.359516e+00 a587=7.920123e+00 a588=9.734190e+00 a589=1.500091e+01 a590=2.288800e+01

a591=6.374424e+00 a592=7.855963e+00 a593=9.768186e+00 a594=1.499776e+01 a595=2.283152e+01

a596=6.389402e+00 a597=7.791853e+00 a598=9.802114e+00 a599=1.499458e+01 a600=2.277512e+01

a601=6.404450e+00 a602=7.727793e+00 a603=9.835971e+00 a604=1.499134e+01 a605=2.271881e+01

a606=6.419569e+00 a607=7.663784e+00 a608=9.869760e+00 a609=1.498807e+01 a610=2.266258e+01

a611=6.434759e+00 a612=7.599825e+00 a613=9.903479e+00 a614=1.498476e+01 a615=2.260644e+01

a616=6.450019e+00 a617=7.535917e+00 a618=9.937128e+00 a619=1.498140e+01 a620=2.255039e+01

a621=6.465348e+00 a622=7.472059e+00 a623=9.970709e+00 a624=1.497800e+01 a625=2.249443e+01

a626=6.480748e+00 a627=7.408252e+00 a628=1.000422e+01 a629=1.497456e+01 a630=2.243855e+01

a631=6.496217e+00 a632=7.344497e+00 a633=1.003766e+01 a634=1.497108e+01 a635=2.238276e+01

a636=6.511756e+00 a637=7.280792e+00 a638=1.007104e+01 a639=1.496756e+01 a640=2.232706e+01

a641=6.527364e+00 a642=7.217138e+00 a643=1.010434e+01 a644=1.496399e+01 a645=2.227144e+01

a646=6.543042e+00 a647=7.153534e+00 a648=1.013758e+01 a649=1.496039e+01 a650=2.221591e+01

a651=6.558788e+00 a652=7.089982e+00 a653=1.017074e+01 a654=1.495675e+01 a655=2.216047e+01

a656=6.574603e+00 a657=7.026481e+00 a658=1.020384e+01 a659=1.495306e+01 a660=2.210512e+01

a661=6.590488e+00 a662=6.963031e+00 a663=1.023687e+01 a664=1.494934e+01 a665=2.204985e+01

a666=6.606440e+00 a667=6.899631e+00 a668=1.026984e+01 a669=1.494558e+01 a670=2.199467e+01

a671=6.622462e+00 a672=6.836282e+00 a673=1.030273e+01 a674=1.494178e+01 a675=2.193958e+01

a676=6.638552e+00 a677=6.772983e+00 a678=1.033556e+01 a679=1.493793e+01 a680=2.188458e+01

a681=6.654710e+00 a682=6.709735e+00 a683=1.036832e+01 a684=1.493405e+01 a685=2.182966e+01

a686=6.670936e+00 a687=6.646536e+00 a688=1.040102e+01 a689=1.493013e+01 a690=2.177483e+01

a691=6.687231e+00 a692=6.583388e+00 a693=1.043364e+01 a694=1.492618e+01 a695=2.172009e+01

a696=6.703594e+00 a697=6.520290e+00 a698=1.046620e+01 a699=1.492218e+01 a700=2.166543e+01

a701=6.720025e+00 a702=6.457240e+00 a703=1.049870e+01 a704=1.491814e+01 a705=2.161086e+01

a706=6.736523e+00 a707=6.394240e+00 a708=1.053113e+01 a709=1.491407e+01 a710=2.155637e+01

a711=6.753090e+00 a712=6.331288e+00 a713=1.056349e+01 a714=1.490996e+01 a715=2.150197e+01

a716=6.769725e+00 a717=6.268385e+00 a718=1.059579e+01 a719=1.490581e+01 a720=2.144766e+01

a721=6.786427e+00 a722=6.205529e+00 a723=1.062802e+01 a724=1.490162e+01 a725=2.139343e+01

a726=6.803197e+00 a727=6.142721e+00 a728=1.066019e+01 a729=1.489739e+01 a730=2.133928e+01

a731=6.820035e+00 a732=6.079959e+00 a733=1.069230e+01 a734=1.489313e+01 a735=2.128522e+01

a736=6.836941e+00 a737=6.017243e+00 a738=1.072434e+01 a739=1.488883e+01 a740=2.123124e+01

a741=6.853915e+00 a742=5.954573e+00 a743=1.075632e+01 a744=1.488449e+01 a745=2.117734e+01

a746=6.870957e+00 a747=5.891947e+00 a748=1.078824e+01 a749=1.488011e+01 a750=2.112353e+01

a751=6.888066e+00 a752=5.829366e+00 a753=1.082010e+01 a754=1.487570e+01 a755=2.106980e+01

a756=6.905244e+00 a757=5.766827e+00 a758=1.085189e+01 a759=1.487124e+01 a760=2.101615e+01

a761=6.922490e+00 a762=5.704331e+00 a763=1.088363e+01 a764=1.486676e+01 a765=2.096257e+01

a766=6.939805e+00 a767=5.641875e+00 a768=1.091531e+01 a769=1.486223e+01 a770=2.090908e+01

a771=6.957187e+00 a772=5.579460e+00 a773=1.094693e+01 a774=1.485767e+01 a775=2.085566e+01

a776=6.974639e+00 a777=5.517084e+00 a778=1.097849e+01 a779=1.485306e+01 a780=2.080233e+01

a781=6.992159e+00 a782=5.454745e+00 a783=1.101000e+01 a784=1.484843e+01 a785=2.074906e+01

a786=7.009749e+00 a787=5.392442e+00 a788=1.104145e+01 a789=1.484375e+01 a790=2.069587e+01

a791=7.027407e+00 a792=5.330174e+00 a793=1.107284e+01 a794=1.483903e+01 a795=2.064276e+01

a796=7.045136e+00 a797=5.267940e+00 a798=1.110419e+01 a799=1.483428e+01 a800=2.058972e+01

a801=7.062934e+00 a802=5.205737e+00 a803=1.113548e+01 a804=1.482949e+01 a805=2.053674e+01

a806=7.080803e+00 a807=5.143564e+00 a808=1.116672e+01 a809=1.482466e+01 a810=2.048384e+01

a811=7.098742e+00 a812=5.081418e+00 a813=1.119791e+01 a814=1.481980e+01 a815=2.043100e+01

a816=7.116753e+00 a817=5.019298e+00 a818=1.122905e+01 a819=1.481489e+01 a820=2.037823e+01

a821=7.134835e+00 a822=4.957202e+00 a823=1.126015e+01 a824=1.480995e+01 a825=2.032552e+01

a826=7.152990e+00 a827=4.895126e+00 a828=1.129120e+01 a829=1.480497e+01 a830=2.027287e+01

a831=7.171218e+00 a832=4.833069e+00 a833=1.132221e+01 a834=1.479995e+01 a835=2.022029e+01

a836=7.189518e+00 a837=4.771028e+00 a838=1.135317e+01 a839=1.479488e+01 a840=2.016775e+01

a841=7.207894e+00 a842=4.708999e+00 a843=1.138410e+01 a844=1.478978e+01 a845=2.011528e+01

a846=7.226344e+00 a847=4.646979e+00 a848=1.141499e+01 a849=1.478464e+01 a850=2.006285e+01

a851=7.244870e+00 a852=4.584965e+00 a853=1.144585e+01 a854=1.477946e+01 a855=2.001047e+01

a856=7.263472e+00 a857=4.522954e+00 a858=1.147668e+01 a859=1.477423e+01 a860=1.995813e+01

a861=7.282153e+00 a862=4.460941e+00 a863=1.150747e+01 a864=1.476897e+01 a865=1.990584e+01

a866=7.300912e+00 a867=4.398922e+00 a868=1.153824e+01 a869=1.476366e+01 a870=1.985358e+01

a871=7.319752e+00 a872=4.336892e+00 a873=1.156899e+01 a874=1.475831e+01 a875=1.980136e+01

a876=7.338673e+00 a877=4.274846e+00 a878=1.159971e+01 a879=1.475291e+01 a880=1.974916e+01

a881=7.357676e+00 a882=4.212780e+00 a883=1.163042e+01 a884=1.474747e+01 a885=1.969699e+01

a886=7.376764e+00 a887=4.150687e+00 a888=1.166111e+01 a889=1.474198e+01 a890=1.964484e+01

a891=7.395939e+00 a892=4.088561e+00 a893=1.169180e+01 a894=1.473645e+01 a895=1.959270e+01

a896=7.415201e+00 a897=4.026395e+00 a898=1.172248e+01 a899=1.473087e+01 a900=1.954057e+01

a901=7.434552e+00 a902=3.964182e+00 a903=1.175316e+01 a904=1.472524e+01 a905=1.948844e+01

a906=7.453996e+00 a907=3.901915e+00 a908=1.178385e+01 a909=1.471956e+01 a910=1.943631e+01

a911=7.473535e+00 a912=3.839585e+00 a913=1.181454e+01 a914=1.471382e+01 a915=1.938416e+01

a916=7.493170e+00 a917=3.777182e+00 a918=1.184525e+01 a919=1.470804e+01 a920=1.933199e+01

a921=7.512905e+00 a922=3.714697e+00 a923=1.187599e+01 a924=1.470220e+01 a925=1.927979e+01

a926=7.532743e+00 a927=3.652120e+00 a928=1.190675e+01 a929=1.469631e+01 a930=1.922756e+01

a931=7.552688e+00 a932=3.589437e+00 a933=1.193755e+01 a934=1.469035e+01 a935=1.917527e+01

a936=7.572742e+00 a937=3.526637e+00 a938=1.196840e+01 a939=1.468434e+01 a940=1.912292e+01

a941=7.592910e+00 a942=3.463705e+00 a943=1.199930e+01 a944=1.467827e+01 a945=1.907050e+01

a946=7.613197e+00 a947=3.400627e+00 a948=1.203026e+01 a949=1.467213e+01 a950=1.901799e+01

a951=7.633607e+00 a952=3.337385e+00 a953=1.206130e+01 a954=1.466593e+01 a955=1.896538e+01

a956=7.654146e+00 a957=3.273960e+00 a958=1.209242e+01 a959=1.465965e+01 a960=1.891266e+01

a961=7.674819e+00 a962=3.210334e+00 a963=1.212364e+01 a964=1.465331e+01 a965=1.885980e+01

a966=7.695634e+00 a967=3.146482e+00 a968=1.215498e+01 a969=1.464689e+01 a970=1.880678e+01

a971=7.716598e+00 a972=3.082382e+00 a973=1.218644e+01 a974=1.464039e+01 a975=1.875359e+01

a976=7.737718e+00 a977=3.018003e+00 a978=1.221805e+01 a979=1.463381e+01 a980=1.870020e+01

a981=7.759004e+00 a982=2.953317e+00 a983=1.224982e+01 a984=1.462714e+01 a985=1.864658e+01

a986=7.780467e+00 a987=2.888289e+00 a988=1.228178e+01 a989=1.462038e+01 a990=1.859270e+01

a991=7.802117e+00 a992=2.822879e+00 a993=1.231396e+01 a994=1.461353e+01 a995=1.853853e+01

a996=7.823969e+00 a997=2.757045e+00 a998=1.234637e+01 a999=1.460657e+01 a1000=1.848403e+01

b1=1.309800e+00 b2=1.309197e+00 b3=1.308596e+00 b4=1.307997e+00 b5=1.307399e+00

b6=1.306803e+00 b7=1.306208e+00 b8=1.305614e+00 b9=1.305023e+00 b10=1.304432e+00

b11=1.303843e+00 b12=1.303256e+00 b13=1.302670e+00 b14=1.302086e+00 b15=1.301503e+00

b16=1.300922e+00 b17=1.300342e+00 b18=1.299763e+00 b19=1.299186e+00 b20=1.298610e+00

b21=1.298036e+00 b22=1.297464e+00 b23=1.296892e+00 b24=1.296323e+00 b25=1.295754e+00

b26=1.295187e+00 b27=1.294622e+00 b28=1.294057e+00 b29=1.293495e+00 b30=1.292933e+00

b31=1.292373e+00 b32=1.291815e+00 b33=1.291258e+00 b34=1.290702e+00 b35=1.290147e+00

b36=1.289594e+00 b37=1.289043e+00 b38=1.288492e+00 b39=1.287943e+00 b40=1.287396e+00

b41=1.286850e+00 b42=1.286305e+00 b43=1.285761e+00 b44=1.285219e+00 b45=1.284678e+00

b46=1.284139e+00 b47=1.283600e+00 b48=1.283064e+00 b49=1.282528e+00 b50=1.281994e+00

b51=1.281461e+00 b52=1.280929e+00 b53=1.280399e+00 b54=1.279870e+00 b55=1.279342e+00

b56=1.278815e+00 b57=1.278290e+00 b58=1.277766e+00 b59=1.277244e+00 b60=1.276722e+00

b61=1.276202e+00 b62=1.275683e+00 b63=1.275166e+00 b64=1.274650e+00 b65=1.274134e+00

b66=1.273621e+00 b67=1.273108e+00 b68=1.272597e+00 b69=1.272087e+00 b70=1.271578e+00

b71=1.271070e+00 b72=1.270564e+00 b73=1.270058e+00 b74=1.269554e+00 b75=1.269052e+00

b76=1.268550e+00 b77=1.268050e+00 b78=1.267550e+00 b79=1.267052e+00 b80=1.266556e+00

b81=1.266060e+00 b82=1.265566e+00 b83=1.265072e+00 b84=1.264580e+00 b85=1.264089e+00

b86=1.263600e+00 b87=1.263111e+00 b88=1.262624e+00 b89=1.262137e+00 b90=1.261652e+00

b91=1.261168e+00 b92=1.260686e+00 b93=1.260204e+00 b94=1.259724e+00 b95=1.259244e+00

b96=1.258766e+00 b97=1.258289e+00 b98=1.257813e+00 b99=1.257338e+00 b100=1.256864e+00

b101=1.256392e+00 b102=1.255920e+00 b103=1.255450e+00 b104=1.254980e+00 b105=1.254512e+00

b106=1.254045e+00 b107=1.253579e+00 b108=1.253114e+00 b109=1.252650e+00 b110=1.252188e+00

b111=1.251726e+00 b112=1.251265e+00 b113=1.250806e+00 b114=1.250347e+00 b115=1.249890e+00

b116=1.249434e+00 b117=1.248978e+00 b118=1.248524e+00 b119=1.248071e+00 b120=1.247619e+00

b121=1.247168e+00 b122=1.246718e+00 b123=1.246269e+00 b124=1.245821e+00 b125=1.245374e+00

b126=1.244928e+00 b127=1.244484e+00 b128=1.244040e+00 b129=1.243597e+00 b130=1.243155e+00

b131=1.242715e+00 b132=1.242275e+00 b133=1.241836e+00 b134=1.241398e+00 b135=1.240962e+00

b136=1.240526e+00 b137=1.240091e+00 b138=1.239658e+00 b139=1.239225e+00 b140=1.238793e+00

b141=1.238362e+00 b142=1.237933e+00 b143=1.237504e+00 b144=1.237076e+00 b145=1.236649e+00

b146=1.236223e+00 b147=1.235799e+00 b148=1.235375e+00 b149=1.234952e+00 b150=1.234530e+00

b151=1.234109e+00 b152=1.233689e+00 b153=1.233270e+00 b154=1.232851e+00 b155=1.232434e+00

b156=1.232018e+00 b157=1.231603e+00 b158=1.231188e+00 b159=1.230775e+00 b160=1.230362e+00

b161=1.229950e+00 b162=1.229540e+00 b163=1.229130e+00 b164=1.228721e+00 b165=1.228313e+00

b166=1.227906e+00 b167=1.227500e+00 b168=1.227095e+00 b169=1.226691e+00 b170=1.226287e+00

b171=1.225885e+00 b172=1.225483e+00 b173=1.225083e+00 b174=1.224683e+00 b175=1.224284e+00

b176=1.223886e+00 b177=1.223489e+00 b178=1.223093e+00 b179=1.222697e+00 b180=1.222303e+00

b181=1.221909e+00 b182=1.221516e+00 b183=1.221125e+00 b184=1.220734e+00 b185=1.220344e+00

b186=1.219954e+00 b187=1.219566e+00 b188=1.219178e+00 b189=1.218792e+00 b190=1.218406e+00

b191=1.218021e+00 b192=1.217637e+00 b193=1.217254e+00 b194=1.216871e+00 b195=1.216490e+00

b196=1.216109e+00 b197=1.215729e+00 b198=1.215350e+00 b199=1.214972e+00 b200=1.214594e+00

b201=1.214218e+00 b202=1.213842e+00 b203=1.213467e+00 b204=1.213093e+00 b205=1.212719e+00

b206=1.212347e+00 b207=1.211975e+00 b208=1.211604e+00 b209=1.211234e+00 b210=1.210865e+00

b211=1.210497e+00 b212=1.210129e+00 b213=1.209762e+00 b214=1.209396e+00 b215=1.209031e+00

b216=1.208666e+00 b217=1.208303e+00 b218=1.207940e+00 b219=1.207578e+00 b220=1.207216e+00

b221=1.206856e+00 b222=1.206496e+00 b223=1.206137e+00 b224=1.205779e+00 b225=1.205421e+00

b226=1.205065e+00 b227=1.204709e+00 b228=1.204354e+00 b229=1.203999e+00 b230=1.203646e+00

b231=1.203293e+00 b232=1.202941e+00 b233=1.202590e+00 b234=1.202239e+00 b235=1.201889e+00

b236=1.201540e+00 b237=1.201192e+00 b238=1.200844e+00 b239=1.200497e+00 b240=1.200151e+00

b241=1.199806e+00 b242=1.199461e+00 b243=1.199117e+00 b244=1.198774e+00 b245=1.198432e+00

b246=1.198090e+00 b247=1.197749e+00 b248=1.197409e+00 b249=1.197069e+00 b250=1.196731e+00

b251=1.196393e+00 b252=1.196055e+00 b253=1.195719e+00 b254=1.195383e+00 b255=1.195047e+00

b256=1.194713e+00 b257=1.194379e+00 b258=1.194046e+00 b259=1.193714e+00 b260=1.193382e+00

b261=1.193051e+00 b262=1.192721e+00 b263=1.192391e+00 b264=1.192062e+00 b265=1.191734e+00

b266=1.191406e+00 b267=1.191080e+00 b268=1.190753e+00 b269=1.190428e+00 b270=1.190103e+00

b271=1.189779e+00 b272=1.189456e+00 b273=1.189133e+00 b274=1.188811e+00 b275=1.188489e+00

b276=1.188169e+00 b277=1.187849e+00 b278=1.187529e+00 b279=1.187211e+00 b280=1.186892e+00

b281=1.186575e+00 b282=1.186258e+00 b283=1.185942e+00 b284=1.185627e+00 b285=1.185312e+00

b286=1.184998e+00 b287=1.184685e+00 b288=1.184372e+00 b289=1.184060e+00 b290=1.183748e+00

b291=1.183437e+00 b292=1.183127e+00 b293=1.182817e+00 b294=1.182508e+00 b295=1.182200e+00

b296=1.181892e+00 b297=1.181585e+00 b298=1.181279e+00 b299=1.180973e+00 b300=1.180668e+00

b301=1.180364e+00 b302=1.180060e+00 b303=1.179757e+00 b304=1.179454e+00 b305=1.179152e+00

b306=1.178850e+00 b307=1.178550e+00 b308=1.178250e+00 b309=1.177950e+00 b310=1.177651e+00

b311=1.177353e+00 b312=1.177055e+00 b313=1.176758e+00 b314=1.176461e+00 b315=1.176166e+00

b316=1.175870e+00 b317=1.175576e+00 b318=1.175282e+00 b319=1.174988e+00 b320=1.174695e+00

b321=1.174403e+00 b322=1.174111e+00 b323=1.173820e+00 b324=1.173530e+00 b325=1.173240e+00

b326=1.172950e+00 b327=1.172662e+00 b328=1.172373e+00 b329=1.172086e+00 b330=1.171799e+00

b331=1.171513e+00 b332=1.171227e+00 b333=1.170941e+00 b334=1.170657e+00 b335=1.170373e+00

b336=1.170089e+00 b337=1.169806e+00 b338=1.169524e+00 b339=1.169242e+00 b340=1.168961e+00

b341=1.168680e+00 b342=1.168400e+00 b343=1.168120e+00 b344=1.167841e+00 b345=1.167563e+00

b346=1.167285e+00 b347=1.167008e+00 b348=1.166731e+00 b349=1.166455e+00 b350=1.166179e+00

b351=1.165904e+00 b352=1.165630e+00 b353=1.165356e+00 b354=1.165082e+00 b355=1.164809e+00

b356=1.164537e+00 b357=1.164265e+00 b358=1.163994e+00 b359=1.163723e+00 b360=1.163453e+00

b361=1.163183e+00 b362=1.162914e+00 b363=1.162645e+00 b364=1.162377e+00 b365=1.162110e+00

b366=1.161843e+00 b367=1.161576e+00 b368=1.161311e+00 b369=1.161045e+00 b370=1.160780e+00

b371=1.160516e+00 b372=1.160252e+00 b373=1.159989e+00 b374=1.159726e+00 b375=1.159464e+00

b376=1.159202e+00 b377=1.158941e+00 b378=1.158680e+00 b379=1.158420e+00 b380=1.158160e+00

b381=1.157901e+00 b382=1.157642e+00 b383=1.157384e+00 b384=1.157127e+00 b385=1.156869e+00

b386=1.156613e+00 b387=1.156357e+00 b388=1.156101e+00 b389=1.155846e+00 b390=1.155591e+00

b391=1.155337e+00 b392=1.155083e+00 b393=1.154830e+00 b394=1.154578e+00 b395=1.154326e+00

b396=1.154074e+00 b397=1.153823e+00 b398=1.153572e+00 b399=1.153322e+00 b400=1.153072e+00

b401=1.152823e+00 b402=1.152574e+00 b403=1.152326e+00 b404=1.152078e+00 b405=1.151831e+00

b406=1.151584e+00 b407=1.151338e+00 b408=1.151092e+00 b409=1.150847e+00 b410=1.150602e+00

b411=1.150357e+00 b412=1.150114e+00 b413=1.149870e+00 b414=1.149627e+00 b415=1.149385e+00

b416=1.149142e+00 b417=1.148901e+00 b418=1.148660e+00 b419=1.148419e+00 b420=1.148179e+00

b421=1.147939e+00 b422=1.147700e+00 b423=1.147461e+00 b424=1.147223e+00 b425=1.146985e+00

b426=1.146748e+00 b427=1.146511e+00 b428=1.146274e+00 b429=1.146038e+00 b430=1.145803e+00

b431=1.145567e+00 b432=1.145333e+00 b433=1.145098e+00 b434=1.144865e+00 b435=1.144631e+00

b436=1.144398e+00 b437=1.144166e+00 b438=1.143934e+00 b439=1.143702e+00 b440=1.143471e+00

b441=1.143241e+00 b442=1.143010e+00 b443=1.142781e+00 b444=1.142551e+00 b445=1.142322e+00

b446=1.142094e+00 b447=1.141866e+00 b448=1.141638e+00 b449=1.141411e+00 b450=1.141184e+00

b451=1.140958e+00 b452=1.140732e+00 b453=1.140507e+00 b454=1.140282e+00 b455=1.140057e+00

b456=1.139833e+00 b457=1.139609e+00 b458=1.139386e+00 b459=1.139163e+00 b460=1.138940e+00

b461=1.138718e+00 b462=1.138497e+00 b463=1.138275e+00 b464=1.138055e+00 b465=1.137834e+00

b466=1.137614e+00 b467=1.137395e+00 b468=1.137176e+00 b469=1.136957e+00 b470=1.136739e+00

b471=1.136521e+00 b472=1.136303e+00 b473=1.136086e+00 b474=1.135869e+00 b475=1.135653e+00

b476=1.135437e+00 b477=1.135222e+00 b478=1.135007e+00 b479=1.134792e+00 b480=1.134578e+00

b481=1.134364e+00 b482=1.134151e+00 b483=1.133938e+00 b484=1.133725e+00 b485=1.133513e+00

b486=1.133301e+00 b487=1.133089e+00 b488=1.132878e+00 b489=1.132668e+00 b490=1.132457e+00

b491=1.132248e+00 b492=1.132038e+00 b493=1.131829e+00 b494=1.131620e+00 b495=1.131412e+00

b496=1.131204e+00 b497=1.130996e+00 b498=1.130789e+00 b499=1.130583e+00 b500=1.130376e+00

b501=1.130170e+00 b502=1.129965e+00 b503=1.129759e+00 b504=1.129555e+00 b505=1.129350e+00

b506=1.129146e+00 b507=1.128942e+00 b508=1.128739e+00 b509=1.128536e+00 b510=1.128333e+00

b511=1.128131e+00 b512=1.127929e+00 b513=1.127728e+00 b514=1.127527e+00 b515=1.127326e+00

b516=1.127126e+00 b517=1.126926e+00 b518=1.126726e+00 b519=1.126527e+00 b520=1.126328e+00

b521=1.126130e+00 b522=1.125931e+00 b523=1.125734e+00 b524=1.125536e+00 b525=1.125339e+00

b526=1.125142e+00 b527=1.124946e+00 b528=1.124750e+00 b529=1.124555e+00 b530=1.124359e+00

b531=1.124164e+00 b532=1.123970e+00 b533=1.123776e+00 b534=1.123582e+00 b535=1.123388e+00

b536=1.123195e+00 b537=1.123003e+00 b538=1.122810e+00 b539=1.122618e+00 b540=1.122426e+00

b541=1.122235e+00 b542=1.122044e+00 b543=1.121853e+00 b544=1.121663e+00 b545=1.121473e+00

b546=1.121284e+00 b547=1.121094e+00 b548=1.120905e+00 b549=1.120717e+00 b550=1.120529e+00

b551=1.120341e+00 b552=1.120153e+00 b553=1.119966e+00 b554=1.119779e+00 b555=1.119593e+00

b556=1.119406e+00 b557=1.119220e+00 b558=1.119035e+00 b559=1.118850e+00 b560=1.118665e+00

b561=1.118480e+00 b562=1.118296e+00 b563=1.118112e+00 b564=1.117929e+00 b565=1.117746e+00

b566=1.117563e+00 b567=1.117380e+00 b568=1.117198e+00 b569=1.117016e+00 b570=1.116835e+00

b571=1.116654e+00 b572=1.116473e+00 b573=1.116292e+00 b574=1.116112e+00 b575=1.115932e+00

b576=1.115752e+00 b577=1.115573e+00 b578=1.115394e+00 b579=1.115216e+00 b580=1.115037e+00

b581=1.114859e+00 b582=1.114682e+00 b583=1.114504e+00 b584=1.114327e+00 b585=1.114151e+00

b586=1.113974e+00 b587=1.113798e+00 b588=1.113623e+00 b589=1.113447e+00 b590=1.113272e+00

b591=1.113097e+00 b592=1.112923e+00 b593=1.112749e+00 b594=1.112575e+00 b595=1.112401e+00

b596=1.112228e+00 b597=1.112055e+00 b598=1.111882e+00 b599=1.111710e+00 b600=1.111538e+00

b601=1.111366e+00 b602=1.111195e+00 b603=1.111024e+00 b604=1.110853e+00 b605=1.110683e+00

b606=1.110512e+00 b607=1.110343e+00 b608=1.110173e+00 b609=1.110004e+00 b610=1.109835e+00

b611=1.109666e+00 b612=1.109498e+00 b613=1.109330e+00 b614=1.109162e+00 b615=1.108994e+00

b616=1.108827e+00 b617=1.108660e+00 b618=1.108494e+00 b619=1.108327e+00 b620=1.108161e+00

b621=1.107996e+00 b622=1.107830e+00 b623=1.107665e+00 b624=1.107500e+00 b625=1.107336e+00

b626=1.107171e+00 b627=1.107007e+00 b628=1.106844e+00 b629=1.106680e+00 b630=1.106517e+00

b631=1.106355e+00 b632=1.106192e+00 b633=1.106030e+00 b634=1.105868e+00 b635=1.105706e+00

b636=1.105545e+00 b637=1.105384e+00 b638=1.105223e+00 b639=1.105062e+00 b640=1.104902e+00

b641=1.104742e+00 b642=1.104582e+00 b643=1.104423e+00 b644=1.104264e+00 b645=1.104105e+00

b646=1.103946e+00 b647=1.103788e+00 b648=1.103630e+00 b649=1.103472e+00 b650=1.103315e+00

b651=1.103158e+00 b652=1.103001e+00 b653=1.102844e+00 b654=1.102688e+00 b655=1.102532e+00

b656=1.102376e+00 b657=1.102220e+00 b658=1.102065e+00 b659=1.101910e+00 b660=1.101755e+00

b661=1.101601e+00 b662=1.101447e+00 b663=1.101293e+00 b664=1.101139e+00 b665=1.100986e+00

b666=1.100833e+00 b667=1.100680e+00 b668=1.100527e+00 b669=1.100375e+00 b670=1.100223e+00

b671=1.100071e+00 b672=1.099919e+00 b673=1.099768e+00 b674=1.099617e+00 b675=1.099466e+00

b676=1.099316e+00 b677=1.099166e+00 b678=1.099016e+00 b679=1.098866e+00 b680=1.098717e+00

b681=1.098567e+00 b682=1.098418e+00 b683=1.098270e+00 b684=1.098121e+00 b685=1.097973e+00

b686=1.097825e+00 b687=1.097678e+00 b688=1.097530e+00 b689=1.097383e+00 b690=1.097236e+00

b691=1.097090e+00 b692=1.096943e+00 b693=1.096797e+00 b694=1.096651e+00 b695=1.096506e+00

b696=1.096360e+00 b697=1.096215e+00 b698=1.096070e+00 b699=1.095926e+00 b700=1.095781e+00

b701=1.095637e+00 b702=1.095493e+00 b703=1.095349e+00 b704=1.095206e+00 b705=1.095063e+00

b706=1.094920e+00 b707=1.094777e+00 b708=1.094635e+00 b709=1.094493e+00 b710=1.094351e+00

b711=1.094209e+00 b712=1.094068e+00 b713=1.093927e+00 b714=1.093786e+00 b715=1.093645e+00

b716=1.093504e+00 b717=1.093364e+00 b718=1.093224e+00 b719=1.093084e+00 b720=1.092945e+00

b721=1.092806e+00 b722=1.092667e+00 b723=1.092528e+00 b724=1.092389e+00 b725=1.092251e+00

b726=1.092113e+00 b727=1.091975e+00 b728=1.091837e+00 b729=1.091700e+00 b730=1.091563e+00

b731=1.091426e+00 b732=1.091289e+00 b733=1.091153e+00 b734=1.091016e+00 b735=1.090880e+00

b736=1.090745e+00 b737=1.090609e+00 b738=1.090474e+00 b739=1.090339e+00 b740=1.090204e+00

b741=1.090069e+00 b742=1.089935e+00 b743=1.089801e+00 b744=1.089667e+00 b745=1.089533e+00

b746=1.089400e+00 b747=1.089266e+00 b748=1.089133e+00 b749=1.089001e+00 b750=1.088868e+00

b751=1.088736e+00 b752=1.088604e+00 b753=1.088472e+00 b754=1.088340e+00 b755=1.088208e+00

b756=1.088077e+00 b757=1.087946e+00 b758=1.087815e+00 b759=1.087685e+00 b760=1.087554e+00

b761=1.087424e+00 b762=1.087294e+00 b763=1.087165e+00 b764=1.087035e+00 b765=1.086906e+00

b766=1.086777e+00 b767=1.086648e+00 b768=1.086519e+00 b769=1.086391e+00 b770=1.086263e+00

b771=1.086135e+00 b772=1.086007e+00 b773=1.085880e+00 b774=1.085752e+00 b775=1.085625e+00

b776=1.085498e+00 b777=1.085372e+00 b778=1.085245e+00 b779=1.085119e+00 b780=1.084993e+00

b781=1.084867e+00 b782=1.084741e+00 b783=1.084616e+00 b784=1.084491e+00 b785=1.084366e+00

b786=1.084241e+00 b787=1.084116e+00 b788=1.083992e+00 b789=1.083868e+00 b790=1.083744e+00

b791=1.083620e+00 b792=1.083496e+00 b793=1.083373e+00 b794=1.083250e+00 b795=1.083127e+00

b796=1.083004e+00 b797=1.082882e+00 b798=1.082759e+00 b799=1.082637e+00 b800=1.082515e+00

b801=1.082394e+00 b802=1.082272e+00 b803=1.082151e+00 b804=1.082030e+00 b805=1.081909e+00

b806=1.081788e+00 b807=1.081668e+00 b808=1.081547e+00 b809=1.081427e+00 b810=1.081307e+00

b811=1.081188e+00 b812=1.081068e+00 b813=1.080949e+00 b814=1.080830e+00 b815=1.080711e+00

b816=1.080592e+00 b817=1.080474e+00 b818=1.080355e+00 b819=1.080237e+00 b820=1.080119e+00

b821=1.080001e+00 b822=1.079884e+00 b823=1.079767e+00 b824=1.079649e+00 b825=1.079532e+00

b826=1.079416e+00 b827=1.079299e+00 b828=1.079183e+00 b829=1.079066e+00 b830=1.078950e+00

b831=1.078835e+00 b832=1.078719e+00 b833=1.078604e+00 b834=1.078488e+00 b835=1.078373e+00

b836=1.078258e+00 b837=1.078144e+00 b838=1.078029e+00 b839=1.077915e+00 b840=1.077801e+00

b841=1.077687e+00 b842=1.077573e+00 b843=1.077460e+00 b844=1.077346e+00 b845=1.077233e+00

b846=1.077120e+00 b847=1.077007e+00 b848=1.076895e+00 b849=1.076782e+00 b850=1.076670e+00

b851=1.076558e+00 b852=1.076446e+00 b853=1.076334e+00 b854=1.076223e+00 b855=1.076112e+00

b856=1.076000e+00 b857=1.075889e+00 b858=1.075779e+00 b859=1.075668e+00 b860=1.075558e+00

b861=1.075447e+00 b862=1.075337e+00 b863=1.075227e+00 b864=1.075118e+00 b865=1.075008e+00

b866=1.074899e+00 b867=1.074790e+00 b868=1.074681e+00 b869=1.074572e+00 b870=1.074463e+00

b871=1.074355e+00 b872=1.074246e+00 b873=1.074138e+00 b874=1.074030e+00 b875=1.073923e+00

b876=1.073815e+00 b877=1.073708e+00 b878=1.073600e+00 b879=1.073493e+00 b880=1.073387e+00

b881=1.073280e+00 b882=1.073173e+00 b883=1.073067e+00 b884=1.072961e+00 b885=1.072855e+00

b886=1.072749e+00 b887=1.072643e+00 b888=1.072538e+00 b889=1.072432e+00 b890=1.072327e+00

b891=1.072222e+00 b892=1.072117e+00 b893=1.072013e+00 b894=1.071908e+00 b895=1.071804e+00

b896=1.071700e+00 b897=1.071596e+00 b898=1.071492e+00 b899=1.071388e+00 b900=1.071285e+00

b901=1.071182e+00 b902=1.071078e+00 b903=1.070976e+00 b904=1.070873e+00 b905=1.070770e+00

b906=1.070668e+00 b907=1.070565e+00 b908=1.070463e+00 b909=1.070361e+00 b910=1.070259e+00

b911=1.070158e+00 b912=1.070056e+00 b913=1.069955e+00 b914=1.069854e+00 b915=1.069753e+00

b916=1.069652e+00 b917=1.069551e+00 b918=1.069451e+00 b919=1.069350e+00 b920=1.069250e+00

b921=1.069150e+00 b922=1.069050e+00 b923=1.068951e+00 b924=1.068851e+00 b925=1.068752e+00

b926=1.068652e+00 b927=1.068553e+00 b928=1.068454e+00 b929=1.068356e+00 b930=1.068257e+00

b931=1.068159e+00 b932=1.068060e+00 b933=1.067962e+00 b934=1.067864e+00 b935=1.067767e+00

b936=1.067669e+00 b937=1.067571e+00 b938=1.067474e+00 b939=1.067377e+00 b940=1.067280e+00

b941=1.067183e+00 b942=1.067086e+00 b943=1.066990e+00 b944=1.066893e+00 b945=1.066797e+00

b946=1.066701e+00 b947=1.066605e+00 b948=1.066509e+00 b949=1.066414e+00 b950=1.066318e+00

b951=1.066223e+00 b952=1.066128e+00 b953=1.066033e+00 b954=1.065938e+00 b955=1.065843e+00

b956=1.065748e+00 b957=1.065654e+00 b958=1.065560e+00 b959=1.065466e+00 b960=1.065372e+00

b961=1.065278e+00 b962=1.065184e+00 b963=1.065091e+00 b964=1.064997e+00 b965=1.064904e+00

b966=1.064811e+00 b967=1.064718e+00 b968=1.064625e+00 b969=1.064532e+00 b970=1.064440e+00

b971=1.064348e+00 b972=1.064255e+00 b973=1.064163e+00 b974=1.064071e+00 b975=1.063980e+00

b976=1.063888e+00 b977=1.063797e+00 b978=1.063705e+00 b979=1.063614e+00 b980=1.063523e+00

b981=1.063432e+00 b982=1.063341e+00 b983=1.063251e+00 b984=1.063160e+00 b985=1.063070e+00

b986=1.062980e+00 b987=1.062890e+00 b988=1.062800e+00 b989=1.062710e+00 b990=1.062621e+00

b991=1.062531e+00 b992=1.062442e+00 b993=1.062353e+00 b994=1.062264e+00 b995=1.062175e+00

b996=1.062086e+00 b997=1.061997e+00 b998=1.061909e+00 b999=1.061820e+00 b1000=1.061732e+00

正在求解方程，当前进度:99.8%

y1=1.005432e-01 y2=7.909468e-02 y3=-3.127517e-03 y4=2.406533e-02 y5=1.591185e-02

y6=8.372796e-02 y7=6.177224e-02 y8=9.914706e-03 y9=3.535806e-02 y10=1.467703e-02

y11=7.305791e-02 y12=5.578011e-02 y13=2.165020e-02 y14=3.997245e-02 y15=1.226295e-02

y16=6.787281e-02 y17=5.375948e-02 y18=2.779379e-02 y19=4.164610e-02 y20=1.085251e-02

y21=6.555497e-02 y22=5.310211e-02 y23=3.052305e-02 y24=4.213610e-02 y25=1.020055e-02

y26=6.449946e-02 y27=5.290328e-02 y28=3.165895e-02 y29=4.217689e-02 y30=9.927661e-03

y31=6.395239e-02 y32=5.285583e-02 y33=3.213459e-02 y34=4.205485e-02 y35=9.817824e-03

y36=6.359894e-02 y37=5.285776e-02 y38=3.235557e-02 y39=4.187635e-02 y40=9.772739e-03

y41=6.331646e-02 y42=5.287599e-02 y43=3.248347e-02 y44=4.167943e-02 y45=9.752156e-03

y46=6.305942e-02 y47=5.290005e-02 y48=3.257912e-02 y49=4.147724e-02 y50=9.740627e-03

y51=6.281153e-02 y52=5.292674e-02 y53=3.266448e-02 y54=4.127419e-02 y55=9.732489e-03

y56=6.256715e-02 y57=5.295510e-02 y58=3.274720e-02 y59=4.107171e-02 y60=9.725738e-03

y61=6.232436e-02 y62=5.298486e-02 y63=3.282985e-02 y64=4.087021e-02 y65=9.719689e-03

y66=6.208254e-02 y67=5.301595e-02 y68=3.291327e-02 y69=4.066980e-02 y70=9.714111e-03

y71=6.184147e-02 y72=5.304835e-02 y73=3.299770e-02 y74=4.047050e-02 y75=9.708929e-03

y76=6.160107e-02 y77=5.308205e-02 y78=3.308321e-02 y79=4.027228e-02 y80=9.704120e-03

y81=6.136132e-02 y82=5.311706e-02 y83=3.316980e-02 y84=4.007511e-02 y85=9.699674e-03

y86=6.112219e-02 y87=5.315337e-02 y88=3.325747e-02 y89=3.987896e-02 y90=9.695588e-03

y91=6.088367e-02 y92=5.319099e-02 y93=3.334619e-02 y94=3.968380e-02 y95=9.691861e-03

y96=6.064575e-02 y97=5.322990e-02 y98=3.343594e-02 y99=3.948960e-02 y100=9.688491e-03

y101=6.040841e-02 y102=5.327011e-02 y103=3.352671e-02 y104=3.929633e-02 y105=9.685476e-03

y106=6.017164e-02 y107=5.331163e-02 y108=3.361849e-02 y109=3.910395e-02 y110=9.682816e-03

y111=5.993543e-02 y112=5.335445e-02 y113=3.371125e-02 y114=3.891244e-02 y115=9.680507e-03

y116=5.969976e-02 y117=5.339856e-02 y118=3.380497e-02 y119=3.872177e-02 y120=9.678551e-03

y121=5.946463e-02 y122=5.344398e-02 y123=3.389965e-02 y124=3.853192e-02 y125=9.676943e-03

y126=5.923002e-02 y127=5.349070e-02 y128=3.399526e-02 y129=3.834284e-02 y130=9.675683e-03

y131=5.899592e-02 y132=5.353872e-02 y133=3.409180e-02 y134=3.815452e-02 y135=9.674770e-03

y136=5.876232e-02 y137=5.358805e-02 y138=3.418924e-02 y139=3.796692e-02 y140=9.674201e-03

y141=5.852922e-02 y142=5.363868e-02 y143=3.428758e-02 y144=3.778003e-02 y145=9.673974e-03

y146=5.829659e-02 y147=5.369061e-02 y148=3.438680e-02 y149=3.759381e-02 y150=9.674089e-03

y151=5.806444e-02 y152=5.374386e-02 y153=3.448688e-02 y154=3.740823e-02 y155=9.674542e-03

y156=5.783275e-02 y157=5.379841e-02 y158=3.458783e-02 y159=3.722328e-02 y160=9.675333e-03

y161=5.760151e-02 y162=5.385428e-02 y163=3.468961e-02 y164=3.703892e-02 y165=9.676459e-03

y166=5.737071e-02 y167=5.391146e-02 y168=3.479223e-02 y169=3.685513e-02 y170=9.677918e-03

y171=5.714034e-02 y172=5.396995e-02 y173=3.489567e-02 y174=3.667189e-02 y175=9.679709e-03

y176=5.691040e-02 y177=5.402977e-02 y178=3.499991e-02 y179=3.648917e-02 y180=9.681828e-03

y181=5.668088e-02 y182=5.409090e-02 y183=3.510496e-02 y184=3.630694e-02 y185=9.684275e-03

y186=5.645177e-02 y187=5.415336e-02 y188=3.521081e-02 y189=3.612519e-02 y190=9.687048e-03

y191=5.622305e-02 y192=5.421715e-02 y193=3.531743e-02 y194=3.594388e-02 y195=9.690143e-03

y196=5.599473e-02 y197=5.428227e-02 y198=3.542483e-02 y199=3.576300e-02 y200=9.693559e-03

y201=5.576679e-02 y202=5.434873e-02 y203=3.553299e-02 y204=3.558253e-02 y205=9.697293e-03

y206=5.553923e-02 y207=5.441652e-02 y208=3.564190e-02 y209=3.540243e-02 y210=9.701344e-03

y211=5.531204e-02 y212=5.448566e-02 y213=3.575157e-02 y214=3.522269e-02 y215=9.705709e-03

y216=5.508521e-02 y217=5.455614e-02 y218=3.586198e-02 y219=3.504328e-02 y220=9.710386e-03

y221=5.485873e-02 y222=5.462798e-02 y223=3.597312e-02 y224=3.486419e-02 y225=9.715372e-03

y226=5.463260e-02 y227=5.470117e-02 y228=3.608500e-02 y229=3.468538e-02 y230=9.720665e-03

y231=5.440681e-02 y232=5.477572e-02 y233=3.619759e-02 y234=3.450684e-02 y235=9.726263e-03

y236=5.418136e-02 y237=5.485164e-02 y238=3.631090e-02 y239=3.432855e-02 y240=9.732163e-03

y241=5.395623e-02 y242=5.492892e-02 y243=3.642492e-02 y244=3.415049e-02 y245=9.738362e-03

y246=5.373143e-02 y247=5.500759e-02 y248=3.653965e-02 y249=3.397263e-02 y250=9.744859e-03

y251=5.350694e-02 y252=5.508764e-02 y253=3.665508e-02 y254=3.379495e-02 y255=9.751650e-03

y256=5.328276e-02 y257=5.516907e-02 y258=3.677120e-02 y259=3.361743e-02 y260=9.758734e-03

y261=5.305887e-02 y262=5.525190e-02 y263=3.688802e-02 y264=3.344005e-02 y265=9.766106e-03

y266=5.283529e-02 y267=5.533613e-02 y268=3.700553e-02 y269=3.326280e-02 y270=9.773766e-03

y271=5.261199e-02 y272=5.542177e-02 y273=3.712372e-02 y274=3.308564e-02 y275=9.781709e-03

y276=5.238898e-02 y277=5.550882e-02 y278=3.724259e-02 y279=3.290857e-02 y280=9.789934e-03

y281=5.216625e-02 y282=5.559729e-02 y283=3.736215e-02 y284=3.273155e-02 y285=9.798438e-03

y286=5.194378e-02 y287=5.568719e-02 y288=3.748238e-02 y289=3.255458e-02 y290=9.807217e-03

y291=5.172159e-02 y292=5.577852e-02 y293=3.760328e-02 y294=3.237762e-02 y295=9.816269e-03

y296=5.149965e-02 y297=5.587130e-02 y298=3.772486e-02 y299=3.220067e-02 y300=9.825591e-03

y301=5.127797e-02 y302=5.596552e-02 y303=3.784711e-02 y304=3.202369e-02 y305=9.835181e-03

y306=5.105654e-02 y307=5.606121e-02 y308=3.797003e-02 y309=3.184668e-02 y310=9.845035e-03

y311=5.083536e-02 y312=5.615836e-02 y313=3.809362e-02 y314=3.166961e-02 y315=9.855150e-03

y316=5.061441e-02 y317=5.625699e-02 y318=3.821788e-02 y319=3.149246e-02 y320=9.865523e-03

y321=5.039369e-02 y322=5.635710e-02 y323=3.834281e-02 y324=3.131521e-02 y325=9.876152e-03

y326=5.017320e-02 y327=5.645870e-02 y328=3.846840e-02 y329=3.113784e-02 y330=9.887034e-03

y331=4.995294e-02 y332=5.656181e-02 y333=3.859466e-02 y334=3.096034e-02 y335=9.898164e-03

y336=4.973289e-02 y337=5.666642e-02 y338=3.872160e-02 y339=3.078269e-02 y340=9.909541e-03

y341=4.951305e-02 y342=5.677256e-02 y343=3.884920e-02 y344=3.060486e-02 y345=9.921161e-03

y346=4.929342e-02 y347=5.688023e-02 y348=3.897747e-02 y349=3.042684e-02 y350=9.933020e-03

y351=4.907399e-02 y352=5.698944e-02 y353=3.910642e-02 y354=3.024860e-02 y355=9.945116e-03

y356=4.885475e-02 y357=5.710020e-02 y358=3.923604e-02 y359=3.007014e-02 y360=9.957445e-03

y361=4.863571e-02 y362=5.721253e-02 y363=3.936634e-02 y364=2.989142e-02 y365=9.970005e-03

y366=4.841684e-02 y367=5.732643e-02 y368=3.949731e-02 y369=2.971243e-02 y370=9.982791e-03

y371=4.819816e-02 y372=5.744191e-02 y373=3.962897e-02 y374=2.953316e-02 y375=9.995800e-03

y376=4.797965e-02 y377=5.755899e-02 y378=3.976131e-02 y379=2.935358e-02 y380=1.000903e-02

y381=4.776131e-02 y382=5.767767e-02 y383=3.989434e-02 y384=2.917368e-02 y385=1.002247e-02

y386=4.754314e-02 y387=5.779798e-02 y388=4.002805e-02 y389=2.899342e-02 y390=1.003613e-02

y391=4.732512e-02 y392=5.791992e-02 y393=4.016247e-02 y394=2.881281e-02 y395=1.005000e-02

y396=4.710725e-02 y397=5.804350e-02 y398=4.029758e-02 y399=2.863181e-02 y400=1.006408e-02

y401=4.688953e-02 y402=5.816874e-02 y403=4.043339e-02 y404=2.845041e-02 y405=1.007835e-02

y406=4.667195e-02 y407=5.829565e-02 y408=4.056991e-02 y409=2.826859e-02 y410=1.009283e-02

y411=4.645451e-02 y412=5.842425e-02 y413=4.070714e-02 y414=2.808634e-02 y415=1.010750e-02

y416=4.623720e-02 y417=5.855454e-02 y418=4.084508e-02 y419=2.790362e-02 y420=1.012236e-02

y421=4.602001e-02 y422=5.868654e-02 y423=4.098375e-02 y424=2.772042e-02 y425=1.013741e-02

y426=4.580294e-02 y427=5.882027e-02 y428=4.112315e-02 y429=2.753673e-02 y430=1.015264e-02

y431=4.558598e-02 y432=5.895575e-02 y433=4.126328e-02 y434=2.735252e-02 y435=1.016805e-02

y436=4.536914e-02 y437=5.909297e-02 y438=4.140415e-02 y439=2.716777e-02 y440=1.018364e-02

y441=4.515239e-02 y442=5.923197e-02 y443=4.154577e-02 y444=2.698247e-02 y445=1.019941e-02

y446=4.493573e-02 y447=5.937276e-02 y448=4.168814e-02 y449=2.679659e-02 y450=1.021533e-02

y451=4.471917e-02 y452=5.951535e-02 y453=4.183127e-02 y454=2.661012e-02 y455=1.023143e-02

y456=4.450269e-02 y457=5.965976e-02 y458=4.197517e-02 y459=2.642304e-02 y460=1.024768e-02

y461=4.428629e-02 y462=5.980600e-02 y463=4.211985e-02 y464=2.623532e-02 y465=1.026409e-02

y466=4.406995e-02 y467=5.995410e-02 y468=4.226531e-02 y469=2.604695e-02 y470=1.028066e-02

y471=4.385369e-02 y472=6.010407e-02 y473=4.241156e-02 y474=2.585790e-02 y475=1.029737e-02

y476=4.363748e-02 y477=6.025593e-02 y478=4.255861e-02 y479=2.566816e-02 y480=1.031422e-02

y481=4.342132e-02 y482=6.040969e-02 y483=4.270647e-02 y484=2.547771e-02 y485=1.033121e-02

y486=4.320520e-02 y487=6.056538e-02 y488=4.285515e-02 y489=2.528651e-02 y490=1.034834e-02

y491=4.298912e-02 y492=6.072301e-02 y493=4.300465e-02 y494=2.509457e-02 y495=1.036560e-02

y496=4.277308e-02 y497=6.088261e-02 y498=4.315500e-02 y499=2.490185e-02 y500=1.038299e-02

y501=4.255706e-02 y502=6.104419e-02 y503=4.330619e-02 y504=2.470833e-02 y505=1.040049e-02

y506=4.234105e-02 y507=6.120778e-02 y508=4.345825e-02 y509=2.451399e-02 y510=1.041812e-02

y511=4.212506e-02 y512=6.137339e-02 y513=4.361117e-02 y514=2.431881e-02 y515=1.043586e-02

y516=4.190906e-02 y517=6.154104e-02 y518=4.376497e-02 y519=2.412277e-02 y520=1.045370e-02

y521=4.169306e-02 y522=6.171077e-02 y523=4.391966e-02 y524=2.392585e-02 y525=1.047165e-02

y526=4.147705e-02 y527=6.188258e-02 y528=4.407526e-02 y529=2.372802e-02 y530=1.048970e-02

y531=4.126102e-02 y532=6.205650e-02 y533=4.423178e-02 y534=2.352927e-02 y535=1.050784e-02

y536=4.104496e-02 y537=6.223256e-02 y538=4.438922e-02 y539=2.332957e-02 y540=1.052607e-02

y541=4.082886e-02 y542=6.241078e-02 y543=4.454760e-02 y544=2.312889e-02 y545=1.054438e-02

y546=4.061272e-02 y547=6.259118e-02 y548=4.470694e-02 y549=2.292722e-02 y550=1.056278e-02

y551=4.039652e-02 y552=6.277379e-02 y553=4.486725e-02 y554=2.272453e-02 y555=1.058124e-02

y556=4.018027e-02 y557=6.295863e-02 y558=4.502854e-02 y559=2.252080e-02 y560=1.059978e-02

y561=3.996394e-02 y562=6.314573e-02 y563=4.519083e-02 y564=2.231601e-02 y565=1.061838e-02

y566=3.974753e-02 y567=6.333511e-02 y568=4.535412e-02 y569=2.211013e-02 y570=1.063704e-02

y571=3.953103e-02 y572=6.352681e-02 y573=4.551845e-02 y574=2.190313e-02 y575=1.065575e-02

y576=3.931443e-02 y577=6.372084e-02 y578=4.568381e-02 y579=2.169499e-02 y580=1.067451e-02

y581=3.909773e-02 y582=6.391724e-02 y583=4.585024e-02 y584=2.148570e-02 y585=1.069332e-02

y586=3.888091e-02 y587=6.411603e-02 y588=4.601774e-02 y589=2.127521e-02 y590=1.071216e-02

y591=3.866396e-02 y592=6.431725e-02 y593=4.618633e-02 y594=2.106351e-02 y595=1.073103e-02

y596=3.844687e-02 y597=6.452093e-02 y598=4.635602e-02 y599=2.085057e-02 y600=1.074993e-02

y601=3.822964e-02 y602=6.472709e-02 y603=4.652685e-02 y604=2.063637e-02 y605=1.076885e-02

y606=3.801224e-02 y607=6.493577e-02 y608=4.669881e-02 y609=2.042087e-02 y610=1.078778e-02

y611=3.779468e-02 y612=6.514699e-02 y613=4.687194e-02 y614=2.020405e-02 y615=1.080672e-02

y616=3.757693e-02 y617=6.536080e-02 y618=4.704626e-02 y619=1.998588e-02 y620=1.082566e-02

y621=3.735900e-02 y622=6.557722e-02 y623=4.722177e-02 y624=1.976634e-02 y625=1.084460e-02

y626=3.714085e-02 y627=6.579629e-02 y628=4.739850e-02 y629=1.954539e-02 y630=1.086353e-02

y631=3.692250e-02 y632=6.601805e-02 y633=4.757648e-02 y634=1.932301e-02 y635=1.088245e-02

y636=3.670391e-02 y637=6.624253e-02 y638=4.775572e-02 y639=1.909917e-02 y640=1.090134e-02

y641=3.648508e-02 y642=6.646977e-02 y643=4.793624e-02 y644=1.887383e-02 y645=1.092021e-02

y646=3.626600e-02 y647=6.669980e-02 y648=4.811808e-02 y649=1.864697e-02 y650=1.093904e-02

y651=3.604665e-02 y652=6.693267e-02 y653=4.830124e-02 y654=1.841855e-02 y655=1.095783e-02

y656=3.582702e-02 y657=6.716842e-02 y658=4.848575e-02 y659=1.818855e-02 y660=1.097657e-02

y661=3.560710e-02 y662=6.740707e-02 y663=4.867164e-02 y664=1.795692e-02 y665=1.099526e-02

y666=3.538686e-02 y667=6.764869e-02 y668=4.885893e-02 y669=1.772365e-02 y670=1.101388e-02

y671=3.516631e-02 y672=6.789331e-02 y673=4.904765e-02 y674=1.748869e-02 y675=1.103244e-02

y676=3.494541e-02 y677=6.814097e-02 y678=4.923783e-02 y679=1.725200e-02 y680=1.105092e-02

y681=3.472416e-02 y682=6.839172e-02 y683=4.942948e-02 y684=1.701357e-02 y685=1.106931e-02

y686=3.450254e-02 y687=6.864560e-02 y688=4.962264e-02 y689=1.677334e-02 y690=1.108762e-02

y691=3.428053e-02 y692=6.890267e-02 y693=4.981734e-02 y694=1.653128e-02 y695=1.110583e-02

y696=3.405812e-02 y697=6.916297e-02 y698=5.001360e-02 y699=1.628735e-02 y700=1.112394e-02

y701=3.383529e-02 y702=6.942655e-02 y703=5.021145e-02 y704=1.604152e-02 y705=1.114193e-02

y706=3.361202e-02 y707=6.969347e-02 y708=5.041093e-02 y709=1.579375e-02 y710=1.115981e-02

y711=3.338829e-02 y712=6.996377e-02 y713=5.061206e-02 y714=1.554400e-02 y715=1.117755e-02

y716=3.316409e-02 y717=7.023752e-02 y718=5.081488e-02 y719=1.529221e-02 y720=1.119516e-02

y721=3.293939e-02 y722=7.051476e-02 y723=5.101942e-02 y724=1.503836e-02 y725=1.121263e-02

y726=3.271418e-02 y727=7.079555e-02 y728=5.122572e-02 y729=1.478240e-02 y730=1.122994e-02

y731=3.248843e-02 y732=7.107996e-02 y733=5.143381e-02 y734=1.452429e-02 y735=1.124709e-02

y736=3.226213e-02 y737=7.136805e-02 y738=5.164372e-02 y739=1.426397e-02 y740=1.126408e-02

y741=3.203524e-02 y742=7.165987e-02 y743=5.185550e-02 y744=1.400141e-02 y745=1.128088e-02

y746=3.180776e-02 y747=7.195550e-02 y748=5.206918e-02 y749=1.373655e-02 y750=1.129750e-02

y751=3.157965e-02 y752=7.225499e-02 y753=5.228481e-02 y754=1.346935e-02 y755=1.131393e-02

y756=3.135089e-02 y757=7.255843e-02 y758=5.250241e-02 y759=1.319975e-02 y760=1.133015e-02

y761=3.112146e-02 y762=7.286587e-02 y763=5.272204e-02 y764=1.292770e-02 y765=1.134615e-02

y766=3.089133e-02 y767=7.317740e-02 y768=5.294374e-02 y769=1.265314e-02 y770=1.136193e-02

y771=3.066047e-02 y772=7.349309e-02 y773=5.316756e-02 y774=1.237603e-02 y775=1.137748e-02

y776=3.042886e-02 y777=7.381301e-02 y778=5.339353e-02 y779=1.209631e-02 y780=1.139279e-02

y781=3.019646e-02 y782=7.413725e-02 y783=5.362171e-02 y784=1.181391e-02 y785=1.140784e-02

y786=2.996326e-02 y787=7.446590e-02 y788=5.385215e-02 y789=1.152877e-02 y790=1.142263e-02

y791=2.972921e-02 y792=7.479904e-02 y793=5.408490e-02 y794=1.124084e-02 y795=1.143715e-02

y796=2.949428e-02 y797=7.513675e-02 y798=5.432001e-02 y799=1.095004e-02 y800=1.145139e-02

y801=2.925845e-02 y802=7.547914e-02 y803=5.455753e-02 y804=1.065630e-02 y805=1.146532e-02

y806=2.902168e-02 y807=7.582630e-02 y808=5.479753e-02 y809=1.035957e-02 y810=1.147896e-02

y811=2.878393e-02 y812=7.617833e-02 y813=5.504007e-02 y814=1.005976e-02 y815=1.149228e-02

y816=2.854516e-02 y817=7.653533e-02 y818=5.528520e-02 y819=9.756798e-03 y820=1.150526e-02

y821=2.830534e-02 y822=7.689742e-02 y823=5.553299e-02 y824=9.450608e-03 y825=1.151791e-02

y826=2.806443e-02 y827=7.726470e-02 y828=5.578350e-02 y829=9.141109e-03 y830=1.153021e-02

y831=2.782238e-02 y832=7.763730e-02 y833=5.603681e-02 y834=8.828214e-03 y835=1.154215e-02

y836=2.757915e-02 y837=7.801533e-02 y838=5.629299e-02 y839=8.511836e-03 y840=1.155370e-02

y841=2.733469e-02 y842=7.839893e-02 y843=5.655210e-02 y844=8.191883e-03 y845=1.156487e-02

y846=2.708896e-02 y847=7.878823e-02 y848=5.681424e-02 y849=7.868259e-03 y850=1.157564e-02

y851=2.684191e-02 y852=7.918336e-02 y853=5.707948e-02 y854=7.540865e-03 y855=1.158599e-02

y856=2.659348e-02 y857=7.958448e-02 y858=5.734791e-02 y859=7.209595e-03 y860=1.159591e-02

y861=2.634361e-02 y862=7.999173e-02 y863=5.761962e-02 y864=6.874339e-03 y865=1.160539e-02

y866=2.609226e-02 y867=8.040528e-02 y868=5.789470e-02 y869=6.534984e-03 y870=1.161441e-02

y871=2.583935e-02 y872=8.082529e-02 y873=5.817325e-02 y874=6.191408e-03 y875=1.162296e-02

y876=2.558483e-02 y877=8.125194e-02 y878=5.845537e-02 y879=5.843485e-03 y880=1.163103e-02

y881=2.532863e-02 y882=8.168541e-02 y883=5.874117e-02 y884=5.491081e-03 y885=1.163859e-02

y886=2.507067e-02 y887=8.212591e-02 y888=5.903077e-02 y889=5.134058e-03 y890=1.164563e-02

y891=2.481088e-02 y892=8.257363e-02 y893=5.932429e-02 y894=4.772267e-03 y895=1.165214e-02

y896=2.454918e-02 y897=8.302879e-02 y898=5.962185e-02 y899=4.405551e-03 y900=1.165810e-02

y901=2.428548e-02 y902=8.349163e-02 y903=5.992359e-02 y904=4.033746e-03 y905=1.166350e-02

y906=2.401970e-02 y907=8.396238e-02 y908=6.022965e-02 y909=3.656676e-03 y910=1.166830e-02

y911=2.375175e-02 y912=8.444131e-02 y913=6.054019e-02 y914=3.274156e-03 y915=1.167250e-02

y916=2.348150e-02 y917=8.492869e-02 y918=6.085535e-02 y919=2.885989e-03 y920=1.167608e-02

y921=2.320887e-02 y922=8.542481e-02 y923=6.117532e-02 y924=2.491963e-03 y925=1.167902e-02

y926=2.293373e-02 y927=8.592999e-02 y928=6.150028e-02 y929=2.091856e-03 y930=1.168130e-02

y931=2.265596e-02 y932=8.644454e-02 y933=6.183042e-02 y934=1.685425e-03 y935=1.168290e-02

y936=2.237542e-02 y937=8.696884e-02 y938=6.216594e-02 y939=1.272416e-03 y940=1.168379e-02

y941=2.209198e-02 y942=8.750326e-02 y943=6.250709e-02 y944=8.525501e-04 y945=1.168396e-02

y946=2.180547e-02 y947=8.804822e-02 y948=6.285408e-02 y949=4.255314e-04 y950=1.168339e-02

y951=2.151572e-02 y952=8.860415e-02 y953=6.320719e-02 y954=-8.961262e-06 y955=1.168204e-02

y956=2.122256e-02 y957=8.917153e-02 y958=6.356669e-02 y959=-4.512764e-04 y960=1.167991e-02

y961=2.092579e-02 y962=8.975089e-02 y963=6.393289e-02 y964=-9.017821e-04 y965=1.167694e-02

y966=2.062521e-02 y967=9.034276e-02 y968=6.430614e-02 y969=-1.360950e-03 y970=1.167324e-02

y971=2.032035e-02 y972=9.094799e-02 y973=6.468663e-02 y974=-1.828798e-03 y975=1.166792e-02

y976=2.001268e-02 y977=9.156557e-02 y978=6.507595e-02 y979=-2.308831e-03 y980=1.166695e-02

y981=1.968982e-02 y982=9.220801e-02 y983=6.546626e-02 y984=-2.780563e-03 y985=1.162856e-02

y986=1.943551e-02 y987=9.279396e-02 y988=6.591529e-02 y989=-3.390653e-03 y990=1.184279e-02

y991=1.866689e-02 y992=9.389536e-02 y993=6.602917e-02 y994=-3.129197e-03 y995=1.030094e-02

y996=2.140479e-02 y997=9.156065e-02 y998=6.851463e-02 y999=-8.960650e-03 y1000=2.085837e-02

四、 附代码

import sympy

import numpy as np

 

\# 计算向量范数

def Norm2(p):

  sum_of_p = sum([i ** 2 for i in p])

  return sum_of_p ** 0.5

 

\#牛顿法

def newton(A, x_0, args, prec=0.00000000000000001, n=200):

  funcs = sympy.Matrix(A)

  args = args

  jacobianMatrix = funcs.jacobian(args) # 计算雅可比矩阵

  x_0 = sympy.Matrix(x_0)

  x = args

  x_pre = x_0

  for k in range(0, n): # 迭代求解

​    b = - funcs.subs(zip(x, x_pre)) # 计算f（x【k】）

​    J = jacobianMatrix.subs(zip(x, x_pre)) # 计算雅可比矩阵Jf（x【k】）

​    deltx = sympy.Matrix(J.inv() * b) # 使用雅可比矩阵的逆求解deltx

​    x_new = x_pre + deltx # 迭代跟新x

​    if Norm2(funcs.subs(zip(x, x_new))) < prec: # 如果满足精度要求则返回

​      return ([x_new, k])

​    x_pre = x_new

  return ("迭代次数较少，无法满足精度要求")

 

\# 弦割法

def string_cut(A, x_0, args, prec=0.000000000001, n=200):

  funcs = sympy.Matrix(A)

  h = 0.1 # 设置h的大小

  e = np.eye(len(A)) * h

  x_0 = sympy.Matrix(x_0)

  x = args

  x_pre = x_0

  c = []

  for k in range(0, n):

​    for p in range(len(A)):

​      c.append([A[p].subs(zip(x, (x_pre + sympy.Matrix(e[j]))))

​           for j in range(0, len(A))])

​    fij = sympy.Matrix(c) # 计算fij【k】

​    \# pprint.pprint(fij)

​    b = funcs.subs(zip(x, x_pre)) # 计算fi【k】

​    z = fij.inv() * b # 求解z

​    deltx = h * z / (sum(z) - 1) # 计算deltx

​    x_new = x_pre + deltx # 更新x

​    if Norm2(funcs.subs(zip(x, x_new))) < prec: # 满足精度要求则返回

​      return ([x_new, k])

​    x_pre = x_new

​    c = []

  return ("迭代次数较少，无法满足精度要求")

 

 

\#列主元高斯消去法

def gauss(a, b):

  m, n = a.shape # 获取矩阵的行数和列数

  c = np.zeros(n) # 根据矩阵的行数构建一个一维0数组

  for k in range(n - 1):

​    tem = 0.0

​    maxn=0

​    for biaoji in range(k, n):

​      if tem < np.fabs(a[biaoji][k]): # fabs是绝对值，将a中绝对值最大的找出来

​        maxn = biaoji

​    exchange = a[k].copy()

​    a[k] = a[maxn]

​    a[maxn] = exchange # 行交换

​    exchange = b[k].copy()

​    b[k] = b[maxn]

​    b[maxn] = exchange # 行交换

​    for i in range(k+1 , n):

​      c[i] = a[i][k] / a[k][k] # 计算出系数

​      for j in range(k+1, n): # 从K开始，减少不必要的计算

​        a[i][j] = a[i][j] - c[i] * a[k][j] # 对矩阵进行高斯消去

​      b[i] = b[i] - c[i] * b[k]

​    print('\r正在求解方程，当前进度:' +str(k/10)+"%", end="")

  y = np.zeros(n)

  y[n - 1] = b[n - 1] / a[n - 1][n - 1] # 解出x[n-1],为回代作准备

  \# 回代求出方程解

  for i in range(n - 2, -1, -1):

​    sum = 0.0

​    for j in range(n - 1, -1, -1):

​      sum = sum + a[i][j] * y[j]

​    y[i] = (b[i] - sum) / a[i][i]

  \#输出计算结果

  for i in range(n):

​    if (i+1)%5==0:

​      print("y" + str(i + 1) + "="+ "%e" % y[i]) # 输出结果

​    else:

​      print("y" + str(i + 1) + "="+ "%e" % y[i],end=" ") # 输出结果

 

\#选主元的Doolittle法

def Doolittle_pivot(A):

  n = len(A)

  LU = A.copy()

  order1 = np.arange(n)

  for r in range(n):

​    ut = LU[:r,r].reshape(r,1)

​    si = A[r:,r] - np.sum(ut*LU[r:,:r].T,axis=0)

​    ir = np.argmax(np.abs(si))

​    if ir!=0:

​      LU[[r,r+ir],:] = LU[[r+ir,r],:]

​      order1[[r,r+ir]] = order1[[r+ir,r]]

​    lt = LU[r,:r].reshape(r,1)

​    LU[r,r:] = LU[r,r:] - np.sum(lt*LU[:r,r:],axis=0)

​    if r==n-1:

​      continue

​    LU[r+1:,r] = (LU[r+1:,r] - np.sum(ut*LU[r+1:,:r].T,axis=0))/LU[r,r]

  U = np.triu(LU)

  L = np.tril(LU) - np.diag(np.diag(LU)) + np.eye(n)

  order0 = []

  [order0.insert(i,np.where(order1==i)[0][0]) for i in range(n)]

  return LU,L,U

 

\#根据L、U、b解线性方程组

def solveLineq(L,U,b):

  rows = len(b)

  y = np.zeros(rows)

  y[0] = b[0]/L[0,0]

  for k in range(1,rows):

​    y[k] = (b[k] - np.sum(L[k,:k]*y[:k]))/L[k,k]

  x = np.zeros(rows)

  k = rows-1

  x[k] = y[k]/U[k,k]

  for k in range(rows-2,-1,-1):

​    x[k] = (y[k] - np.sum(x[k+1:]*U[k,k+1:]))/U[k,k]

  return x

 

 

if __name__ == "__main__":

  \#-------------------------创建A和b矩阵----------------------------

  A=np.zeros((1000,1000))

  for i in range(0,999):

​    A[i][i+1]=10

​    A[i+1][i]=1

​    if i!=0:

​      A[i+1][i-1]=10

  b=np.zeros((1000,1))

   \# ---------------------------计算xk-------------------------------

  x1, x2, x3 ,x4,x5 = sympy.symbols("x1 x2 x3 x4 x5") # 定义自变量

  args = sympy.Matrix([x1, x2, x3, x4, x5]) # 定义自变量参数矩阵

  \# 定义非线性方程组列表

  x_0 = [2.0, 2.0, 2.01, 2.01, 2.0] # 设置初始向量

  prec = 0.000000000001 # 设置精度

  n = 200 # 设置迭代次数

  for k in range(1,201):

​    tk=1+0.001*(k-1)

​    x_0 = [5.0, 15.0, 5.0, 15.0, 29.0]

​    \# 定义非线性方程组列表

​    A1 = [sympy.exp(-x1)+sympy.exp(-2*x2)+x3-2*x4+tk*x5-5.3,

​       sympy.exp(-2*x1)+sympy.exp(-x2)-2*x3+tk*x4-x5+25.6,

​       tk*x1+3*x2+sympy.exp(-x3)-3*x5+37.8,

​       2*x1+tk*x2+x3-sympy.exp(-x4)+2*sympy.exp(-2*x5)-31.3,

​       x1-2*x2-3*tk*x3+sympy.exp(-2*x4)+3*sympy.exp(-x5)+42.1]

​    p = newton(A1, x_0, args, prec, n) # 牛顿法求解

​    \#p = string_cut(A1, x_0, args, prec, n) # 弦割法求解

​    for i in range(1,6):

​      A[5*(k-1)+i-1][5*(k-1)+i-1]=p[0][i-1]

​      if i==5:

​        print("a" + str(5 * (k - 1) + i)+ "=" + "%e" % p[0][i - 1])

​      else:

​        print("a" + str(5 * (k - 1) + i)+"=" + "%e" % p[0][i - 1],end=" ")

  \# ---------------------------计算bk-------------------------------

  bk = sympy.symbols("bk") # 定义自变量

  args = sympy.Matrix([bk]) # 定义自变量参数矩阵

  bk_0 = [1.3] # 设置初始向量

  prec = 0.000000000001 # 设置精度

  n = 200 # 设置迭代次数

  for k in range(1, 1001):

​    tk = 1 + 0.001 * (k - 1)

​    \# 定义非线性方程组列表

​    A1 = [sympy.exp(-tk * bk) - tk * sympy.ln(bk)]

​    p = newton(A1, bk_0, args, prec, n) # 牛顿法求解

​    \# p = string_cut(A1, bk_0, args, prec, n) # 弦割法求解

​    if k%5==0 :

​      print("b" + str(k)+ "=" + "%e" % p[0][0])

​    else:

​      print("b" + str(k)+ "=" + "%e" % p[0][0],end=" ")

​    b[k-1]=p[0][0]

  x = np.linalg.solve(A, b)

  \# ---------------------------解线性方程组Ay=b-------------------------------

  \#y=gauss(A,b) #列主元高斯消去法解线性方程组

  LU,L,U=Doolittle_pivot(A) #选主元的Doolittle法解线性方程组

  y=solveLineq(L, U, b)

  for i in range(1000):

​    if (i + 1) % 5 == 0:

​      print("y" + str(i + 1) + "=" + "%e" % x[i]) # 输出结果

​    else:

​      print("y" + str(i + 1) + "=" + "%e" % x[i], end=" ") # 输出结果