function [y] = funcion_headingok(u1,u2)
%u1: coordenada x del avion
%u2: coordenada y del avion
%u3: coordenada z del avion
%c1,c2;constantes multiplicadoras
c1=0.01;
c2=0.99;
%r: vector con origen el avion y fin el origen de pista, sin coordenada z
r=[u2 abs(u1)];
%p: vector con origen el avion y fin el punto mas cercano de la linea que
%contiene el eje central de la pista
p=[0 abs(u1)];
%hd: V
hd=c1*r+c2*p;
%calculo el angulo que forma el vector suma con el vector perpendicular
cosang=dot(hd,p)/(norm(hd)*norm(p));
ang=(acos(cosang))*180/pi;
%if u1>0
% heading=360-ang;
%else
% heading=180+ang;
%end
y(1)=-ang;
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
