matlab中用高斯消元法接方程组程序

function [X]=gaussian(A,B,n)C=[A B];for i=2:n [r,c]=max(abs(C(i-1:n,i-1))); VEC=C(i-1,:);C(i-1,:)=C(c+i-2,:);C(c+i-2,:)=VEC;if C(i-1,i-1)==0 returnendfor k=(i-1):(n-1)C(k+1,:)=C(i-1,:)*-C(k+1,i-1)/C(i-1,i-1)+C(k+1,:);endendD=C;B2=C(:,n+1);D(:,n+1)=[];X=zeros(1,n);for i=1:nX(n+1-i)=(B2(n+1-i)-D(n+1-i,:)*X')/C(n+1-i,n+1-i);end运行:>>A=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];>>B=[6,25,-11,15]';>>gaussian(A,B,4)ans =1。0000 2。0000 -1。0000 1。0000