西交大計(jì)算方法上機(jī)報(bào)告.doc

計(jì)算方法（B）實(shí)驗(yàn)報(bào)告姓名：學(xué)號：學(xué)院：專業(yè)：實(shí)驗(yàn)一 三對角方程組的求解一、 實(shí)驗(yàn)?zāi)康恼莆杖龑欠匠探M求解的方法。二、 實(shí)驗(yàn)內(nèi)容求三對角方程組的解，其中： ， 三、 算法組織設(shè)系數(shù)矩陣為三對角矩陣則方程組稱為三對角方程組。設(shè)矩陣T非奇異，T可分解為T=LU，其中L為下三角矩陣，U為單位上三角矩陣，記可先依次求出中的元素后，令，先求解下三角方程組得出y，再求解上三角方程組。追趕法的算法組織如下：1.輸入三對角矩陣和右端向量；2.將壓縮為四個一維數(shù)組，是T的三對角線性方程組的三個對角，是右端向量。將分解矩陣壓縮為三個一維數(shù)組。3.對做分解（也可以用分解）導(dǎo)出追趕法的計(jì)算步驟如下： 4.回代求解 5. 停止，輸出結(jié)果。四、 MATLAB 程序MATLAB程序見附件1.五、 結(jié)果及分析實(shí)驗(yàn)結(jié)果為：實(shí)驗(yàn)二 Jacobi迭代和Gauss-Seidel迭代解線性方程組一、 實(shí)驗(yàn)?zāi)康恼莆認(rèn)acobi迭代和Gauss-Seidel迭代解線性方程組的方法。二、 實(shí)驗(yàn)內(nèi)容用Jacobi迭代和Gauss-Seidel迭代解電路電流方程組，使各部分電流誤差均小于。三、 算法組織形如的方程組，用Jacobi迭代求解，算法組織如下：1. 將系數(shù)矩陣分解成對角元、下三角部分元和上三角部分元,于是.2. 由。3. 從而構(gòu)成形如迭代格式：其中4. 選取初始向量進(jìn)行迭代計(jì)算。5. 當(dāng)?shù)蟮慕鉂M足題中的約束條件時迭代停止。形如的方程組，用GaussSeidel迭代求解，算法組織如下：1. 將系數(shù)矩陣分解成對角元、下三角部分元和上三角部分元,于是.2. 由。3. 從而構(gòu)成形如迭代格式：其中4. 選取初始向量進(jìn)行迭代計(jì)算。5. 當(dāng)?shù)蟮慕鉂M足題中的約束條件時迭代停止。四、 MATLAB 程序MATLAB程序見附件2，其中1為Jacobi迭代，2為GaussSeidel迭代。五、 結(jié)果及分析Jacobi迭代結(jié)果:方程組的解為迭代次數(shù)GaussSeidel迭代結(jié)果：方程組的解為迭代次數(shù)由以上結(jié)果可知,達(dá)到相同的計(jì)算精度，GaussSeidel迭代比Jacobi迭代的速度快，GaussSeidel迭代比Jacobi迭代次數(shù)少。實(shí)驗(yàn)三多項(xiàng)式插值及誤差計(jì)算一、實(shí)驗(yàn)?zāi)康恼莆斩囗?xiàng)式插值的原理和基本方法。二、實(shí)驗(yàn)內(nèi)容已知，對a. 計(jì)算函數(shù)在點(diǎn)處的值;b. 求插值數(shù)據(jù)點(diǎn)的插值多項(xiàng)式和三次樣條插值多項(xiàng)式;c. 對，計(jì)算和相應(yīng)的函數(shù)值;d. 計(jì)算，解釋所得到結(jié)果。三、 算法組織（一）本題第一問是簡單的用matlab程序可以計(jì)算，算法很簡單。（二） 本題在算法上第二問中的Newton插值多項(xiàng)式和三次樣條插值多項(xiàng)式。計(jì)算兩種插值多項(xiàng)式的算法如下：1. 求Newton插值多項(xiàng)式，算法組織如下：Newton插值多項(xiàng)式的表達(dá)式如下：其中每一項(xiàng)的系數(shù)ci的表達(dá)式如下：根據(jù)上述公式，為了得到系數(shù)需計(jì)算：1) 一階差商2) 二階差商 3) n階差商4) n+1階差商2. 求三次樣條插值多項(xiàng)式，算法組織如下：所謂三次樣條插值多項(xiàng)式是一種對區(qū)間進(jìn)行分段的分段函數(shù)，然后在每一段上進(jìn)行分析，即它在節(jié)點(diǎn)分成的每個小區(qū)間上是3次多項(xiàng)式，其在此區(qū)間上的表達(dá)式如下：因此，只要確定了的值，就確定了整個表達(dá)式，的計(jì)算方法如下：令：則滿足如下n-1個方程：方程中有n+1個未知量，則令和分別為零，則由上面的方程組可得到的值，可得到整個區(qū)間上的三次樣條插值多項(xiàng)式。（三） 第三問和第四問的算法與第二問的算法類似，不再贅述。四、 MATLAB 程序MATLAB程序見附件3。五、 結(jié)果及分析第一問當(dāng)n=5時,各節(jié)點(diǎn)及f(x)值為：x(0)=-1,y(0)=3.846154e-02x(1)=-6.000000e-01,y(1)=1.000000e-01x(2)=-2.000000e-01,y(2)=5.000000e-01x(3)=2.000000e-01,y(3)=5.000000e-01x(4)=6.000000e-01,y(4)=1.000000e-01x(5)=1,y(5)=3.846154e-02當(dāng)n=10時,各節(jié)點(diǎn)及f(x)值為：x(0)=-1,y(0)=3.846154e-02x(1)=-8.000000e-01,y(1)=5.882353e-02x(2)=-6.000000e-01,y(2)=1.000000e-01x(3)=-4.000000e-01,y(3)=2.000000e-01x(4)=-2.000000e-01,y(4)=5.000000e-01x(5)=0,y(5)=1x(6)=2.000000e-01,y(6)=5.000000e-01x(7)=4.000000e-01,y(7)=2.000000e-01x(8)=6.000000e-01,y(8)=1.000000e-01x(9)=8.000000e-01,y(9)=5.882353e-02x(10)=1,y(10)=3.846154e-02當(dāng)n=20時,各節(jié)點(diǎn)及f(x)值為：x(0)=-1,y(0)=3.846154e-02x(1)=-9.000000e-01,y(1)=4.705882e-02x(2)=-8.000000e-01,y(2)=5.882353e-02x(3)=-7.000000e-01,y(3)=7.547170e-02x(4)=-6.000000e-01,y(4)=1.000000e-01x(5)=-5.000000e-01,y(5)=1.379310e-01x(6)=-4.000000e-01,y(6)=2.000000e-01x(7)=-3.000000e-01,y(7)=3.076923e-01x(8)=-2.000000e-01,y(8)=5.000000e-01x(9)=-1.000000e-01,y(9)=8.000000e-01x(10)=0,y(10)=1x(11)=1.000000e-01,y(11)=8.000000e-01x(12)=2.000000e-01,y(12)=5.000000e-01x(13)=3.000000e-01,y(13)=3.076923e-01x(14)=4.000000e-01,y(14)=2.000000e-01x(15)=5.000000e-01,y(15)=1.379310e-01x(16)=6.000000e-01,y(16)=1.000000e-01x(17)=7.000000e-01,y(17)=7.547170e-02x(18)=8.000000e-01,y(18)=5.882353e-02x(19)=9.000000e-01,y(19)=4.705882e-02x(20)=1,y(20)=3.846154e-02第二問牛頓插值算法當(dāng)n=5時Newton插值多項(xiàng)式的系數(shù)分別為：c0=3.846154e-02c1=1.538462e-01c2=1.057692e+00c3=-1.923077e+00c4=1.201923e+00c5=-5.551115e-16當(dāng)n=10時Newton插值多項(xiàng)式的系數(shù)分別為：c0=3.846154e-02c1=1.018100e-01c2=2.601810e-01c3=7.918552e-01c4=2.686652e+00c5=-6.363122e+00c6=-1.767534e+01c7=8.484163e+01c8=-1.679157e+02c9=2.209417e+02c10=-2.209417e+02當(dāng)n=20時Newton插值多項(xiàng)式的系數(shù)分別為：c0=3.846154e-02c1=8.597285e-02c2=1.583710e-01c3=2.860070e-01c4=5.335952e-01c5=1.037751e+00c6=2.001902e+00c7=2.796775e+00c8=-7.543931e+00c9=-1.011991e+02c10=-6.439941e+01c11=2.152780e+03c12=-7.267934e+03c13=1.139374e+04c14=-3.538429e+03c15=-2.830744e+04c16=8.669152e+04c17=-1.592293e+05c18=2.237536e+05c19=-2.601786e+05c20=2.601786e+05三次樣條插值算法當(dāng)n=5時，取邊界條件為自然樣條的三次樣條插值多項(xiàng)式的系數(shù)分別為：m = 0 4.1296 -3.8259 -3.8259 4.1296 0當(dāng)n=10時，取邊界條件為自然樣條的三次樣條插值多項(xiàng)式的系數(shù)分別為：m = 0 0.4101 1.4820 2.4856 18.5755 -46.7878 18.5755 2.4856 1.4820 0.4101 0當(dāng)n=20時，取邊界條件為自然樣條的三次樣條插值多項(xiàng)式的系數(shù)分別為：m = 0 0.3615 0.4543 0.7514 1.2681 2.2177 4.3438 7.7810 15.3016 -4.3719 -57.8141 -4.3719 15.3016 7.7810 4.3438 2.2177 1.2681 0.7514 0.4543 0.3615 0第三問當(dāng)n=5時，給定點(diǎn)xi的f(xi),N(xi),S(xi),分別為:x(1)=-9.800000e-01,f=3.998401e-02,N=1.369250e-02,S=3.604615e-02x(2)=-9.600000e-01,f=4.159734e-02,N=-6.920000e-03,S=3.371336e-02x(3)=-9.400000e-01,f=4.330879e-02,N=-2.359981e-02,S=3.154575e-02x(4)=-9.200000e-01,f=4.512635e-02,N=-3.656615e-02,S=2.962591e-02x(5)=-9.000000e-01,f=4.705882e-02,N=-4.603365e-02,S=2.803644e-02x(6)=-8.800000e-01,f=4.911591e-02,N=-5.221231e-02,S=2.685992e-02x(7)=-8.600000e-01,f=5.130836e-02,N=-5.530750e-02,S=2.617895e-02x(8)=-8.400000e-01,f=5.364807e-02,N=-5.552000e-02,S=2.607611e-02x(9)=-8.200000e-01,f=5.614823e-02,N=-5.304596e-02,S=2.663401e-02x(10)=-8.000000e-01,f=5.882353e-02,N=-4.807692e-02,S=2.793522e-02x(11)=-7.800000e-01,f=6.169031e-02,N=-4.079981e-02,S=3.006235e-02x(12)=-7.600000e-01,f=6.476684e-02,N=-3.139692e-02,S=3.309798e-02x(13)=-7.400000e-01,f=6.807352e-02,N=-2.004596e-02,S=3.712470e-02x(14)=-7.200000e-01,f=7.163324e-02,N=-6.920000e-03,S=4.222510e-02x(15)=-7.000000e-01,f=7.547170e-02,N=7.812500e-03,S=4.848178e-02x(16)=-6.800000e-01,f=7.961783e-02,N=2.398769e-02,S=5.597733e-02x(17)=-6.600000e-01,f=8.410429e-02,N=4.144635e-02,S=6.479433e-02x(18)=-6.400000e-01,f=8.896797e-02,N=6.003385e-02,S=7.501538e-02x(19)=-6.200000e-01,f=9.425071e-02,N=7.960019e-02,S=8.672308e-02x(20)=-6.000000e-01,f=1.000000e-01,N=1.000000e-01,S=1.000000e-01x(21)=-5.800000e-01,f=1.062699e-01,N=1.210925e-01,S=1.148885e-01x(22)=-5.600000e-01,f=1.131222e-01,N=1.427415e-01,S=1.312696e-01x(23)=-5.400000e-01,f=1.206273e-01,N=1.648156e-01,S=1.489844e-01x(24)=-5.200000e-01,f=1.288660e-01,N=1.871877e-01,S=1.678737e-01x(25)=-5.000000e-01,f=1.379310e-01,N=2.097356e-01,S=1.877783e-01x(26)=-4.800000e-01,f=1.479290e-01,N=2.323415e-01,S=2.085393e-01x(27)=-4.600000e-01,f=1.589825e-01,N=2.548925e-01,S=2.299974e-01x(28)=-4.400000e-01,f=1.712329e-01,N=2.772800e-01,S=2.519935e-01x(29)=-4.200000e-01,f=1.848429e-01,N=2.994002e-01,S=2.743686e-01x(30)=-4.000000e-01,f=2.000000e-01,N=3.211538e-01,S=2.969636e-01x(31)=-3.800000e-01,f=2.169197e-01,N=3.424463e-01,S=3.196192e-01x(32)=-3.600000e-01,f=2.358491e-01,N=3.631877e-01,S=3.421765e-01x(33)=-3.400000e-01,f=2.570694e-01,N=3.832925e-01,S=3.644763e-01x(34)=-3.200000e-01,f=2.808989e-01,N=4.026800e-01,S=3.863595e-01x(35)=-3.000000e-01,f=3.076923e-01,N=4.212740e-01,S=4.076670e-01x(36)=-2.800000e-01,f=3.378378e-01,N=4.390031e-01,S=4.282397e-01x(37)=-2.600000e-01,f=3.717472e-01,N=4.558002e-01,S=4.479184e-01x(38)=-2.400000e-01,f=4.098361e-01,N=4.716031e-01,S=4.665441e-01x(39)=-2.200000e-01,f=4.524887e-01,N=4.863540e-01,S=4.839577e-01x(40)=-2.000000e-01,f=5.000000e-01,N=5.000000e-01,S=5.000000e-01x(41)=-1.800000e-01,f=5.524862e-01,N=5.124925e-01,S=5.145385e-01x(42)=-1.600000e-01,f=6.097561e-01,N=5.237877e-01,S=5.275466e-01x(43)=-1.400000e-01,f=6.711409e-01,N=5.338463e-01,S=5.390243e-01x(44)=-1.200000e-01,f=7.352941e-01,N=5.426338e-01,S=5.489717e-01x(45)=-1.000000e-01,f=8.000000e-01,N=5.501202e-01,S=5.573887e-01x(46)=-8.000000e-02,f=8.620690e-01,N=5.562800e-01,S=5.642753e-01x(47)=-6.000000e-02,f=9.174312e-01,N=5.610925e-01,S=5.696316e-01x(48)=-4.000000e-02,f=9.615385e-01,N=5.645415e-01,S=5.734575e-01x(49)=-2.000000e-02,f=9.900990e-01,N=5.666156e-01,S=5.757530e-01x(50)=0,f=1,N=5.673077e-01,S=5.765182e-01x(51)=2.000000e-02,f=9.900990e-01,N=5.666156e-01,S=5.757530e-01x(52)=4.000000e-02,f=9.615385e-01,N=5.645415e-01,S=5.734575e-01x(53)=6.000000e-02,f=9.174312e-01,N=5.610925e-01,S=5.696316e-01x(54)=8.000000e-02,f=8.620690e-01,N=5.562800e-01,S=5.642753e-01x(55)=1.000000e-01,f=8.000000e-01,N=5.501202e-01,S=5.573887e-01x(56)=1.200000e-01,f=7.352941e-01,N=5.426338e-01,S=5.489717e-01x(57)=1.400000e-01,f=6.711409e-01,N=5.338463e-01,S=5.390243e-01x(58)=1.600000e-01,f=6.097561e-01,N=5.237877e-01,S=5.275466e-01x(59)=1.800000e-01,f=5.524862e-01,N=5.124925e-01,S=5.145385e-01x(60)=2.000000e-01,f=5.000000e-01,N=5.000000e-01,S=5.000000e-01x(61)=2.200000e-01,f=4.524887e-01,N=4.863540e-01,S=4.839577e-01x(62)=2.400000e-01,f=4.098361e-01,N=4.716031e-01,S=4.665441e-01x(63)=2.600000e-01,f=3.717472e-01,N=4.558002e-01,S=4.479184e-01x(64)=2.800000e-01,f=3.378378e-01,N=4.390031e-01,S=4.282397e-01x(65)=3.000000e-01,f=3.076923e-01,N=4.212740e-01,S=4.076670e-01x(66)=3.200000e-01,f=2.808989e-01,N=4.026800e-01,S=3.863595e-01x(67)=3.400000e-01,f=2.570694e-01,N=3.832925e-01,S=3.644763e-01x(68)=3.600000e-01,f=2.358491e-01,N=3.631877e-01,S=3.421765e-01x(69)=3.800000e-01,f=2.169197e-01,N=3.424463e-01,S=3.196192e-01x(70)=4.000000e-01,f=2.000000e-01,N=3.211538e-01,S=2.969636e-01x(71)=4.200000e-01,f=1.848429e-01,N=2.994002e-01,S=2.743686e-01x(72)=4.400000e-01,f=1.712329e-01,N=2.772800e-01,S=2.519935e-01x(73)=4.600000e-01,f=1.589825e-01,N=2.548925e-01,S=2.299974e-01x(74)=4.800000e-01,f=1.479290e-01,N=2.323415e-01,S=2.085393e-01x(75)=5.000000e-01,f=1.379310e-01,N=2.097356e-01,S=1.877783e-01x(76)=5.200000e-01,f=1.288660e-01,N=1.871877e-01,S=1.678737e-01x(77)=5.400000e-01,f=1.206273e-01,N=1.648156e-01,S=1.489844e-01x(78)=5.600000e-01,f=1.131222e-01,N=1.427415e-01,S=1.312696e-01x(79)=5.800000e-01,f=1.062699e-01,N=1.210925e-01,S=1.148885e-01x(80)=6.000000e-01,f=1.000000e-01,N=1.000000e-01,S=1.000000e-01x(81)=6.200000e-01,f=9.425071e-02,N=7.960019e-02,S=8.672308e-02x(82)=6.400000e-01,f=8.896797e-02,N=6.003385e-02,S=7.501538e-02x(83)=6.600000e-01,f=8.410429e-02,N=4.144635e-02,S=6.479433e-02x(84)=6.800000e-01,f=7.961783e-02,N=2.398769e-02,S=5.597733e-02x(85)=7.000000e-01,f=7.547170e-02,N=7.812500e-03,S=4.848178e-02x(86)=7.200000e-01,f=7.163324e-02,N=-6.920000e-03,S=4.222510e-02x(87)=7.400000e-01,f=6.807352e-02,N=-2.004596e-02,S=3.712470e-02x(88)=7.600000e-01,f=6.476684e-02,N=-3.139692e-02,S=3.309798e-02x(89)=7.800000e-01,f=6.169031e-02,N=-4.079981e-02,S=3.006235e-02x(90)=8.000000e-01,f=5.882353e-02,N=-4.807692e-02,S=2.793522e-02x(91)=8.200000e-01,f=5.614823e-02,N=-5.304596e-02,S=2.663401e-02x(92)=8.400000e-01,f=5.364807e-02,N=-5.552000e-02,S=2.607611e-02x(93)=8.600000e-01,f=5.130836e-02,N=-5.530750e-02,S=2.617895e-02x(94)=8.800000e-01,f=4.911591e-02,N=-5.221231e-02,S=2.685992e-02x(95)=9.000000e-01,f=4.705882e-02,N=-4.603365e-02,S=2.803644e-02x(96)=9.200000e-01,f=4.512635e-02,N=-3.656615e-02,S=2.962591e-02x(97)=9.400000e-01,f=4.330879e-02,N=-2.359981e-02,S=3.154575e-02x(98)=9.600000e-01,f=4.159734e-02,N=-6.920000e-03,S=3.371336e-02x(99)=9.800000e-01,f=3.998401e-02,N=1.369250e-02,S=3.604615e-02當(dāng)n=10時，給定點(diǎn)xi的f(xi),N(xi),S(xi),分別為:x(1)=-9.800000e-01,f=3.998401e-02,N=1.230317e+00,S=4.022710e-02x(2)=-9.600000e-01,f=4.159734e-02,N=1.804385e+00,S=4.200907e-02x(3)=-9.400000e-01,f=4.330879e-02,N=1.958952e+00,S=4.382384e-02x(4)=-9.200000e-01,f=4.512635e-02,N=1.845845e+00,S=4.568782e-02x(5)=-9.000000e-01,f=4.705882e-02,N=1.578721e+00,S=4.761740e-02x(6)=-8.800000e-01,f=4.911591e-02,N=1.240213e+00,S=4.962900e-02x(7)=-8.600000e-01,f=5.130836e-02,N=8.880811e-01,S=5.173901e-02x(8)=-8.400000e-01,f=5.364807e-02,N=5.604444e-01,S=5.396383e-02x(9)=-8.200000e-01,f=5.614823e-02,N=2.802176e-01,S=5.631987e-02x(10)=-8.000000e-01,f=5.882353e-02,N=5.882353e-02,S=5.882353e-02x(11)=-7.800000e-01,f=6.169031e-02,N=-1.007243e-01,S=6.149562e-02x(12)=-7.600000e-01,f=6.476684e-02,N=-2.012964e-01,S=6.437461e-02x(13)=-7.400000e-01,f=6.807352e-02,N=-2.496082e-01,S=6.750337e-02x(14)=-7.200000e-01,f=7.163324e-02,N=-2.546027e-01,S=7.092479e-02x(15)=-7.000000e-01,f=7.547170e-02,N=-2.261963e-01,S=7.468173e-02x(16)=-6.800000e-01,f=7.961783e-02,N=-1.743355e-01,S=7.881707e-02x(17)=-6.600000e-01,f=8.410429e-02,N=-1.083152e-01,S=8.337369e-02x(18)=-6.400000e-01,f=8.896797e-02,N=-3.631685e-02,S=8.839447e-02x(19)=-6.200000e-01,f=9.425071e-02,N=3.487332e-02,S=9.392228e-02x(20)=-6.000000e-01,f=1.000000e-01,N=1.000000e-01,S=1.000000e-01x(21)=-5.800000e-01,f=1.062699e-01,N=1.553586e-01,S=1.066700e-01x(22)=-5.600000e-01,f=1.131222e-01,N=1.987262e-01,S=1.139730e-01x(23)=-5.400000e-01,f=1.206273e-01,N=2.292273e-01,S=1.219491e-01x(24)=-5.200000e-01,f=1.288660e-01,N=2.471535e-01,S=1.306384e-01x(25)=-5.000000e-01,f=1.379310e-01,N=2.537555e-01,S=1.400810e-01x(26)=-4.800000e-01,f=1.479290e-01,N=2.510206e-01,S=1.503172e-01x(27)=-4.600000e-01,f=1.589825e-01,N=2.414494e-01,S=1.613870e-01x(28)=-4.400000e-01,f=1.712329e-01,N=2.278397e-01,S=1.733307e-01x(29)=-4.200000e-01,f=1.848429e-01,N=2.130858e-01,S=1.861883e-01x(30)=-4.000000e-01,f=2.000000e-01,N=2.000000e-01,S=2.000000e-01x(31)=-3.800000e-01,f=2.169197e-01,N=1.911589e-01,S=2.149065e-01x(32)=-3.600000e-01,f=2.358491e-01,N=1.887784e-01,S=2.314509e-01x(33)=-3.400000e-01,f=2.570694e-01,N=1.946178e-01,S=2.502767e-01x(34)=-3.200000e-01,f=2.808989e-01,N=2.099145e-01,S=2.720276e-01x(35)=-3.000000e-01,f=3.076923e-01,N=2.353466e-01,S=2.973471e-01x(36)=-2.800000e-01,f=3.378378e-01,N=2.710241e-01,S=3.268788e-01x(37)=-2.600000e-01,f=3.717472e-01,N=3.165048e-01,S=3.612664e-01x(38)=-2.400000e-01,f=4.098361e-01,N=3.708328e-01,S=4.011534e-01x(39)=-2.200000e-01,f=4.524887e-01,N=4.325960e-01,S=4.471834e-01x(40)=-2.000000e-01,f=5.000000e-01,N=5.000000e-01,S=5.000000e-01x(41)=-1.800000e-01,f=5.524862e-01,N=5.709536e-01,S=5.597038e-01x(42)=-1.600000e-01,f=6.097561e-01,N=6.431626e-01,S=6.242233e-01x(43)=-1.400000e-01,f=6.711409e-01,N=7.142283e-01,S=6.909440e-01x(44)=-1.200000e-01,f=7.352941e-01,N=7.817471e-01,S=7.572512e-01x(45)=-1.000000e-01,f=8.000000e-01,N=8.434074e-01,S=8.205306e-01x(46)=-8.000000e-02,f=8.620690e-01,N=8.970803e-01,S=8.781675e-01x(47)=-6.000000e-02,f=9.174312e-01,N=9.409023e-01,S=9.275474e-01x(48)=-4.000000e-02,f=9.615385e-01,N=9.733459e-01,S=9.660558e-01x(49)=-2.000000e-02,f=9.900990e-01,N=9.932776e-01,S=9.910782e-01x(50)=0,f=1,N=1.000000e+00,S=1x(51)=2.000000e-02,f=9.900990e-01,N=9.932776e-01,S=9.910782e-01x(52)=4.000000e-02,f=9.615385e-01,N=9.733459e-01,S=9.660558e-01x(53)=6.000000e-02,f=9.174312e-01,N=9.409023e-01,S=9.275474e-01x(54)=8.000000e-02,f=8.620690e-01,N=8.970803e-01,S=8.781675e-01x(55)=1.000000e-01,f=8.000000e-01,N=8.434074e-01,S=8.205306e-01x(56)=1.200000e-01,f=7.352941e-01,N=7.817471e-01,S=7.572512e-01x(57)=1.400000e-01,f=6.711409e-01,N=7.142283e-01,S=6.909440e-01x(58)=1.600000e-01,f=6.097561e-01,N=6.431626e-01,S=6.242233e-01x(59)=1.800000e-01,f=5.524862e-01,N=5.709536e-01,S=5.597038e-01x(60)=2.000000e-01,f=5.000000e-01,N=5.000000e-01,S=5.000000e-01x(61)=2.200000e-01,f=4.524887e-01,N=4.325960e-01,S=4.471834e-01x(62)=2.400000e-01,f=4.098361e-01,N=3.708328e-01,S=4.011534e-01x(63)=2.600000e-01,f=3.717472e-01,N=3.165048e-01,S=3.612664e-01x(64)=2.800000e-01,f=3.378378e-01,N=2.710241e-01,S=3.268788e-01x(65)=3.000000e-01,f=3.076923e-01,N=2.353466e-01,S=2.973471e-01x(66)=3.200000e-01,f=2.808989e-01,N=2.099145e-01,S=2.720276e-01x(67)=3.400000e-01,f=2.570694e-01,N=1.946178e-01,S=2.502767e-01x(68)=3.600000e-01,f=2.358491e-01,N=1.887784e-01,S=2.314509e-01x(69)=3.800000e-01,f=2.169197e-01,N=1.911589e-01,S=2.149065e-01x(70)=4.000000e-01,f=2.000000e-01,N=2.000000e-01,S=2.000000e-01x(71)=4.200000e-01,f=1.848429e-01,N=2.130858e-01,S=1.861883e-01x(72)=4.400000e-01,f=1.712329e-01,N=2.278397e-01,S=1.733307e-01x(73)=4.600000e-01,f=1.589825e-01,N=2.414494e-01,S=1.613870e-01x(74)=4.800000e-01,f=1.479290e-01,N=2.510206e-01,S=1.503172e-01x(75)=5.000000e-01,f=1.379310e-01,N=2.537555e-01,S=1.400810e-01x(76)=5.200000e-01,f=1.288660e-01,N=2.471535e-01,S=1.306384e-01x(77)=5.400000e-01,f=1.206273e-01,N=2.292273e-01,S=1.219491e-01x(78)=5.600000e-01,f=1.131222e-01,N=1.987262e-01,S=1.139730e-01x(79)=5.800000e-01,f=1.062699e-01,N=1.553586e-01,S=1.066700e-01x(80)=6.000000e-01,f=1.000000e-01,N=1.000000e-01,S=1.000000e-01x(81)=6.200000e-01,f=9.425071e-02,N=3.487332e-02,S=9.392228e-02x(82)=6.400000e-01,f=8.896797e-02,N=-3.631685e-02,S=8.839447e-02x(83)=6.600000e-01,f=8.410429e-02,N=-1.083152e-01,S=8.337369e-02x(84)=6.800000e-01,f=7.961783e-02,N=-1.743355e-01,S=7.881707e-02x(85)=7.000000e-01,f=7.547170e-02,N=-2.261963e-01,S=7.468173e-02x(86)=7.200000e-01,f=7.163324e-02,N=-2.546027e-01,S=7.092479e-02x(87)=7.400000e-01,f=6.807352e-02,N=-2.496082e-01,S=6.750337e-02x(88)=7.600000e-01,f=6.476684e-02,N=-2.012964e-01,S=6.437461e-02x(89)=7.800000e-01,f=6.169031e-02,N=-1.007243e-01,S=6.149562e-02x(90)=8.000000e-01,f=5.882353e-02,N=5.882353e-02,S=5.882353e-02x(91)=8.200000e-01,f=5.614823e-02,N=2.802176e-01,S=5.631987e-02x(92)=8.400000e-01,f=5.364807e-02,N=5.604444e-01,S=5.396383e-02x(93)=8.600000e-01,f=5.130836e-02,N=8.880811e-01,S=5.173901e-02x(94)=8.800000e-01,f=4.911591e-02,N=1.240213e+00,S=4.962900e-02x(95)=9.000000e-01,f=4.705882e-02,N=1.578721e+00,S=4.761740e-02x(96)=9.2000