% *************
% * crame law *
% *************
% input : augmented matrix C
% output: row vector X

function [x] = cramelaw (c)

N=size(c);
N(2)=[];
x=zeros(1,N);
jacobimatrix=det(c(:,1:N));
for j=1:N
    x(j)=det([c(:,1:j-1),c(:,N+1),c(:,j+1:N)])/jacobimatrix;
end