用平方根法解线性方程（c++)

#include <iostream>
#include <iomanip>
#include <cmath>
int k;
const int n=4;
using namespace std;
double q_s1(double array[][n],int i,int j);
double q_s2(double array[][n],int j);
double q_s3(double array[][n],double Y[],int i);
double q_s4(double array[][n],double X[],int i);
int main()
{
    int i,j;
 double a[n][n];
 cout<<"please input a[i][j]:"<<endl;
 for(i=0;i<n;i++)                         //输入a[i][j]
 {
  for(j=0;j<n;j++)
   cin>>a[i][j];
 }
 double b[n];
 cout<<"please input b[i]:"<<'/n';
 for(i=0;i<n;i++)                            //输入b[i]
  cin>>b[i];
 double l[n][n];
 for(j=0;j<n;j++)                    //求l[i]
 {
  for(i=j;i<n;i++