【正文】 S()= S()=S()= S()= S()=S()=習(xí)題五 重積分的計算23（上機題）重積分的計算題目：給定積分。=39。[39。k4。(j+1)39。到39。39。 s[2]=M[j]/2。} //輸出插值多項式的系數(shù) for(j=0。]39。M39。i) { M[i]=(yy[i]v[i]*M[i+1])/t[i]。 yy[i]=d[i]l[i]*yy[i1]。 //消元過程 for(i=1。i++) {b[i]=2。\t39。]39。d39。 } coutendl。39。39。 for(j=0。=39。[39。j=n。j++) {d[j]=6*((y[j+1]y[j])/h[j](y[j]y[j1])/h[j1])/(h[j]+h[j1])。j++) { u[j]=h[j1]/(h[j1]+h[j])。 //求出u[j]和v[j]的初值 v[0]=1。=39。[39。j++) { h[j]=x[j+1]x[j]。i2。i++)ciny[i]。i++)cinx[i]。 cinn。 float s1[10]。//存放數(shù)組v[j] float d[11]。//存放數(shù)組x[j] float y[11]。數(shù) 值 分 析 上 機 題姓名：戴載星 學(xué)號：040139習(xí)題 438.（上機題）3次樣條插值函數(shù)（1）編制求第一型3次樣條插值函數(shù)的通用程序。 } return k。 k=(fabs(m[0]))。\t39。i9。\t39。i++){m[i]=x1[i]x0[i]。j++){sum=sum+a[i][j]*x0[j]。ji。}for(i=0。 w=t/50。t=99。 float a[9][9]={{31,13,0,0,0,10,0,0,0},{13,35,9,0,11,0,0,0,0},{0,9,31,10,0,0,0,0,0},{0,0,10,79,30,0,0,0,9},{0,0,0,30,57,7,0,5,0},{0,0,0,0,7,47,30,0,0},{0,0,0,0,0,30,41,0,0},{0,0,0,0,5,0,0,27,2}, {0,0,0,9,0,0,0,2,29}}。程序為：includeincludedefine eps //迭代誤差void main(void){ int i,j,l。x[i]endl。i++){cout39。j++)sum3+=a[i][j]*x[j]。i){x[n]=a[n][n+1]/a[n][n]。39。39。39。39。39。39。39。39。39。39。39。39。39。39。cout39。i++){for(j=1。q++)sum2+=a[i][q]*a[q][k]。i=n。qk。}for(j=k。j=n+1。i=n。s1[i]=a[i][k]sum。i++){float sum=0。//for(k=2。39。39。39。39。39。39。39。39。39。39。39。39。39。39。39。39。}//輸出進行第一步變換的數(shù)組acout39。j++)couta[i][j]39。 //第一列除以a[1][1]for(i=1。a[t][j]=b。}for(j=1。for(i=2。for(i=1。39。39。39。39。39。39。39。39。39。39。39。39。39。\t39。i=n。i++)for(j=1。cinn。習(xí)題335．（上機題）列主元三角分解法對于某電路的分析，歸結(jié)為求解線性方程組RI=V。對于不同得初始值收斂于不同的根, 在（∞，1）內(nèi)收斂于，在（，）內(nèi)收斂于，在（1，+∞）內(nèi)收斂于，但在內(nèi)（，1）和（－1，）均可能收斂于和。 flag=judgement(delta)。 int flag=1。 x0=x1。\t39。 type1=fabs(x1x0)。 float x0,x1,type,type1。 df=x*x1。 f=x*x*x/3x。 //輸出迭代的次數(shù)和根值}（2）計算迭代序列收斂于根的盡可能大的的函數(shù)為：includeincludevoid delay(int n) //定義延時函數(shù){for(n=10000。 } while(fabs(a)eps)。 do { a=f(x0)/df(x0)。}void main(void){ float x0,x1,a。}float df(float x) //定義函數(shù)df(x)，計算f(x)的導(dǎo)函數(shù){ float df。（3）通過本上機題，你明白了什么？解：（1）編制的通用程序：includeincludedefine eps /給定容許誤差float f(float x) //定義函數(shù)f(x){ float f。（2）給定方程，易知其有三個根。}從大到小的順序的值從小到大的順序的值精確值有效位數(shù)從大到小從小到大654436通過本上機題，看出按兩種不同的順序計算的結(jié)果是不相同的，按從大到小的順序計算的值與精確值有較大的誤差，而按從小到大的順序計算的值與精確值吻合。j=2。 sum+=s。（編制程序時用單精度）（4）通過本上機題，你明白了什么？第 6 頁 共 16 頁按從大到小的順序計算的通用程序為：includefloat sum(float N){ float j,s,sum=0。數(shù)值分析上機題姓名：陳作添 學(xué)號：040816習(xí)題120．（上機題）舍入誤差與有效數(shù)設(shè)，其精確值為。（3）按兩種順序分別計算，并指出有效位數(shù)。j++) { s=1/(j*j1)。for(j=N。}return sum。習(xí)題220．（上機題）Newton迭代法（1）給定初值及容許誤差，編制Newton法解方程根的通用程序。2．試取若干初始值，觀察當(dāng)，時Newton序列是否收斂以及收斂于哪一個根。 return(f)。 return (df)。 cinx0。 x0=x1。x0。}define eps float f(float x) //定義函數(shù)f(x){ float f。}float df(float x) //定義函數(shù)df(x)，計算f(x)的導(dǎo)函數(shù){ float df。}int judgement(float z){ int count=5。 type=fabs(x1)。type=type39。 if(fabs(x1x0)eps) return 1。 }void main(void){ float delta=0。 delta+=eps。 //輸出收斂的區(qū)間值}，程序計算出的δ的為δ=，即在區(qū)間（，）內(nèi)迭代序列收斂于0。產(chǎn)生上述結(jié)果的原因是區(qū)間不滿足大范圍收斂的條件。cout請輸入n的值:。i=n。 //給矩陣a賦值for(i=1。j++)couta[i][j]39。} //輸出數(shù)組acout39。39。39。39。39。39。39。39。39。39。39。39。