2、A+temp);}3.n拉格朗日插值计算公式Ln(x)doubleLn(doublexx,intn,double*x,double*y);其中xx代表待求节点的横坐标值；n代表n次插值；x，y表示插值节点的横纵坐标值的存放地址。4.n次牛顿向前插值计算公式Nnbefore(x)doubleNnbefore(doublexx,intn,doublex0,doubleh,double*y);其中xx代表待求节点的横坐标值；n代表n次插值；y表示插值节点的纵坐标值的存放地址；x0表示初始插值节点坐标；h表示步长。5.n次牛顿向前插值计算公式Nnbehind(x)；double
3、Nnbehind(doublexx,intn,doublex0,doubleh,double*y);其中xx代表待求节点的横坐标值；n代表n次插值；y表示插值节点的纵坐标值的存放地址；x0表示初始插值节点坐标；h表示步长。6.杜力特尔三角分解求解线性方程组voidDLTR(intn,double*A,constdouble*b,double*x);n表示方程组的阶次；A表示存放系数矩阵的地址（如果是二维指针将首行地址传入即可）；b表示b向量；x是用来存放方程的解。通过该函数矩阵A将变成三角分解的紧凑格式，指针x中存放计算结果。7.乔立斯基三角分解求解对称正定方程组voi
4、dQLSJ(intn,double*A,constdouble*b,double*x);其中n表示阶次；A表示存放系数矩阵的地址（如果是二维指针将首行地址传入即可）；b表示b向量；x是用来存放方程的解。通过该函数矩阵变成了乔立斯基分解的对称三角阵，x指针存放解。8.输入数据函数，将数据输入到指针A所指向的空间，其中A为普通矩阵voidInputDataToMatrixA(intn,double*A);-15-吉林大学机械学院机械设计专业廉伟东学号20044121131.输入向量函数，将数据输入到指针b所指向的空间voidInputDataVector(intn,doub
5、le*b);2.输出矩阵A中的数据按照n行n列形式输出voidOutputDataMatrix(intn,constdouble*A);3.输出向量函数voidOutputDataVector(intn,constdouble*x);4.输入对称矩阵函数，将数据输入到指针A所指向的空间voidInputDataToSymMatrixA(intn,double*A);只需要输入上三角矩阵即可，按照从左向右从上向下的顺序输入数据5.输出对称矩阵A中的数据按照n维n列形式输出voidOutputDataSymMatrix(intn,double*A);6.高斯赛德尔迭代法vo
6、idGSSDE(double*a,double*b,intn,double*x,doubleeps);其中a表示系数矩阵；b表示向量b；其中x用来存放方程的解；eps设置误差的大小，并且该函数按照来作为循环的出口。经过该函数并不该表矩阵a，和向量b的值7.雅可比迭代法voidYKB(double*a,double*b,intn,double*x,doubleeps);其中a表示系数矩阵；b表示向量b；其中x用来存放方程的解；eps设置误差的大小，并且该函数按照来作为循环的出口。经过该函数并不该表矩阵a，和向量b的值程序的源代码//数值方法.cpp:Definesthee
7、ntrypointfortheconsoleapplication.#include"stdafx.h"#include"iostream.h"#include"math.h"inlinedouble&Data(double*A,intn,intx,inty){//获得矩阵的x行j列，完全的矩阵return*(A+x*n+y);}inlinedouble&DataSym(double*A,intn,intx,inty){//获得矩阵的x行j列，对称的矩阵if(x>y){intt;t=x;x=y;y=t;}inttemp=0;for(int
