2.全部源程序
#include <iostream.h>
#include <iomanip.h>
#include <math.h>
#include <stdlib.h>
#include <fstream.h>
#define e 0.000001; //定义精度水平
double vect_x[11]; //题目中已知的两组xy值
double vect_y[21];
double vect_x2[8];
double vect_y2[5];
double tempG[2000]; //系数矩阵
double res[11*21]; //拟合结果p(x,y)
int res_num;
double insert_valuez[11*21]; //通过方程及插值得到z(x,y)
struct one_root //定义解结构体
{
double t,u,v,w;
};
one_root root[11*21];
//一些矩阵及向量的初等计算
//矩阵与矩阵乘法
void MM_mul(double a[],double b[],double c[],int dim)
{
int i,j,z;
double k=0;
for(i=0;i<dim;i++)
{
for(j=0;j<dim;j++)
{
for(z=0;z<dim;z++)
{
k=k+a[i*dim+z]*b[z*dim+j];
}
c[i*dim+j]=k;
k=0;
}
}
}
//矩阵与向量乘法
void MV_mul(double a[],double b[],double c[],int dim)
{
int i,j;
double k=0;
for(i=0;i<dim;i++)
{
