circle.m

function[Vx,Vy,x_hc,y_hc,phi_hc]=circle(i,phi_hFIX,check,ome,t,Tetha,x0,y0)

global r x0 y0



Vh=1;
Vy=Vh*cos(Tetha/57.29);
Vx= Vh*sin(Tetha/57.29);
v_c=r*cos(Tetha/57.29);

if Tetha<=180
    y_hc=y0+sqrt(r.^2-(v_c^2));
else
    y_hc=y0-sqrt(r.^2-(v_c^2));
end

if (check==1)
    phi_hc=ome*t*57.29;
else
    phi_hc=phi_hFIX;
end
x_hc=x0-r*cos((Tetha./57.29));