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