[DLAM,DS,GAM1C,SINPI]=addsol(DLAM,DS,GAM1C,SINPI,COEFFL,COEFFS,COEFFG,COEFFP,P,Q,R,S,CO,SI)

function [DLAM,DS,GAM1C,SINPI]=addsol(DLAM,DS,GAM1C,SINPI,COEFFL,COEFFS,COEFFG,COEFFP,P,Q,R,S,CO,SI)

% function [DLAM,DS,GAM1C,SINPI]=addsol(DLAM,DS,GAM1C,SINPI,COEFFL,COEFFS,COEFFG,COEFFP,P,Q,R,S,CO,SI)
%
% subroutine of function moon
% algorithms by Montenbruck/Pfleger: "Astronomy with the Personal Computer"
%
% 05.05.04 M.Penzkofer

[x,y] = moonterm(P,Q,R,S,CO,SI);
DLAM  = DLAM +COEFFL*y;
DS    = DS   +COEFFS*y;
GAM1C = GAM1C+COEFFG*x;
SINPI = SINPI+COEFFP*x;