
function X=cramer(A,b)
%-----------------------------------------------------------------------
%             -------------                             -------
%             function file                             cramer.m
%             -------------                             -------
%-----------------------------------------------------------------------
%                        ----------------------------
%                        Zeros Reais de Sistemas Lineares
%                        ----------------------------
%                     ---------------------------------------
%                             Regra de Cramer
%                     ---------------------------------------
%
%            Chamada do Procedimento  Cramer
%------------------------------------------
%X =cramer(A,B)
%
% Entrada: A (N x N) matriz dos coeficientes (não singular);
%          B (N x 1) matriz dos termos independentes;
% Saída: X (N x 1) solução de AX=B.
%
%Veja Tambem: gauss gausspiv jacobi sor
%--------------------------------------------------------
[N N] = size(A);
n1=size(A);
n2=size(b);

if n1(1)~=n1(2)
    beep
    error('a');
end
if n2(2)>n2(1)
    b=b';
    n2=size(b);
end
if n1(1)~=n2(1)
    beep
    error('inconsistencia');
end

Delta=det(A);
DeltaX=[];

for i=1:n1(1)
    Ax=A;
    Ax(:,i)=b;
    DeltaX=[DeltaX,det(Ax)];
end


%Verificacao do sistema:


X=DeltaX/Delta;


    
