用LU分解法求解n阶线性代数方程AX=b c++实现

#include<iomanip.h>
#include<iostream.h>
#include<stdlib.h>
void main()
{
double a[10][10];
double b[10];
int n;
int i;
int k;
int j;
cout<<"请输入矩阵a的阶数 "<<endl;
cin>>n;
int max=n;

    W:cout<<"请以行的顺序输入矩阵a "<<endl;
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
cin>>a[i][j];
cout<<endl;
}
cout<<"请输入矩阵b"<<endl;
for( i=0;i<n;i++)
cin>>b[i];

for(i=1;i<n;i++)
a[i][1]=a[i][1]/a[0][0];
for(k=1;k<n;k++)
{