cramer.m

来自「programas matlab metodos numericos」· M 代码 · 共 22 行

% resolucion por Cramer
n=length(A)-1;
d=det(A(1:n,1:n));
if d==0
    %no es compatible determinado
    A(1:n,1)=A(1:n,n+1);
    if det(A(1:n,1:n))==0
        'Sistema Compatible Indeterminado'
    else
        'Sistema Incompatible'
    end
else
    % compatible determinado
    for i=1:n %Calcula los det(Xi) y calcula Xi
        B(1:n,1:n)=A(1:n,1:n);
        B(1:n,i)=A(1:n,n+1);
        x(i)=det(B(1:n,1:n))/d;   
    end
    'Sistema Compatible Determinado'
    x % muestra solucion
end