数值分析 各个算法的C语言程序
#include<stdlib.h> #include<math.h> void main() {
double A[3][3]={2,-1,1,-1,-2,3,1,3,1},x[3],y[3],b[3]={4,5,6},d[3],L[3][3],U[3][3],f1=0,f2=0; int i,j,k,n=3; for(i=0;i<n;i++)
for(j=0;j<n;j++) {U[i][j]=0;L[i][j]=0;} d[0]=A[0][0];L[0][0]=1; for(i=1;i<n;i++) {
L[i][i]=1;
for(j=0;j<=i-1;j++) {for(k=0;k<=j-1;k++) f1+=U[i][k]*L[j][k];
U[i][j]=A[i][j]-f1;L[i][j]=U[i][j]/d[j];f2+=U[i][j]*L[i][j];f1=0;} d[i]=A[i][j]-f2;f2=0;} y[0]=b[0];
for(i=1;i<n;i++) {for(j=0;j<=i-1;j++) f1+=L[i][j]*y[j];y[i]=b[i]-f1;f1=0;} x[n-1]=y[n-1]/d[n-1];
for(i=n-2;i>=0;i--) {for(j=i+1;j<n;j++) f2+=L[j][i]*x[j];x[i]=y[i]/d[i]-f2;f2=0;} printf("输出L矩阵：\n"); for(i=0;i<n;i++)
{for(j=0;j<n;j++) printf("%f ",L[i][j]); printf("\n");} printf("输出U矩阵：\n"); for(i=0;i<n;i++)
{for(j=0;j<n;j++) printf("%f ",U[i][j]); printf("\n");} printf("输出求解结果：\n");
for(i=0;i<n;i++) printf("%f ",x[i]); printf("\n");
}
九、 追赶法
#include<stdio.h> #include<stdlib.h> #include<math.h> void main()