在c++下编写高斯列主元消去解方程组

//本程序用文件输入输出
#include<iostream.h>
#include<fstream.h>
#include<math.h>
#include<stdlib.h>

void init(double[100][100],int &);
void work(double[100][100],int);
void pr(double[100][100],int);

void main()
{
double x[100][100];
int n;
init(x,n);
work(x,n);
pr(x,n);
}

void init(double x[100][100],int &n)
{
ifstream infile;
infile.open("equations.in",ios::in);
infile>>n;
for(int i=0;i<n;i++)
for(int j=0;j<n+1;j++)
infile>>x[i][j];
infile.close();
}

void work(double x[100][100],int n)
{
int q[100],p,i,j,k;
double min,l;
for(i=0;i<n;i++)
q[i]=0;
for(k=0;k<n;k++)
{
p=-1;
min=10000;
for (i=0;i<n;i++)
if (!q[i] && x[i][k]!=0 && min-1e-10>fabs(x[i][k]))
{
min=fabs(x[i][k]);