數(shù)值分析上機(jī)報(bào)告-wenkub.com
2025-08-01 00:46 本頁(yè)面
【正文】 能干的人，不在情緒上計(jì)較，只在做事上認(rèn)真；無(wú)能的人！不在做事上認(rèn)真，只在情緒上計(jì)較。寧可累死在路上，也不能閑死在家里！寧可去碰壁，也不能面壁。printf(\b\n)。jN+1。}n++。if(n==50||n==90||n==170||n==200)/*得到所需的數(shù)據(jù)*/{for(i=0。k[2][3]=h*(10001000*(y[1]+k[2][0]k[2][1]+k[2][2])100*(y[2]+k[2][0]k[2][1]+k[2][2]))。k[1][2]=h*(y[2]+(1/)*k[1][0]+k[1][1])。k[0][1]=h。i++)/*計(jì)算每一個(gè)k值*/k[l][i]=0。do{for(l=0。一、分析如下：二、 程序如下：includeincludedefine N 3float y[N]={0,0,0},k[N][N+1]。 y=pow(3,x)*pow(x,)*(5*x+7)*sin(x*x)。 } if(abs(t[0][k+1]t[0][k])eps) break。k=8。 } t[k][1]=t[k1][1]+h*tt/pow(2,k1)。i=pow(2,k1)。 t[0][1]=(ya+yb)*h/2。 tt=0。 a=1。（2）計(jì)算新分點(diǎn)的值時(shí)，這些數(shù)值的個(gè)數(shù)成倍增加。應(yīng)用Romberg算法時(shí)，可以先給出所要求的誤差ε，當(dāng)就停機(jī)。第五題：用Romberg算法求（允許誤差ε= 01）。printf(“%f”,x2)。x2=x1f/f’。一、 分析如下：對(duì)滿足下列條件的函數(shù)：i. f(a)f(b)0；ii. f”(x)在區(qū)間[a,b]上不變號(hào)；iii. f’(x)≠0；iv. |f(c)|/(ba)≤|f’(c)|，其中c 是a,b中使min(|f’(a)|,|f’(b)|)達(dá)到的一個(gè)。 else printf(f=%,f1=%,f[9],f1[9])。 else if(7=aamp。 else if(5=aamp。 else if(3=aamp。 else if(1=aamp。 \*f(x)表達(dá)式*\ f1[i]=x[i]*pow((axm[i1]),2)/2+c1[i]x[i1]*pow((xm[i]a),2)/2。i10。 \*用追趕法求解三樣條插值法中的M*\printf(please input:)。}x[9]=c[9]/b[9]。i10。d[9]=6*((y[9]y[8])/h)/h。 \*輸入矩陣值*\d[0]=6*((y[1]y[0])/h1)/h。i9。i10。double b[10],l[10]。double y[10]={0,},xm[10]={1,2,3,4,5,6,7,8,9,10}。求的方法為：這里最終歸結(jié)為求解一個(gè)三對(duì)角陣的解。 /* 主函數(shù)返回z 告訴系統(tǒng) 程序結(jié)束 */ }三、運(yùn)行結(jié)果：四．問(wèn)題討論：題中矩陣B是對(duì)稱正定陣，且是三對(duì)角的，所以選擇合適的松馳因子ω，收斂速度是很快的。 printf( x%d=%,i,x[i])。i) x[i]=(y[i]t[i][i+1]*x[i+1])/u[i]。 y[i]=d[i]q[i]*y[i1]。 for(i=1。 printf(x%d=%,i,x1[i])。 x1[i]=x1[i]*(1w)+w*s。 for(j=0。h9。} for(i=0。i++)
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1