C语言，C++，数值分析实验，Jacobi迭代，高斯赛德尔迭代

#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;
 }

实验结果