double precision s,c,n,h,j,sNew,cNew

	n=1 000 000 000
	h=1.0d0/n
	s=0
	c=1

		do 17 j=1,n
		sNew=s+c*h
		cNew=c-s*h
		s=sNew
		c=cNew
17		continue

	print *,s,c
	end