Gauss

Gauss33 , Gauss44, Gauss55를 아우르는 일반적인 가우스 일반해

모든 크기의 행렬에 대해 사용 가능

function x=gauss(A,c)

% A: 계수행렬

% c: 상수행렬

% x: 연립방정식의 해

[m1,k1]=size(A);

% 상삼각행렬로 변환

for i=1:m1-1

m=-A(i+1:m1,i)/A(i,i);

A(i+1:m1,:)=A(i+1:m1,:)+m*A(i,:);

c(i+1:m1,1)=c(i+1:m1,1)+m*c(i,1);

end

x(m1,1)=c(m1,1)./A(m1,m1);

% 후진대입법을 사용하여 근 x를 구함

for i=m1-1:-1:1

x(i,1)=((c(i,1))-A(i,i+1:m1)*x(i+1:m1,1))./A(i,i);

end