	real*8 pi,d
	integer*8 i
	equivalence (pi,i)
	pi=0
	do k=1,4001,4
	d=1.d0/k-1.d0/(k+2)
	pi=pi+d
	print*,4*real(pi),real(d)
	enddo
	i=i+2_8**53
	print*,'end=',pi
	end
	