🚀 go-pugleaf

In article <1708@sousa.ltn.dec.com> perkins@sousa.ltn.dec.com (Eric Perkins) 
writes:
>There was an article posted a while back asking for an infinite 
>PI program. Did anyone have such a program ?

Here goes:

	program pi
        implicit double precision (a-h,o-z)
c
        a=1.d0
        b=1.d0/dsqrt(2.d0)
        t=.25d0
        x=1.d0
c
        do 10 i=1,10
         y=a
         a=.5d0*(a+b)
         b=dsqrt(b*y)
         t=t-x*(y-a)*(y-a)
         x=2.d0*x
         p0=p1
         p1=.25d0*(a+b)*(a+b)/t
         write(6,*)i,p1
 10     continue
c
        stop
        end