北航大學(xué)數(shù)值分析計(jì)算實(shí)習(xí)報(bào)告(文件)
【正文】 否判斷B_K1是否為0？ wucha=(fabs(B_k1B_k0))/(fabs(B_k1))。 } else B_k0 = B_k1。 double b[501] = {0}。 get_y(y,u)。 B_k1=NeiJi(y,u)。 else if(i10000) { printf(迭代次數(shù)超長(zhǎng)，請(qǐng)更改初始向量\n)。 return (1/B_k1)。i = n。t++) { e+=C[it+s+11][t1]*y[t1]。i=1。t++) e += C[it+s+11][t1]*u[t1]。i++) { b[i] = y[i]。k = n。 for(t = max(max(1,kr),js)。 } if(k == n) break。 for(t = max(max(1,ir),ks)。 } } }int min(int a,int b){ if(ab) return a。} 附錄：部分實(shí)驗(yàn)程序截圖 另，其余都為0 ，其余為0 ，其余為01，其余為0 25。}int max(int a,int b){ if(a=b) return a。t++) { e = e + C[it+s+11][t1]*C[tk+s+11][k1]。i = min(k+r,n)。t++) { e = e + C[kt+s+11][t1]*C[tj+s+11][j1]。j = min(k+s,n)。 double e。 }}void copy(double b[501],double y[501]){ for(int i=0。 for(t = i+1。 } u[n1] = y[n1]/C[s+11][n1]。 for(t = max(1,ir)。 double e。 } else B_k0 = B_k1。//get wucha // wucha=(fabs(B_k1B_k0))/(fabs(B_k1))。//保護(hù)y向量 Doolittle_getx(C,y,u,501,2,2)。 int i=0。 return B_k1。 else if(i10000) { printf(迭代次數(shù)超長(zhǎng)，請(qǐng)更改初始向量\n)。 get_u(u,y,a)。 double wucha。 } for(i = 1。 } for(i = 1。 for(i = 2。i++) u[i]=c*y[i2]+b*y[i1]+a[i]*y[i]+b*y[i+1]+c*y[i+2]。 }}void get_u(double *u,double *y,double *a){ u[0]=a[0]*y[0]+b*y[1]+c*y[2]。i501。 b++。 for(int i=0。i=501。 else *p=1。 p++。 return e。 double e =1。 B[i1] +=Uk[i1]。j501。i++) {
試題試卷相關(guān)推薦
鄂ICP備17016276號(hào)-1