function [c,s] = sincos(theta,flag)
format long
if flag == 1,
   c=cos(theta);
   s=sin(theta);
else
   a=tan(theta/2);
   aa=a*a;
   a1=1+aa;
   s=2*a/a1;    %   sine
   c=(1-aa)/a1; % cosine
end
