高斯迭代法解矩阵方程

高斯迭代法解矩阵方程function[a,b]=gaussj(a,n,np,b,m,mp)forj=1:nipiv(j)=0;endfori=1:nbig=0;forj=1:nif(ipiv(j)~=1)fork=1:nif(ipiv(k)==0)if(abs(a(j,k))>=big)big=abs(a(j,k));irow=j;icol=k;endelseif(ipiv(k)>1)pause,dips