很久没有编程了,请看一下这个程序的指针调用如何改正.

LittlesnowLinux

#include<iostream.h>

/* This function uses Gauss-elimination method
   to solve multi-variable  */
void Gauss(double *A,         /* A is nXn order matrix */
           double *b,         /* b is nX1 order matrix */
           double *x,         /* x is nX1 order matrix */
           int n              /* row and column of A */
           )
{
  int i,j,k;                  /* indices */
  double temp;

  for(i=0;i<n;i++)              /* convert A to top-triangle matrix */
    {
      for(k=i+1;k<n;k++)
        for(j=i;j<n;j++)
          {
            temp=A[k]/(A);
            A[k][j]=(A[k][j])-temp*(A[j]);
            b[k]=b[k]-temp*b;
          }
    }

  for(i=n-1;i>=0;i--)                /* get the solution of x */
    {
      if(i-n-1) x=b/A;
      else
        {
          for(j=i+1;j<n;j++)
            temp+=A[j]*x[j];
          x=(b-temp)/A;
        }
    }
  
}

int main()        /* main part of program */
{
  int n;        /* row and column of A or row of b */
  int i,j;        /* indexes */

  cout<<"请输入n:";
  cin>>n;

  double A[n][n],b[n],x[n],(*pa)[n],*pb,*px;        /* define matrixes */
  pa=A;
  pb=b;
  px=x;  

  cout<<"请输入A:\n";
  for(i=0;i<n;i++)
    for(j=0;j<n;j++)
      cin>>A[j];

  cout<<"请输入b:"<<endl;
  for(i=0;i<n;i++)
    cin>>b;

  Gauss(pa,pb,px,n);

  for(i=0;i<n;i++)
    cout<<x<<endl;
}

pupilzeng

(*pa)[n] -->*pa[n]

LittlesnowLinux

不行.关键在函数调用那里出错了.那个二维数组怎么调用啊.

pupilzeng

Sorry, 我刚刚在回答这个问题的时候，也没有动脑子
#define GET(A,i,j,n)     A[i*n + j]
定义的时候应该是
void Gauss(double *A, /* A is nXn order matrix */
double *b, /* b is nX1 order matrix */
double *x, /* x is nX1 order matrix */
int n /* row and column of A */
)
然后将这个函数定义内的A[j]全部改成GET(A,i,j,n)

调用的时候: Gauss(A, b, x, n)

或者使用gcc的扩充
这样定义函数
void Gauss(int n; /*forward declaration http://gcc.gnu.org/onlinedocs/gc ... tml#Variable-Length
*/
double A[n][n], /* A is nXn order matrix */
double *b, /* b is nX1 order matrix */
double *x, /* x is nX1 order matrix */
int n /* row and column of A */
)

crquan

所有楼上各位大哥：能不能用[ CODE ][ /CODE ]把代码格式化好看一点啊？会不会用啊？