function x=gausselim(A,b)
n=size(A,1);
for i=1:n-1
[foo,p]=max(abs(A(i:n,i)));
p=p+i-1;
A([i,p],:)=A([p,i],:);
b([i,p])=b([p,i]);
for j=i+1:n
m=A(j,i)/A(i,i);
for k=i:n
A(j,k)=A(j,k)-m*A(i,k);
end
b(j)=b(j)-m*b(i);
end
end
x=b;
for i=n:-1:1
for j=i+1:n
x(i)=x(i)-A(i,j)*x(j);
end
x(i)=x(i)/A(i,i);
end
慕雪6442864
相关分类