%% ECI
% –––––––––––––––––––––––––––––––––––––––––––––---------------------------
% UTN -FRH
% MECANICA ORBITAL
% SISTEMAS DE COORDENADAS Y MATRIZ DE ROTACION
% –––––––––––––––––––––––––––––––––––––––––––––---------------------------
% Esta programa calcula la ascencion recta y la declinacion a partir del
% vector posicion en sistema ecuatorial geocentrico
%
% Variables utilizadas: 
% r       - Vector posicion
% l, m, n - direccion de los cosenos de r
% ra      - ascension recta (grados)
% dec     - declinacion (grados)
%
% Función M del usuario requerida: ra_and_dec_from_r
% Sub-función requerida: ninguna
% –––––––––––––––––––––––––––––––––––––––––––––---------------------------
%% Limpieza
close all   % Cierra todas las ventanas abiertas
clear       % Borra todas las variables guardadas en el Workspace
clc         % Limpia el Command Window

%% Datos
r = [-5368 -1784 3691];

%% Desarrollo
[ra dec] = ra_and_dec_from_r(r);

%% Resultado
fprintf('- - - - - - - - - - - - - - - - - - - - \n')
fprintf('\nAscension recta y declinacion\n')
fprintf('\n r = [%g %g %g] (km)', r(1), r(2), r(3))
fprintf('\n ascension recta = %g deg', ra)
fprintf('\n declinacion = %g deg', dec)
fprintf('\n- - - - - - - - - - - - - - - - - - - - \n')