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

数值分析 各个算法的C语言程序

f2=0; int i,j,k;

for(i=0;i<3;i++) for(j=0;j<3;j++)

{U[i][j]=0;L[i][j]=0;} for(i=0;i<3;i++)

{U[0][i]=A[0][i];L[i][0]=A[i][0]/U[0][0];L[i][i]=1;} for(i=1;i<3;i++) for(j=i;j<3;j++)

{for(k=0;k<=i-1;k++){f1=f1+L[i][k]*U[k][j];f2+=L[j][k]*U[k][i];} U[i][j]=A[i][j]-f1;L[j][i]=(A[j][i]-f2)/U[i][i];f1=0;f2=0;} y[0]=b[0];

for(i=1;i<3;i++)

{for(j=0;j<=i-1;j++) f1+=L[i][j]*y[j];y[i]=b[i]-f1;f1=0;} x[2]=y[2]/U[2][2]; for(i=1;i>=0;i--)

{for(j=i+1;j<3;j++) f2+=U[i][j]*x[j];x[i]=(y[i]-f2)/U[i][i];f2=0;} printf("输出L矩阵：\n"); for(i=0;i<3;i++)

{for(j=0;j<3;j++)

printf("%f ",L[i][j]); printf("\n");}

printf("输出U矩阵：\n"); for(i=0;i<3;i++)

{for(j=0;j<3;j++)

printf("%f ",U[i][j]); printf("\n");}

printf("输出求解结果：\n"); for(i=0;i<3;i++)

printf("%f ",x[i]); printf("\n");

}

八、 改进的平方法

#include<stdio.h>

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