function [pos_ou] = f_ou(xg,yg,xa,ya)

ae=6378.137*10^3;
rg=ae+20000*10^3;

teta1=atan(yg/xg);
teta2=atan(ya/xa);

if ( xg<0 )
   teta1=teta1+pi;
end

if ( xa<0 )
   teta2=teta2+pi;
end

if (teta1<teta2)
   pos_ou=1;
else
   pos_ou=0;   
end