function [ra, dec] = ra_and_dec_from_r(r)
% –––––––––––––––––––––––––––––––––––––––––––––---------------------------
% Esta funcion 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: ninguna
% Sub-función requerida: ninguna
% –––––––––––––––––––––––––––––––––––––––––––––---------------------------
l = r(1) / norm(r);
m = r(2) / norm(r);
n = r(3) / norm(r);
dec = asind(n);
if m > 0
    ra = acosd(l / cosd(dec));
else
    ra = 360 - acosd(l / cosd(dec));
end
end % ra_and_dec_from_r