用c++写一个万能的高斯-赛德尔迭代法怎样写啊 哪位高手能告诉我啊

来源:百度知道 编辑:UC知道 时间:2024/09/23 01:33:25
之所为万能,就是无论是几维数组它都能数出来

这是 C 程序,你把头文件换成 c++,把 printf 改 cout 就可以了。

无论是几维数组 ?-- 数组大小还是维数?

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define N 100
float *GauseSeidel(float *a,int n)
{
int i,j,nu=0;
float *x,dx;
x=(float *)malloc(n*sizeof(float));
for(i=0;i<=n-1;i++)
x[i]=0.0;
do {
for(i=0;i<=n-1;i++) {
float d=0.0;
for(j=0;j<=n-1;j++)
d+=*(a+i*(n+1)+j)*x[j];
dx=(*(a+i*(n+1)+n)-d)/(*(a+i*(n+1)+i));
x[i]+=dx;
}
if(nu>=N)
{
printf("iter divergence\n");
exit(1);
}
nu++;
}
while(fabs(dx)>1e-6);
return x;
}
void main()
{
int i;
float *x;
float c[12]={5,2,1,8,2,8,-3,21,1,-3,-6,1};
float *GauseSeidel(