数值分析插值mtlab程序1.docVIP

%高斯赛德尔迭代解线性方程组a=input(请输入系数矩阵a:) %输入系数矩阵a disp(a)[m,n]=size(a)if m~=n %若a不是方阵则显示错误disp(false)endde=det(a);if de==0 %判断矩阵a是否为奇异矩阵disp(矩阵奇异)endb=input(请输入b矩阵:) %输入矩阵bdisp(b)x=input(请输入迭代初始值:)disp(x)k=x;jindu=10e-7 %迭代精度p=jindu+1;while p=jindu %迭代次数控制条件s=zeros(m,1);te=zeros(m,1);for i=1:m %迭代计算for j=1:i-1s(i)=s(i)+a(i,j)*k(j);endfor j=i+1:nte(i)=te(i)+a(i,j)*x(j);x=k;endk(i)=(b(i)-s(i)-te(i))/a(i,i)endif norm(k(i)-x(i),inf)jindu %迭代终止条件的控制p=abs(k(i)-x(i));elsebreak 拉格朗日算法 #includeiostream.h int main() { float x;//插值 float p[10][2];//已知(x0,y0),(x1,y1)... int n