#include<iostream>
using namespace std;
void main()
{
int i=0,j=0,k,n,sum=0;
cout<<"输入未知数个数:";
cin>>n;
double **a=new double*[n+1];
for(i=0;i<=n;i++)
{
a[i]=new double[n+1];
}
double *b=new double[n+1];
cout<<"输入线性方程的系数矩的每一行和等号右侧结果:"<<endl;
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
{
cin>>a[i][j];
}
cin>>b[i];
}
k=1;
while(1)
{
for(i=k+1;i<=n;i++)
{
a[i][k]=a[i][k]/a[k][k];
}
for(i=k+1;i<=n;i++)
{
for(j=k+1;j<=n;j++)
{
a[i][j]=a[i][j]-a[i][k]*a[k][j];
}
b[i]=b[i]-a[i][k]*b[k];
}
if(k!=(n-1)) k++;
else break;
}
b[n]=b[n]/a[n][n];
for(i=n-1;i>=1;i--)
{
for(j=i+1;j<=n;j++)
{
b[i]-=a[i][j]*b[j];
}
b[i]=(1/a[i][i])*b[i];
sum=0;
}
for(i=1;i<=n;i++)
{
cout<<"x"<<i<<" = "<<b[i]<<endl;
}
for(i=0;i<=n;i++)
{
delete []a[i];
}
delete []a;
delete []b;
}
n维线性方程组,哪里错了求大神指点