数值分析 算法C语言程序(7)

数值分析 各个算法的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()

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说教育文库数值分析 算法C语言程序(7)在线全文阅读。