
function [g1,g2] = makez(g1,g2,k);
%
% make g2(k) = 0, by rotation on (g1,g2);
%

t = atan(g2(k)/g1(k));
v1 = cos(t)*g1+sin(t)*g2;
v2 = sin(t)*g1-cos(t)*g2;
g1 = v1;
g2 = v2;

if max(g1)+min(g1) < 0, g1 = -g1; end
if max(g2)+min(g2) < 0, g2 = -g2; end

