function cgeig_test()
[A,b]=matrix5p(200);
toler =1e-6;
itmax = 500;
x0=zeros(size(b));
x=cgeig(A,b,x0,toler,itmax);