超定方程组的最小二乘解

#include<iostream.h>
#include<iomanip.h>
#include<math.h>
#define m 4//n<m
#define n 2
void line_fun(double a[n][n],double b[n],double p[n])//列主高斯消元求解线性方程组
{
double x[n];//线性方程组的解  
double mm[n];//辅助乘数
double T[n],BT=0;//存储换行临时变量
double temp=0;
double t;//最大列主元
int tap1=0,tap2=0;//最大列主元下标
int i,j,s;
for(s=0;s<n;s++)
{
t=a[s][s];//赋初值
for(int p=s;p<n;p++)//选取最大列主元
{
if(fabs(a[p][s])>t)
{
t=a[p][s];
tap1=p;
tap2=s;
}
}
if(t==0)
{
cout<<"失败"<<endl;
break;