北航数值分析第一次大作业（高斯gauss lu分解）(2)

#define N 20 void main() { long double a[N][N],l[N][N],u[N][N]; long double b[N],x[N],y[N],s[N],temp; long double max; int i,j,k,t; int M[N];

//*************************原始方程的计算和输出************************ //方程组原始条件的代入 for(i=0;i

printf(\printf(\选主元的Doolittle消去法*****************\\n\printf(\//输出原方程组系数矩阵

printf(\原方程组系数矩阵A：\\n\for(i=0;i

printf(\第%d列至第%d列：\\n\ printf(\ if((j+1)%5==0) printf(\ } }

//输出原方程组右端列向量

printf(\原方程组右端列向量b：\\n\for(i=0;i

//*****************************分解过程******************************** for(k=0;k

{ s[i]=a[i][k]; for(t=0;t<=k-1;t++) s[i]-=l[i][t]*u[t][k]; } max=0; for(i=k;i

//输出LU矩阵

printf(\矩阵：\\n\for(i=0;i

}

}

}

printf(\第%d列至第%d列：\\n\if(j=i) printf(\if((j+1)%5==0) printf(\

//求Qb

for(k=0;k

//求解y和x y[0]=b[0];

for(i=1;i

x[N-1]=y[N-1]/u[N-1][N-1]; for(i=N-2;i>=0;i--) { x[i]=y[i]; for(t=i+1;t

//输出结果

printf(\方程组的解**************************\\n\for(i=0;i