高斯赛德尔迭代法matlab编程(共2页).doc

精选优质文档-倾情为你奉上functionx,k=GaussSeidel(A,b,x0,eps,M)%高斯赛德尔迭代法求方程组的解（矩阵公式求解）%A为方程组的系数矩阵；b为方程组的右端项%x为线性方程组的解了；x0为迭代初值%eps为误差限；M为迭代的最大次数if nargin=3 eps= 1.0e-6;%默认精度 M = 10000;%参数不足时默认后两个条件elseif nargin =4 M = 10000;%参数的默认值elseif nargin3 error(参数不足); returnendn,m=size(A);nb=length(b);%当方程组行与列的维数不相等时，停止计算，并输出出错信息if n=m error(矩阵A行数和列数必须相等!); return;end%当方程组与右端项的维数不匹配时，停止计算，并输出出错信息if n=nb error(矩阵A的行数必须和b的长度相等!);