Cálculo del ángulo de Azimut y Elevación en MATLAB

%Cálculo del ángulo de elevación y azimut
%
%Datos del Satélite Geoestacionario
Los=30; %Longitud en Grados sexagesimales (ejm: Hispasat 1D   30.0 W )
%
%Datos de la posición de la Antena Parabólica Ricardo Palma
%Lop=76.65; %Longitug de la ant.en Grados sexagesimales
%Lap=11.95; %Latitud de la ant. en Grados sexagesimales
%
%Datos de la posición de la Antena Parabólica
Lop=75; %Longitug de la ant.en Grados sexagesimales
Lap=11.85; %Latitud de la ant. en Grados sexagesimales
%
R=6378; %Radio de la tierra en el Ecuador en Km
h=35768;%Altura del satélite con respecto a la superficie de la Tierra en Km
%
%Cálculo de angulos en radianes
Loss=Los*pi/180;
Lopp=Lop*pi/180;
Lapp=Lap*pi/180;
%
%Cálculo de x
DLo=Lopp-Loss;
x=acos(cos(DLo)*cos(Lapp));
%
%Cálculo de Azimut
v=(h+R)*cos(x);
r=(h+R)*sin(x);
m=v*tan(Lapp);
y=v*sec(Lapp);
ley1=(R+h)*(R+h)+y*y-2*(R+h)*y*cos(DLo);
s=sqrt(ley1);
ley2=(r*r+m*m-s*s)/2/r/m;
az=acos(ley2);
azimut=az*180/pi
if (Los<Lop)
        disp('          Este');
    elseif (Los>Lop)
        disp('          Oeste');
else
      disp('            Ubicación Perpendicular')
end
%
%Cálculo de Ángulo de Elevación
mm=(v-R)/r;
elev=atan(mm);
elevacion=elev*180/pi
%
%
%NOTA: Los textos y símbolos en color verde son comentarios
%   que pueden ser omitidos.