function [y] = funcion_heading(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, fin el punto mas cercano de la linea que
%contiene el eje central de la pista y que es perpendicular a esta
p=[0 abs(u1)];

%hd: vector resultado de la suma ponderada de r y p
hd=c1*r+c2*p;

%calculo del angulo que forma el vector suma con el vector perpendicular
cosang=dot(hd,p)/(norm(hd)*norm(p));
ang=(acos(cosang))*180/pi;

%Salida
y(1)=-ang;