function [flag,x,y,angle_rad,phi] = g_ini(phi_l_rad,r_leo,r1,xl,yl)


phi=asin((r_leo/r1)*sin(phi_l_rad));
flag = 1;

if (isreal(phi) == 1)
   
alpha=phi-phi_l_rad;

teta=atan(yl/xl);

	if ( xl < 0 )
 	  teta=teta+pi;
	else
	   teta=teta;
   end
   
   angle_rad = teta+alpha;
   x=r1*cos(angle_rad);
   y=r1*sin(angle_rad);

   
else
   
   flag = 0;
   phi=0;
   angle_rad=0;
   x=xl;
   y=yl;
   
end