function r=maxrad(x,g0)
g=x(1)+i*x(2);
r=max(abs(g-g0));

