[xt,yt,zt]=ecl2equ(T,x,y,z)

function [xt,yt,zt]=ecl2equ(T,x,y,z)

% function [xt,yt,zt]=ecl2equ(T,x,y,z)
%
% transformation from ecliptical coordinates to equatorial coordinates
% algorithms by Montenbruck/Pfleger: "Astronomy with the Personal Computer"
%
% 05.05.04 M.Penzkofer

EPS = 23.43929111-(46.8150+(0.00059-0.001813*T)*T)*T/3600;
C = cos(EPS*pi/180);
S = sin(EPS*pi/180);
xt = x;
yt = C*y-S*z;
zt = S*y+C*z;