clsin.m

function  result = clsin(ar,degree,argument);
%Clenshaw summation of sinus of argument

% Written by Kai Borre
% December 20, 1995

% See also WGS2UTM

cos_arg = 2*cos(argument);
hr1 = 0;
hr = 0;
for t = degree:-1:1
   hr2 = hr1;
   hr1 = hr;
   hr = ar(t)+cos_arg*hr1-hr2;
end;
result = hr*sin(argument);
%%%%%%%%%%%%%%%%%%%%%%% end clsin.m  %%%%%%%%%%%%%%%%%%%%%