#include<stdio.h>
int main()
{ int n,i,j,k,b;
double num1;
double num2;
printf("输入方程组的维数:\n");
scanf("%d",&n);//方程组维数
double arrD[n][n]={0};//系数矩阵D
double arrB[n][n]={0};//矩阵D-1*(l+U)
printf("输入方程组各个系数:\n");
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{ scanf("%lf",&num1);
arrD[i][j]=num1;
}
}
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{ if(i==j)
{
arrB[i][j]=0;
}
else{
arrB[i][j]=-arrD[i][j]/arrD[i][i];
}
}
}
printf("输入矩阵b:\n");
double arrf[n]={0};//矩阵f
for(b=0;b<n;b++)
{
scanf("%lf",&num2);
arrf[b]=num2/arrD[b][b];
}
double sum=0;
double arrx[n]={0};
double arrxd[n]={0};//初始X向量
int ddcs=0;//迭代次数
printf("请输入迭代次数:\n");
scanf("%d",&ddcs);
int dd;
int m;
for(dd=0;dd<ddcs;dd++)
{
for(k=0;k<n;k++)
{
for(j=0;j<n;j++)
{
sum=sum+arrB[k][j]*arrx[j];
}
arrxd[k]=sum+arrf[k];
sum=0;
}
for(m=0;m<n;m++)
{
arrx[m]=arrxd[m];
}
printf("Jacobi迭代%d次后的结果为:",dd+1);
for(k=0;k<n;k++)
{
printf("%lf ",arrxd[k]);
if(k==(n-1))
{
printf("\n");
}
}
}
return 0;
}
实验结果
#include<stdio.h>
int main()
{ int n,i,j,k,b;
double num1;
double num2;
printf("输入方程组的维数:\n");
scanf("%d",&n);//方程组维数
double arrD[n][n]={0};//系数矩阵D
double arrB[n][n]={0};//矩阵D-1*(l+U)
printf("输入方程组各个系数:\n");
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{ scanf("%lf",&num1);
arrD[i][j]=num1;
}
}
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{ if(i==j)
{
arrB[i][j]=0;
}
else{
arrB[i][j]=-arrD[i][j]/arrD[i][i];
}
}
}
printf("输入矩阵b:\n");
double arrf[n]={0};//矩阵f
for(b=0;b<n;b++)
{
scanf("%lf",&num2);
arrf[b]=num2/arrD[b][b];
}
double sum=0;
double arrx[n]={0};
double arrxd[n]={0};//初始X向量
int ddcs=0;//迭代次数
printf("请输入迭代次数:\n");
scanf("%d",&ddcs);
int dd;
int m;
for(dd=0;dd<ddcs;dd++)
{
for(k=0;k<n;k++)
{
for(j=0;j<n;j++)
{
sum=sum+arrB[k][j]*arrxd[j];
}
arrxd[k]=sum+arrf[k];
sum=0;
}
printf("高斯赛德尔迭代迭代%d次后的结果为:",dd+1);
for(k=0;k<n;k++)
{
printf("%lf ",arrxd[k]);
if(k==(n-1))
{
printf("\n");
}
}
}
return 0;
}
实验结果