[rec,dec,r]=moonequ(T)

function [rec,dec,r]=moonequ(T)

% function [rec,dec,r]=moonequ(T)
%
% equatorial coordinates of the moon                                     
% algorithms by Montenbruck/Pfleger: "Astronomy with the Personal Computer"
%
% 05.05.04 M.Penzkofer

% ecliptical coordinates
[l,b,r] = moon(T);
% mean equinox of the date
[x,y,z] = pol2cart(r,b,l);
% transformation in equatorial coordinates
[x,y,z] = ecl2equ(T,x,y,z);
% nutation
[x,y,z] = nutequ(T,x,y,z);
% equatorial coordinates
[r,dec,rec]= cart2pol(x,y,z);