[X,Y]=moonterm(P,Q,R,S,CO,SI)

function [X,Y]=moonterm(P,Q,R,S,CO,SI)

% function [X,Y]=moonterm(P,Q,R,S,CO,SI)
%
% subroutine of function moon
% algorithms by Montenbruck/Pfleger: "Astronomy with the Personal Computer"
%
% 05.05.04 M.Penzkofer

I(1)=10+P; I(2)=10+Q; I(3)=10+R; I(4)=10+S;
X=1.0; Y=0.0;

for K = 1:4
   if (I(K) ~= 10)
      [X,Y]=addtheo(X,Y,CO(I(K),K),SI(I(K),K));
   end
end