function theta  = compute_theta( t );

theta = sign(t(:,2)).*acos( t(:,1));

