ni hao.cpp

用C++实现的高斯塞德尔的数值分析实验作业

#include <iostream>
#include <math.h>
using namespace std;
float *GauseSeidel(float *a,int n)
{
 
 int   i,j,k=0;
 float *x,y,e;
 cout<<"input e:"<<endl;
 cin>>e;
 x =new float[n];
 for( i=0;i<=n-1;i++ )
  x[i]=0;
 do
 { 
  for( i=0;i<=n-1;i++ )
  {
  y=0;
   for( j=0;j<=n-1;j++ )     
    y+=*(a+i*(n+1)+j)*x[j];   
   y=(*(a+i*(n+1)+n)-y)/(*(a+i*(n+1)+i));
   x[i]+=y;
  }  
k++;
  cout<<k<<"\t"; 
for( i=0;i<3;i++ )
  cout<<"x["<<i<<"]="<<x[i]<<"\t";
cout<<endl;
 }
 while(fabs(y)>e);
 return x;
}
void main()
{
 int i;
 float *x;
 float c[12]={8,-3,2,20,4,11,-1,33,6,3,12,36}; 
 float *GauseSeidel(float *,int);
 x=GauseSeidel(c,3);
}