#include <stdio.h>
#include <math.h>
#include <stdlib.h>

int rgauss(n,a,b)
int n;
double a[],b[];
{
int *js,L,i,j,k,is,p,q;
double d,t;
js=malloc(n*sizeof(int));
L=1;
for (k=0;k<=n-2;k++)
{ d=0.0;
for (i=k;i<=n-1;i++)
for (j=k;j<=n-1;j++)
{t=fabs(a[i*n+j]);
if (t>d){d=t; js[k]=j; is=i; }
}
if (d+1.0== 1.0) L=0;
else
{ if (js[k]!=k)
for (i=0;i<=n-1;i++)
{p=i*n+k; q=i*n+js[k];
t=a[p];a[p]=a[q];a[q]=t;
}
if (is!=k)
{ for(j=k;j<=n-1;j++)
{ p=k*n+j; q=is*n+j;
t=a[p]; a[p]=a[q]; a[q]=t;
}
t=b[k]; b[k]=b[is]; b[is]=t;
}
}
}
if (L==0)
{free(js); printf("fail\n");
return (0);
}
d=a[k*n+k];
for (j=k+1;j<=n-1;j++)
{ p=k*n+j; a[p]=a[p]/d; }
b[k]=b[k]/d;
for (i=k+1; i<=n-1;i++)
{for (j=k+1;j<=n-1;j++)
{p=i*n+j;
a[p]=a[p]-a[i*n+k]*a[k*n+j];
}
b[i]=b[i]-a[i*n+k]*b[k];
}
}
d=a[(n-1)*n+n-1];
if (fabs(d)+1.0==1.0)
{ free(js); printf("fail\n");
return (0);
}
b[n-1]=b[n-1]/d;
for (i=n-2;i>=0;i--)
{t=0.0;
for(j=i+1; j<=n-1;j++)
t=t+a[i*n+j]*b[j];
b[i]=b[i]-t;
}
js[n-1]=n-1;
for(k=n-1; k>=0;k--)
if(js[k]!=k)
{ t=b[k]; b[k]=b[js[k]]; b[js[k]]=t; }
free(js);
return (1);


上面是我在做数值方法是看到的一个高斯消元的消元过程程序,在开头
(
int rgauss(n,a,b)
int n;
double a[],b[];
)
里,int n;为什么要放在int rgauss(n,a,b)和double a[],b[];之间啊?

这是以前Ｃ语言函数的写法(旧标准)
返回值类型　function(参数名字列表)
参数类型　变量１；
参数类型　变量２；
……

上面那位高手,请问,可以把int n; double a[],b[];放在int rgauss(n,a,b)前面吗?

这样写就可以了：
int rgauss(int n,double a[],double[])

[此贴子已经被作者于2006-11-5 19:51:02编辑过]

int rgauss(n,a,b)
int n;
double a[],b[];
是一种很古老的写法

int n;
double a[],b[];
是不是全局变量。

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
int c,d;
int rgauss(n,a,b)
int n;
double a[],b[];
{
//函数体；
}
c,d是全局变量 而 n ,a[],b[] 是rgauss()函数的局部变量

这种写法少见.

把
int n;
double a[],b[];
放在里面如:
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
int c,d;
int rgauss(n,a,b)
{
int n;
double a[],b[];
……
}
可以吗？