Question

哪位大哥帮我看看这个Matlab程序，是关于高斯消去法球逆的

我是在原矩阵的基础上改变，结果输出时候出现其他不知怎样的矩阵，最后才出现逆矩阵，下面是代码：
function A=gaosi_inv(A)
n=size(A);
if n(1)~=n(2)
error('不是方阵');
end
n=n(1);
flag=1:n;
for k=1:n
t=find(abs(A(k:n,k))==max(abs(A(k:n,k))));
t=t(1)+k-1;
flag(k)=t;
if t~=k
p=A(k,:);A(k,:)=A(t,:);A(t,:)=p;
end
if A(k,k)==0
error('矩阵不可逆');
end
h=1/A(k,k);A(k,k)=h;
A([1:k-1 k+1:n],k)=A([1:k-1 k+1:n],k)*(-h)
for i=1:n
if i~=k
A(i,[1:k-1 k+1:n])=A(i,[1:k-1 k+1:n])+A(k,[1:k-1 k+1:n])*A(i,k);
end
end
A(k,[1:k-1 k+1:n])=A(k,[1:k-1 k+1:n])*h;
end
for k=n-1:-1:1
t=flag(k);
if k~=t
p=A(:,t);A(:,t)=A(:,k);A(:,k)=p;
end
end

Answer 1

函数定义有问题，输入参数和输出参数相同。编译不会出现问题。
建议function Ainv=gaosi_inv(A)
然后程序的最末尾加上一句Ainv=A;
OK

Answer 2

A([1:k-1 k+1:n],k)=A([1:k-1 k+1:n],k)*(-h) 这行后面没分号

这种事matlab的editor都会提示的