function [C]=cramer(m)
  ms = size(m);
  ml = ms(1,1);
  mc = ms(1,2);
  for a = 1:ml
    for b = 1:ml
      for c = 1:ml
        A(b,c)=m(b,c);
        B(b,c)=m(b,c);
      endfor
    endfor
    for b = 1:ml
      B(b,a)=m(b,mc);
    endfor
    C(a,1)=det(B)/det(A);
  endfor
endfunction