[vx,vy,vz]=aberrat(T)

function [vx,vy,vz]=aberrat(T)

% function [vx,vy,vz]=aberrat(T)
%
% stellar aberration - velocity vector of the Eearh
% (equatorial, in units of the speed of light)
% algorithms by Montenbruck/Pfleger: "Astronomy with the Personal Computer"
%
% 05.05.04 M.Penzkofer

% constants
P2 = 2*pi;

L = P2*frac(0.27908+100.00214*T);
CL= cos(L); 
vx = -0.994E-4*sin(L);
vy = +0.912E-4*CL;
vz = +0.395E-4*CL;