mgauss.m

来自「various matlab programs to slove various」· M 代码 · 共 7 行

function [Mk,Mkinv]=MGauss(A,k,n)
Mk=eye(n);  
for i=k+1:n  
   Mk(i,k)=-A(i,k)/A(k,k);
end 
Mkinv=2*eye(n)-Mk; 
return