function x = CG(A,b)
  x = zeros(size(A,1),1);
  x = CGx(A,b,x);