#include<stdio.h>
#include<math.h>
main()
{float a[10][10],b[10],s,t,e,sum;
int i,j,k,n,m;
printf("The top exp is ");
scanf("%d",&n);
for(i=0;i<n;i++)
for(j=0;j<n;j++)
scanf("%f",&a[i][j]);
for(i=0;i<n;i++)
scanf("%f",&b[i]);
scanf("%f",&e);
k=0;
do{t=a[k][k];
for(i=k;i<n;i++)
{if(fabs(t)<fabs(a[i][k]))
{t=a[i][k];
m=i;
}
else m=k;
}
if(fabs(t)<e)
printf("det A = 0\n");
else {if(m!=k)
{for(j=0;j<n;j++)
{s=a[m][j];
a[m][j]=a[k][j];
a[k][j]=s;
}
s=b[m];
b[m]=b[k];
b[k]=s;
}
for(i=k+1;i<n;i++)
for(j=k+1;j<n;j++)
{a[i][k]=a[i][k]/a[k][k];
a[i][j]=a[i][j]-a[i][k]*a[k][j];
b[i]=b[i]-a[i][k]*b[k];
}
}
k++;
}while(k<n-2);
if(fabs(a[n-1][n-1])<e)
printf("det A = 0\n");
else {b[n-1]=b[n-1]/a[n-1][n-1];
for(i=n-2;i>=0;i--)
{sum=0;
for(k=i+1;k<n;k++)
{sum+=a[k][j]*b[j];}
b[i]=(b[i]-sum)/a[i][i];
}
}
for(i=0;i<n;i++)
printf("%f\n",b[i]);
}
数值分析程序_c语言.rar_数值分析_数值分析程序_c语言_龙哥库塔法C
版权申诉
127 浏览量
2022-09-20
23:35:42
上传
评论
收藏 8KB RAR 举报
数值分析程序_c语言.rar （22个子文件） 
新建文件夹 
牛顿下山.C 483B 复合梯形法.c 630B 龙格库塔方法.C 450B 系统.c 318B 弦割法.c 389B 秦九韶.c 369B 高斯消去法.c 807B 杜氏分解法.C 907B 二分法.c 353B 牛顿插值多项式.c 492B 龙贝格算法.c 1017B 雅克比迭代.C 830B 三对角线追赶法.C 734B 综合测评.c 653B 埃特肯.c 333B 复合辛普森.c 291B 分段线性插值.c 752B I显示方法.c 206B 列主元元素消元.C 1KB 牛顿迭代.c 357B 改进欧拉法.C 344B 简单迭代.c 342B
评论0